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PREFACE 


The National Aeronautics and Space Administration through 
its Tropospheric Chemistry Program has sought and utilized the 
advice of the scientific community in the formulation of its re- 
search efforts. The NASA Working Group on Tropospheric Program 
Planning, led by Dr. John Seinfeld, established a foundation to 
guide program development (NASA RP-1062, 1981). One of four 
areas recommended for expanded activities was instrument develop- 
ment aimed at improving our capability to measure the important 
trace gases and aerosols which are the key species in the major 
atmospheric biogeochemical cycles. To assist in guiding our 
utilization and development of instrumentation and to focus more 
clearly on specific needs, the Instrument Workshop for HxOy Tro- 
pospheric Species was conducted in Palo Alto, California, August 
16-20, 1982. The workshop participants were asked to provide 
assessments of existing sensors and of new technologies deemed 
to have potential for measuring Hj^Oy species at levels character- 
istic of the nonurban troposphere. 

The availability of adeguate instrumentation is a preregui- 
site for the implementation of a program to address the pertinent 
scientific issues of trace constituents within the global tropo- 
sphere and the factors that control their concentrations. The 
efforts of this workshop provide meaningful guidance for future 
instrument developments. The assistance of all who contributed 
to this workshop is sincerely appreciated. Special thanks are 
extended to James M. Hoell, Jr, Workshop Chairman; Malcolm J. 
Campbell, Coordinator, and Warren D. Hypes and Richard J. Bendura, 
NASA Rapporteurs. The support of t)ie Bionetics Corporation, and 
particularly the assistance of Howard J. Curfman, Jr., and Helen 
Ann Thompson, is also gratefully acknowledged. 


Robert J. McNeal, Manager 

Tropospheric Chemistry Program 

Office of Space Science and Applications 
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EXECUTIVE SUMMARY 


The National Aeronautics and Space Administration through its 
Tropospheric Chemistry Program has sought and utilized the advice 
of the scientific community in the formulation of its research ef- 
forts. The NASA Working Group on Tropospheric Program Planning 
established a foundation to guide program development (Seinfeld et 
al., 1981). One of the four areas recommended for expanded activi- 
ties was instrument development, which is aimed at improving our 
capability to measure the trace gases and aerosols which are the 
key species in the major atmospheric biogeochemical cycles. To 
focus more clearly on specific needs, the Instrumentation Workshop 
for H^Oy Tropospheric Species was conducted in Palo Alto, California, 
August 16-20, 1982. The objectives of this workshop were: 

1. To provide an assessment of the capability of the existing 
sensors of H^Oy species at levels characteristic of the 
nonurban troposphere 

2. To identify those techniques adaptable to real-time mea- 
surement (sampled or continuous) onboard an aircraft plat- 
form 

3. To address the concern for intercalibration and intercom- 
parison of techniques for measuring specific species 

4. To recommend promising technologies for research and de- 
velopment for measuring important species for which a 
capability does not now exist or is presently of inade- 
quate sensitivity 

Workshop attendees included approximately 16 invited participants 
from government, universities, and industry. Each participant was 
selected based on a demonstrated expertise in either the develop- 
ment or application of HxOy monitoring techniques for global trop- 
ospheric research. 

As input into the workshop, NASA prepared an overview of its 
current Global Tropospheric Program, highlighting ongoing activi- 
ties as well as long-range plans. Included in these discussions 
was a brief summary of HxOy species chemistry in the troposphere 
and required measurement capabilities for various HxOy species. 

These measurement capabilities (horizontal and vertical spatial 
resolution and required concentration levels), combined with the 
experiences of the workshop participants, became the basis for 
discussion of the various techniques. This workshop report is a 
summary of those discussions and written material prepared by the 
participants to describe the techniques discussed. 

The species discussed in this workshop document include OH, 

H 202 f HO 2 and CH 3 OOH. The techniques encompassed direct measure- 
ment concepts utilizing various approaches to the laser-induced 



fluorescence technique, absorption spectroscopy, wet chemical 
techniques and several indirect approaches which depend upon mea- 
surements of reactant products or rate of removal of global tra- 
cers. The major emphasis at the workshop, however, centered on 
an assessment of techniques for the measurement of OH. That this 
was so recognizes the fact that OH is the single most important 
constituent in the troposphere and that major differences current- 
ly exist in the understanding of its global distribution. More- 
over, efforts to develop techniques suitable for measuring tropo- 
spheric OH span more than a decade and have included university, 
private industry, and government laboratories. The few OH mea- 
surements that have been reported by any approach are, at best, 
questionable. Further, it was recognized that of the various ap- 
proaches to measure OH, the ones that have received the most at- 
tention and indeed generated the most controversy are those that 
employ the LIF concept. Accordingly, a large fraction of the 
workshop's discussion was in fact devoted to an assessment of 
existing and proposed LIF sensors for measuring OH. The workshop 
focused on assessment of the more prominent LIF concept, namely, 
a single photon excitation process. Techniques utilizing multi- 
ple photon excitation steps were not evaluated and hence will not 
be discussed in this document. The assessment activity included 
(1) developing a generalized or "central" LIF equation from which 
the working expression (e.g., relating OH concentration to observ- 
ed signal) for any LIF approach could be derived, (2) developing 
a common formulization for the interference due to laser-generated 
OH, and (3) using the results from (1) and (2) to compare the ex- 
pected or achieved performance of various LIF approaches. 

The formulation of a "central" LIF equation was, in fact, 
initiated at an ad hoc meeting of OH investigators convened at 
the 2nd Symposium on the Composition of the Nonurban Troposphere, 
held in Williamsburg, VA, May 1982. The rationale for such an 
effort stems from the fact that all LIF approaches for OH measure- 
ments are based upon a common physical phenomenon which can be 
mathematically formulated to explicitly relate the OH concentra- 
tion to an observed signal for a given laser flux. Such a for- 
malization can be cast in terms of standard spectroscopic quanti- 
ties and can be general enough to encompass the various approach- 
es for implementation of the LIF technique. The working expres- 
sion for each LIF instrument would, in effect, be a specialized 
case of the central equation. The advantage of such an approach 
is that comparisons of current and future LIF systems can be 
implemented from a common basis. 

The OH measurement techniques discussed in the document are 
at various stages of development, ranging from concepts (on paper 
only) to operational in-the-field instruments. For techniques 
existing as a concept or as laboratory instruments, discussion fo- 
cused on anticipated measurement capabilities and future research 
and/or hardware development required to obtain anticipated capabi- 
lities. For operational techniques, discussion focused on current 
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capabilities, known instrument interferences and shortcomings, and 
research required to improve instrument cabilities. A brief de- 
scription of each technique discussed at the workshop is included 
in the document. 

Conclusions, observations, and the consensus of opinion aris- 
ing from the workshop are briefly noted in the following summary 
material . 

1 . Alternate techniques for OH measurement should continue 
to be explored and supported if found to offer signifi- 
cant advantages. These techniques include; 

a) LIF measurement with sample expansion to low pres- 
sures 

b) LIF measurement using higher repetition rates 

c) LIF measurement using narrower band detectors 

d) LIF measurement using 308-nm excitation in the 
0-0 band 

e) LIF measurements using short laser pulse widths 

f) Measurement techniques analogous to the radiochemical 
method but using other, possibly nonradioactive, 
reagent species 

2. A better understanding, at least at the empirical level, 
is needed of the dynamics of the OH radical, in order 

to assess the improvements to be expected in reducing 
interference in the LIF method from ozone photolysis. 

3. NASA should encourage the assembly of carefully reviewed 
compilations of data, such as quenching constants for OH* 
and H0*2, of interest to experimentalists. 

4. There should be coordinated studies of the relationship 
of the nonresonant fluorescent background to LIF measure- 
ments and the results of such studies should be dissemi- 
nated. Features of interest include meteorology, on 
which some data may already be available, spectrum, and 
two-point autocorrelation function. 

5. Development of OH calibration techniques should be expe- 
dited, including exploration of common calibration sys- 
tems such as that proposed by the Ford group and improve- 
ments in the open-air long-path absorption technique. 

6. Intercomparison of OH instruments at an early time should 
be encouraged. NASA should endeavor to see that the best 
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calibration feasible is utilized in each field expedition, 
and both the relative calibration factors and calibration 
stabilities of all instruments should be determined and 
treated as part of the record of the field expedition. 

7. Through this report and other means, NASA should request 
the community to consider possible techniques for air- 
craft measurement of HO 2 , and CH^OOH. 

INTRODUCTION 

In its continuing efforts to direct its applications programs 
toward relevant national needs, NASA is conducting the Tropospher- 
ic Chemistry Program. The long-range objective of this program 
is to apply NASA's space technology to assess and predict human 
impact on the troposphere, particularly on the regional to global 
scale. The increasing importance of pollution on these scales 
and the synoptic view afforded from satellites suggest that space 
observations can play a unique and critical role in satisfying the 
objective. A NASA-sponsored working group of scientists prepared 
an overview of the scientific problems that need to be addressed 
in order to understand the large-scale troposphere (Seinfeld et 
al., 1981). The group recommended that NASA undertake expanded 
efforts to develop space applications for tropospheric air quality 
monitoring in the areas of instrument development, modeling, lab- 
oratory studies, and field measurement activities. 

To help NASA in formulating a detailed implementation plan, 
assistance has been sought from workers in this field to identify 
the high-priority research required in each of the four components. 
Such studies included specific regional scientific concerns (Levine 
and Schryer, 1978); the results from and the needs of tropospheric 
modeling (National Aeronautics and Space Administration, 1981); 
multiphase processes, including heterogeneous catalysis (Schryer, 
1982); passive remote sensing (Keafer, 1982); and instrumentation 
for NxOy tropospheric species (National Aeronautics and Space 
Administration, 1983). 

This volume reports the results of the Instrumentation Work- 
shop for HxOy Tropospheric Species. The workshop was sponsored by 
the NASA Office of Space Science and Applications and conducted un- 
der the direction of the Global Tropospheric Experiment Project Of- 
fice of NASA Langley Research Center. 

OBJECTIVES AND ORGANIZATION 

The objectives of this workshop were; 

(1) To provide an assessment of the capability of existing 
sensors of H^Oy species at levels characteristic of the 
nonurban troposphere 
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(2) To identify those techniques adaptable to real-time mea- 
surement (sampled or continuous) onboard an aircraft 
platform 

(3) To address the concern for intercalibration and inter- 
comparison of techniques for measuring specific species 

(4) To recommend promising technologies for research and 
development for measuring important species for which a 
capability does not now exist or is presently of inade- 
quate sensitivity 

All workshop participants were actively involved in the de- 
velopment of instruments or the utilization of promising tech- 
niques for the measurement of gaseous H^Oy species. The planning 
for the workshop and the overall leadership was the responsibility 
of the workshop chairman. A coordinator was selected from the 
participants to conduct the workshop, make writing assignments, 
and assist the chairman in this final documentation. Two NASA 
rapporteurs also provided assistance in coordinating the documen- 
tation assignments. Those contributing to the discussions and 
documentation are as follows: 


James M. Hoell, Jr., Workshop Chairman 
Warren D. Hypes, NASA Rapporteur 
Richard J. Bendura, NASA Rapporteur 
Malcolm J. Campbell, Coordinator 


PARTICIPANTS 


James G. Anderson 
John Bradshaw 
David R. Crosley 
Douglas D. Davis 
L. I. Davis 
Carl Farmer 
Thomas M. Hard 


William S. Heaps 
I. Stuart McDermid 
Matthias Rateike 
Michael O. Rodgers 
Hanwant B. Singh 
Richard M. Stimpfle 
Charles C. Wang 


MEASUREMENT REQUIREMENTS FOR HxOy TROPOSPHERIC SPECIES 


Highly reactive free-radical molecules are recognized as cri- 
tical to maintaining atmospheric homeostasis. Because free radi- 
cals have an unpaired electron in their outer shell they can act 
as strong oxidizers of reduced gases emitted to the atmosphere 
from the biosphere. The H^Oy species, and in particular the OH 
radical, have been identified as pivotal in tropospheric photo- 
chemistry (Levy, 1971). A large literature documents hypothesiz- 
ed photochemical reaction pathways for H^Oy species in the tropo- 
sphere (e.g., Wofsy et al., 1972; Crutzen, 1973; Seinfeld et al., 
1981; National Aeronautics and Space Administration, 1981), and 
the detailed scientific rationale for the importance of these 
species to understanding atmospheric chemical processes need not 
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be repeated here. The critical fact to this workshop is that few, 
if any, widely accepted measurements exist for HxOy species in non- 
urban air. Mathematical models have been used to calculate con- 
centrations of species such as OH. Such calculations suggest an 
average value of 7 ± 3 x 10^ cm~3 in the remote troposphere, with 
highest levels in the tropics. Direct measurement of gas phase 
OH, HO 2 , and H 2 O 2 is a top priority for further progress in under- 
standing the chemistry of the troposphere. 

In contrast to OH and HO 2 , the most interesting aspects of 
the chemistry of H 2 O 2 may occur in the liquid phase, especially 
in cloud droplets. Chameides and Davis (1982) have proposed a 
mechanism by which H 20 y radicals generate H 2 O 2 in solution, with 
the H 2 O 2 subseguently becoming involved in heterogeneous oxida- 
tion of species like SO 2 . It is important to pursue the develop- 
ment of direct measurement techniques for H 2 O 2 in both gas and 
liquid phases. 

A summary of requirements identified in existing literature 
for instrument sensitivity, spatial resolution, and estimated 
residence times for OH, HO 2 , and H 2 O 2 are given in table 1. 

TABLE 1. REQUIREMENTS FOR H^Oy MEASUREMENTS 


Residence Spatial Resolution, km Detection 


Species 

Time 



Vertical 

Horizontal 

Requirement 

OH <1 

sec - 

1 

day 

0.5 

200 

5 X 1 0 5 cm ^ 

HO 2 <1 

sec - 

1 

day 

0.5 

200 

1 X lo"^ m"^ 

H 2 O 2 <1 

day - 

1 

week 

0.5 

200 

200 pptv 


SPECIES AND TECHNIQUES CONSIDERED 

Table 2 lists the techniques considered by the workshop par- 
ticipants along with the respective species. Most of the workshop 
discussion centered on techniques for detection of OH. This re- 
flects the uncertain state of OH measurements and the long history 
associated with efforts to develop technigues to measure tropo- 
spheric OH coupled with the scarcity of techniques for measuring 
the other H^Oy species. Among the OH technigues listed in table 2, 
the laser-induced fluorescence (LIF) method received the most con- 
sideration. The workshop discussion addressed only LIF approaches 
utilizing a single photon excitation step. 


6 




The major emphasis was placed on a detailed assessment of the 
status of LIF systems that have undergone some level of either 
field or laboratory testing. The approach taken for the assess- 
ment focused on the development of a generalized central eguation 
to relate OH concentration to observed signal, laser flux and stan- 
dard spectroscopic guantities. The working expression used by each 
LIF system to relate the OH concentration to observed signal was 
obtained from the central eguation as a specialized case. A major 
concern for implementation of the LIF technigue for OH measurements 
is the potential interference effects produced by laser-generated 
OH. The various approaches for implementing the LIF concept are, 
in fact, fostered in large measure by a need to minimize the ser- 
iousness of this interference. Accordingly, it was felt by all 
participants that an assessment of the OH LIF systems reguired an 
evaluation of the magnitude of the laser-generated OH as well as 
a common formulation of this phenomenon for several of the current 
and proposed techniques. The actual status assessment culminated 
in a comparison of the signal-to-noise ratio expected from several 
LIF systems "operating" in two different model atmospheres. 

In the following sections, a brief description of the various 
approaches to measure OH is given. The next section will summarize 
the workshop evaluation of the\ LIF concepts discussed. 


TABLE 2 . TECHNIQUES CONSIDERED 


Detection Principle 


Technique 

Species 

Laser-induced 

fluorescence 

(1) 

Single Photon 

a) Atmospheric pressure 
in situ 

b) Low pressure in situ 

c) Lidar 

OH 


(2) 

Photo fragment at ion 

H2O2 ,CH300H 

Radiochemical 

(1) 

l^CO/OH reaction 

OH 

Global Tracers 

(1) 

Emission vs. Removal 
via OH reaction 

OH 

Absorption 

Spectroscopy 

(1) 

Long-path UV 

OH 

Wet Chemical 

(1) 

Condensation sampling 

H 2 O 2 
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MEASUREMENT OF OH 


Laser-Induced Fluorescence 

Laser-induced fluorescence (LIF) is a measurement technique 
that, under laboratory conditions, has demonstrated measurement 
sensitivities suitable for detection of single atoms and molecules. 
This demonstrated sensitivity, coupled with the fact that LIF re- 
resents one of the few techniques for direct detection of OH, has 
been the impetus for developing this technique for ambient OH mea- 
surements . 

It should be recognized at the outset that there has been 
and currently exists a wide diversity in the research and develop- 
ment efforts directed toward developing field-compatible LIF sys- 
tems. Because of this diversity, it is easy to lose sight of the 
fact that all the LIF approaches are based upon a single funda- 
mental phenomenon. Simply stated, this phenomenon involves reso- 
nant excitation of a molecule either with a single photon at 
or with two sequential photons at X^ and then X 2 , followed by 
detection of the resulting fluorescence. In general, the single- 
photon excitation scheme discussed at the workshop has typically 
been initiated using excitation wavelengths of either 282 nm or 
308 nm. The 282-nm excitation scheme, illustrated in figure 1„ is 
the only one which to date has been employed for atmospheric mea- 
surements. As shown in figure 1, the OH radical is initially ex- 
cited into the v' = 1 manifold of the A^ e electronic state by 

R-R COLLISIONAL 
RELAXATION 



R-R COLLISIONAL 
RELAXATION 


Figure 1. Laser-induced fluorescence, 282-nm excitation. 
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pulsed laser excitation in the region of 281.9 to 282.5 nm. This 
corresponds to either ( 1 ) , 0i ( 2 ) , P][(l) or Pi(2) transitions and 
will be referred to here as 282-nm excitation. Rotational and vi- 
brational relaxation of this excited state causes a significant 
part of the OH fluoresence to occur from the (A^Z, v' = 0) ^ 

(X^n, v' = 0) ) transition around X = 309.5 nm. However, due to 
collisions with N 2 r 02 / and H 2 O, significant electronic quenching 
occurs resulting in attenuated fluorescence from rotational states 
in the v' = 1 and v' = 0 manifold of the E state. The rela- 
tively large red shift in the fluorescence signal (e.g., 309 nm) 
relative to the pump wavelength (e.g., 282 nm) is one of the advan- 
tages of this particular excitation scheme. In particular, the OH 
fluorescence can be detected with minimal interference from the 
pumping freguency. The disadvantage relative to the 308-nm excita- 
tion scheme discussed below is the increase in photolytic interfer- 
ence (e.g., laser-generated OH, discussed below). For the scheme 
shown in figure 1, the fluorescence emission centered at 309 nm is 
typically measured using a narrow bandpass filter backed by a high- 
efficiency photomultiplier tube and photon-counting electronics. 

The LIF scheme just outlined is deceptively simple, particu- 
larly when considered for field measurements of a trace atmospheric 
constituent such as OH. Successful implementation of this tech- 
nique requires development of state-of-the-art laser excitation and 
detection systems that optimize the return signal while minimizing 
the contribution from a number of noise sources. It is important 
to note that the various LIF approaches currently being consid- 
ered for OH measurements differ primarily in the rationale and 
procedures used to minimize the noise contribution. The major 
noise sources that must be overcome are: ( 1 ) nonresonant fluores- 

cence (NRF) from other atmospheric constituents; (2) Rayleigh, 

Raman and Mie scattering; (3) photolytic interference, resulting 
from photolysis of other molecules to produce OH; and (4) external 
background radiation such as solar background. The Rayleigh, 

Raman and Mie scattering, as well as contributions from external 
background radiation, represent well-understood and characterized 
phenomena which, while not insignificant, can be minimized with 
straightforward design precautions and use of narrowband filters. 

Nonresonant fluorescence is generally attributed to broadband 
fluorescence from aerosols. This noise source has not been fully 
characterized and several instrumental approaches involving combi- 
nations of spectral, temporal and baric filtering are under consid- 
eration for minimizing its effect. To date photolytic interference 
has been the major source of uncertainty in the measurement of 
tropospheric levels of OH via the LIF approach. This interference 
originates from the production of OH via photolysis of O 3 by laser 
radiation producing 0 (^ 0 ) followed by its subsequent reaction with 
ambient H 2 O. Like the LIF phenomenon itself, photolytic interference 
is a fundamental process inherent in the LIF measurement approach 
and represents a major noise source for which a number of instru- 
ment approaches have been employed to minimize its consequences. 
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An excitation and detection scheme which has been suggested as 
a method for minimizing the photolytic interference is illustrated 
in figure 2. Conceptually this is similiar to the 282-nm excita- 
tion scheme. Here, however, the excitation is with X = 30.8 nm. 

The OH fluorescence originates from the same manifold, resulting 
in a shift of only a few tenths of a nanometer. Thus, the prim.e 
disadvantage of this excitation scheme is the need for a high 
resolution filter to discriminate between signal and the back- 
ground originated from the excitation pulse. An important advan- 
tage is that 0^ absorption at 308 nm is significantly less than at 
282 nm, resulting in less photodissociation of O 3 and therefore 
lower OH generated from the excitation laser. A more detailed dis- 
cussion of this excitation/detection scheme will be given below. 



Figure 2. Laser-induced fluorescence, 308-nm excitation. 


There are a number of LIF approaches in various stages of de- 
velopment. In general, the distinguishing feature between those 
approaches is the method employed to minimize contributions from 
the noise sources noted above while maximizing the SNR. Table 3 
identifies the major LIF techniques currently under consideration, 
along with the more obvious distinguishing features. The status 
column provides a general assessment of the field readiness of the 
approach while the comment column indicates the major advantage or 
objective of the given approach. The source column identifies one 
or more of the laboratories that support a given technique. Each 
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TABLE 3. SALIENT CHARACTERISTICS OF LASER-INDUCED FLUORESCENCE (LIF) TECHNIQUES FOR 
MEASUREMENT OF TROPOSPHERIC OH 


TECHNIQUE 

^on 

(nm) 

^det 

(nm) 

STATUS 

COMMENTS 

SOURCE 

(1) 

(2) 

Lidar 
In Situ 

282 

309 

Field Tested 

Large sanple volume 
to minimize laser - 
generated OH and 
absorption saturation 

Ford, 

Goddard 


(a) Atmospheric 
Pressure 

282 

309 

Field Tested 

Short pulse to minimize 
laser-generated OH; 
near-simultaneous Xon 
& Xoff to freeze non- 
resonant fluorescence 

Ga. Tech 



282 

309 

Proposed 

Low pulse energy, 
high repetition rate 

Harvard 



308 

309 

Proposed 

Excimer laser with 

high-resolution 

filtering 

JPL 


(b) Low Pressure 
(Baric 
filtering) 

308 

309 

Under develop- 
ment 

Minimize laser-generated 
& background noise 

Portland 
State/I nstitut 
fCfr Chemie 


of these approaches was discussed at the workshop. A brief de- 
scription of each technique will follow the discussion of the 
dominant noise sources. 

Interferences . Of the four major sources of noise or interference 
noted in the previous section, two are considered well understood. 
These are Rayleigh, Raman, and Mie scattering and external radia- 
tion, particularly the solar background. These phenomena are 
well characterized and, while not insignificant, can be minimized 
through various design considerations and filtering techniques. 
This leaves the nonresonant fluorescence and photolytic interfer- 
ence as the two remaining major sources of noise. To appreciate 
the impact of these two noise sources, one must realize that the 
raw measured signal in an LIF system is, in general, the sum of 
the desired atmospheric OH fluorescence plus contributions due to 
fluorescence associated with the photolyt ically generated OH and 
the NRF. The NRF is generally a slowly varying function of both 
the excitation and fluorescence spectrum. Accordingly, correction 
for the NRF contribution can be implemented through simultaneous 
measurements in a spectral channel close to the OH fluorescence 
channel or by temporally separated measurements in the single 
OH fluorescence channel with the NRF excited by wavelength Xq^ 
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and then ^off* account for the contributions from photoly- 

tically generated OH, the sampling geometry and operating para- 
meters are designed to minimize the amount of OH created. In 
addition, the signal contribution due to the laser-generated OH 
can be estimated and mathematical corrections implemented. 

Nonresonant Fluorescence . The nonresonant fluorescence phenomenon 
concerns a large and often variable nonresonant, broadband fluor- 
escence background signal. In some environments, it has been ob- 
served that the uncertainty in this background signal does not 
follow simple Poisson statistics (Wang and Davis, 1982; Rodgers et 
al., 1982). These observed fluctuations appear to be related to 
inhomogeneities in the levels of aerosol species. The net effect 
of this nonstatistical noise source is that the SNR in ,a given 
system does not increase with the square root of the integration 
time. Hence, a significant deterioration occurs between the cal- 
culated or expected detection limit and that actually obtained 
during field measurements. In principle, the above problem could 
be solved if the switching time for tuning the laser "on" and 
"off" the OH resonance frequency is made short compared to the 
rate at which changes occur in the environment. It is known for 
example that at time intervals as short as 1 ms, atmospheric mo- 
tion is unimportant. What is not yet known is whether the "on"/ 
"off" cycling time for an LIF system needs to be this fast. 

It is clear, from the measurements which have been made, that 
NRF can dominate the noise in all but the lidar-type instruments 
with unterminated beams. The noise in unterminated lidar systems 
is generally dominated by solar scattering. It is less clear, but 
nevertheless true, that the flux of air through the sample volume 
interrogated by the instrument (a parameter partly determined, 
for example, by the speed of the aircraft on which the instrument 
is mounted) has a major influence on the impact of NRF. At this 
point our knowledge of the causes and properties of the NRF is 
very limited. A crucial characteristic of NRF is the autocorre- 
lation function describing the NRF signal, both in time and space. 
Optimum instrument design to minimize NRF requires that signal 
and background measurements be spaced temporally (or spectrally) 
close enough so that the NRF signal on the two occasions is essen- 
tially unchanged. Unpublished data from the Ford group suggest 
that the autocorrelation function (at ground level, in moderate 
winds) is close to unity for times less than 1 second. These 
limited data are encouraging. The universality of this figure 
and the validity of relating correlations in space to correlations 
in time in this context remain to be established. 

It is reasonable to assume that the NRF signal is proportion- 
al to the spectral bandwidth of the detector. Typical values, 
measured in units of equivalent hydroxyl radical fluorescences 
per unit volume, range from lO^ to lO^ cm“3 for a 2.7 nm band- 
width. These values seem to apply to both the boundary layer and 
the mid-troposphere over land, with the lower figure more common 
for the mid-troposphere over oceans. Much lower values occur oc- 
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casionally, but apparently not for a sufficiently large propor- 
tion of the time to permit representative hydroxyl data to be 
collected . 

In summary the impact of the NRF is twofold. Not only does 
it decrease the s ignal-to-noise ratio by raising the background 
count, but it also causes an additional degradation unless the 
signal and background measurements are close enough in time for 
the background to remain constant. 

Photolytic Interference . As noted above, photolytic interference 
has, in general, represented a major source of uncertainty in mea- 
surements of OH with LIF techniques. This interference originates 
from laser-induced production of the OH radical which is then de- 
tected and interpreted as ambient OH. As discussed by Wang et al. 
(1976), Hanabusa et al. (1977), Davis et al. (1981a, b) , Ortgies 
et al. (1980), and Wang and Davis (1982), at wavelengths in the 
region of 282 nm, the laser photolysis of atmospheric O3 results 
in the production of metastable atomic oxygen 0(^D). And, even 
though only a small fraction of this 0(^D) goes on to react with 
atmospheric H2O to produce OH radicals, the artificial OH so pro- 
duced can prove to be quite troublesome. In fact, depending on 
the laser energy density and the laser pulse width, the signal 
produced from laser-generated OH can overwhelm the signal from 
natural atmospheric OH. Potential solutions to the LIF artifi- 
cial OH problem include the use of expanded laser beams (Wang et 
al., 1981), the use of very short laser pulses (Davis et al., 

1981a, b; Rodgers et al., 1982), the use of sampling manifolds 
that operate at reduced pressure (Hard et al., 1979, 1980), the 
use of high-repetition-rate low-peak-power laser excitation 
(Stimpfle and Anderson, 1982), or the use of 308-nm excitation. 

The physical process for photolytic production of OH is understood, 
and perhaps more importantly, the severity of the interference can 
be estimated and in some cases measured for a given experimental 
configuration. The purpose of this section is to describe the phy- 
sical process which gives rise to laser-induced OH interference and 
to present the "working expressions" for the in situ configuration 
employed by Georgia Tech and Harvard Univ. and the lidar configu- 
ration employed by Wayne State Univ. /Ford Motor Co. 

As noted earlier the single-photon excitation schemes can 
employ either 282 nm or 308 nm wavelength. Excitation of atmo- 
spheric OH at 282 nm introduces a potentially more serious compli- 
cation than might be experienced at 308 nm. This stems from the 
larger O3 absorption cross section at 282 nm than at 308 nm. The 
photolysis of O3 at 282 nm produces atomic oxygen in the ■'■D ex- 
cited state which reacts with ambient H2O to yield OH: 


O3 + hv(282 

nm) ->■ 0(^D) + 02 (^A) 

( 1 ) 

0(^D) + M > 

0(3p) + M 

( 2 ) 

0(^D) + H 2 O 

^ 20H* 

(3) 
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OH* + M ^ OH + M 


(4) 


where OH* signifies vibrat ionally and/or rotationally excited OH 
generated as a result of process (3). Interference then results 
when the excited OH generated in this manner is relaxed through 
collisions into the particular rotational level from which exci- 
tation to the state of OH takes place. The key time con- 

stants for the relevant molecular processes as a function of al- 
titude are summarized in figure 3. Note that the 0(^D) lifetime 


throughout the troposphere is less than 10 nsec. The extremely 
rapid R-R relaxation rate of OH further implies that the OH formed 
from reaction (3) in the v=0 level (approximately 70%) appears ro- 
tationally thermalized to the laser pulse, although the nascent 
rotational distribution is "hot". Hydroxyl formed by reaction (3) 
in the v=l and 2 levels is not seen by the laser because the v-v 
transfer rate is insufficiently fast (see figure 3) in the X^n 
state. Thus, within the troposphere a significant fraction of the 
0(^D) can form OH during the laser pulse. The atmospheric pres- 
sure and laser pulse width play critically important roles, since 
they define the time constants of the system. 


Prior to writing down the LIF-lidar expression used by the 
Ford group and the LIF in situ expressions used by the Georgia 
Tech and Harvard groups for estimating the laser-induced OH, two 
different quantities with which the results are cast are defined; 



r, LIFETIME (nsec) 

Figure 3. Key time constants in the ozone photolysis process. 
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= the average laser-induced OH density which one 
would derive by converting an observed laser- 
induced OH count rate to an effective "steady 
state" density via knowledge of the absolute OH 
calibration of the instrument; in general, 
[ 0 H] 1 ^®®^ is never in steady state during the 
course of a laser pulse 

[OH] laser = the total amount of laser-induced OH formed per 
total pulse, allowing for all 0 ( 1 d) formed to decay; 
in general [OH] laser does not exceed 1/2 
[OH] laser 
total 

The format used here for comparison of the three expression 
is to first give the actual expression as cited in the literature 
along with definitions of the specific parameters used in each 
case. The similarities and differences are then noted. 

The expression used by the Ford group to estimate this inter- 
ference for their lidar system is given by 

[OH) laser = F(x)[OHlJ®tal = *o total '5> 


laser 

[OH] total " 


ao (|) E [O3] 


k[H90] 

k[H20] + k„[N2] + ko[02] 


( 6 ) 


where : 

2 

Oq = absorption cross section of ozone, cm /molec. 
(}) = quantum yield of 0(^D) , unitless 

E = integrated energy flux per pulse, J/cm^ 


k = rate constant for generation of OH from the reac- 
tion of metastable oxygen atoms with water, 
cm^/molec/s 

[H 2 O] , [O 3 ], [N 2 ], [O 2 ] = concentration of water vapor, 

ozone, nitrogen, and oxygen, respectively, molec/cm^ 

knfko = rate constants for collisional de-excitation of 
0 ( 1 d) due to nitrogen and oxygen respectively 


A' = the ratio of the laser OH density residing in the 
quantum state pumped by the laser to that expected 
if thermal equilibrium is attained; in principle, 
Aq may take on values, 0 ^ A<-. £ 1 
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F(X) = 1 - f 



= the deviation of from the "steady state 

solution" due to the pulse width being of finite 
width: primarily, it provides the means to predict 
the amount of OC^D) which forms OH during the 
pulse, 0 £ F{ x) £ 1 

X = At 

t = laser pulse length 

A = + k^lOj] 

Ac = At F(x) 

The corresponding expressions used by the Georgia Tech group 
to estimate the OH interference for their in situ sampling config- 
urations are given by: 


laser ^interf Tl 
[OH] laser = [oHl^otal ; 2 

^nat 


exp(-i7x) - 1 


laser 

[OH] 


total 



2k3 [H 2 O] 

k3[H20] + 1<2[N2] + 


( 8 ) 


where E = laser energy per pulse, J 

a = absorption cross section of ozone, cm^/molec 
(|) = quantum yield of 0(^0), unitless 
d = laser beam diameter, cm 
h = Planck's constant, Js 
c = speed of light, cm/s 

X = excitation wavelength (e.g., 282 nm) , nm 


k 3 = rate constant for generation of OH from the reaction 
of metastable oxygen atoms with water, cm^/molec/s 

^2/^*2 = rate constants for collisional de-excitation of 0(^D) 
due to nitrogen and oxygen respectively, cm^/molec/s 



16 



T = laser pulse width, s 
K = k3[H20] + k2[N2] + k'2[02l 

Finterf “ fraction of the that resides in the 

quantum state pumped by the laser 

fnat ~ fraction of ambient OH (in thermal equilibrium) 

that resides in the quantum state pumped by the 
laser 

The corresponding expression used by the Harvard group to 
estimate the OH interference is derived from a computer simulation 
of the problem by numerically solving the differential equations 
that describe the interference: 

[O^] [H 2 O] E 

[OHjlaser = 3.3 x 10"2 f(x) (9) 

[M] NpA 

where [M] = molec cm“2 

E = average power, mW 


Np = pulse repetition frequency, Hz 

A = beam area, cm^ 


f(x) = a polynomial in x which gives the fraction of 

[OH] laser that is actually detected by the laser 
total 

pulse; the separate effects of 


(1) 0(1d) relaxation 

(2) R-R transfer in OH 

(3) v-v transfer in OH 

are contained within the f(x) function 
f(x) = -4.225 + 7.79750X - 5.7498x2 + 1.9579x3 - .30916x4 
+ .01831x3 

X = 3 + log(3 X 10"11[M]t) 

X = FWHM laser pulse width, sec 
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with the present parameterization, f(x) in principle has the 
limits 0 £ f(x). £ 0.91. The valid limits on x and f(x) which can 
be used in this formulation are 1.75 £ x £ 4.35 and .0166 £ f(x) 

£ .651 because of the range of data fitted and the polynomial used. 

Thus 


[OH] laser = 2 x 0.91 x 3.8 x 10"2 
total 


[O 3 ] [H 2 O] 
[M] 


E 

NpA 


( 10 ) 


Table 4 identifies the corresponding notation used in each 
of the expressions given above. The working expressions derived 
by the Ford and Georgia Tech Group are in effect identical. For 
the Harvard Group and the Ford Group the functional dependence of 
f(x) and F(x) on [M] and t(= At) is essentially identical. 

This result is not surprising since in both approaches the time 
dependence of the interference is controlled largely by 0 ( 1 d) re- 
laxation; i.e., the state-to-state calculation of the Harvard 
group employing rapid R-R relaxation is similar in effect to im- 
posing a fixed thermal R distribution on the v = 0 OH quantum 
states, whereas the Ford group's approach treats only 0(1 d) re- 
laxation explicitly. This correspondence is borne out by direct 
comparison of 


f(x) F(X) = 0.7 F(X) 

as illustrated in figure 4. 


TABLE 4. 

COMPARISON OF 

WORKING EXPRESSIONS FOR 

LASER-INDUCED OH 

WAYNE STATE 
FORD MOTOR 

UNIVERSITY/ 

CO 

GEORGIA TECH 

HARVARD UNIVERSITY 



a 

— 

■)> 


■!> 


At 


T 

T 

A 


K 

— 

X 


¥X 

— 

A* 

C 


Finterf 

Fnat 

— 

F(X) 

2 

1 _ 1 /i - 1 - exp(-TrX)-l\ 1 

2 Ttx^ ttx yj 

f (x) 

E 


E 

E 



(^=f) ■ 14 )" 

NpA 

k 


"3 

— 

k 

n 


CM 

— 

k 

0 


to - 

— 
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Figure 4. Comparison of f (x) and F(x) for the Ford and 
Harvard formulations of laser-induced OH. 


Choosing = 0.7 excludes the laser OH formed in the v = 1 

and 2 levels from contributing to the interference. Thus, if 
=0.7 the expressions of the Harvard, Ford, and Georgia Tech groups 
are in substantial agreement. 

The point of contention between the two theories at present 
is the correct value for A^ or the rate of R-R relaxation in the 
ground electronic state of OH which has yet to be measured de- 
finitively in the laboratory. The Harvard group has adopted the 
viewpoint that R-R relaxation is rapid; Ap = 0.7. The Ford group 
favors the viewpoint that R-R relaxation is slow, or A^ =0.15 
(Wang and Davis, 1982). 

Several general comments are in order. First, the OH interfer- 
ence is directly proportional to the concentration of H 2 O and O 3 . 
Consequently, measurements in the boundary layer with high levels 
of H 2 O can be severely impacted. Second, the [OH] given in 
each of the expressions is dependent upon the laser energy flux. 

From this dependence stems the rationale for the lidar approach 
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being pursued by the Ford and Goddard groups and the approach pro- 
posed by the Harvard group. Namely, the lidar configuration lends 
itself to a reduced energy density at the sample volume because of 
the large volume that is probed, thereby reducing E. The Harvard 
group proposes to reduce the energy flux by utilizing low peak 
pulse energies at high repetition rate. The dependence of 

on the laser pulse width is also included in each ex- 
pression via the X parameter, and more subtly in the lifetime of 
0(^D). The dependence of [OH] l^ser on pressure is also through 
the lifetime of 0(^D). The excitation pulse widths and 0(^D) 
lifetime become comparable for excitation pulses of 1 nsec or 
less at atmospheric pressure. For short pulses and/or low pres- 
sures, the buildup and decay of 0(^D) and consequently the OH in- 
terference lag the laser pulse shape (and consequently the OH in- 
terference) thereby minimizing its effects. Operation with short 
pulses is the basis of the approach that the Georgia Tech group 
has taken while operation at low pressures is the approach pre- 
ferred by the groups from the Kernforschungsanlage (KFA) and 
Portland State University. 

As a final comment concerning the working equations given 
above, it is clear that [OH] l^ser made as small as desired 

through manipulation of various system parameters (e.g., energy 
density). In general, however, this is done at the expense of 
real signal. Therefore, a comparison between the various instru- 
mentation approaches to determine which minimizes the OH interfer- 
ence problem is meaningless if overall system performance is not 
also considered. Because of this, we now resist the temptation 
to utilize the above working equations for absolute numerical com- 
parisons between the one lidar and two in situ approaches. This 
comparison is in fact presented along with system performance 
expectations for current and projected systems in a later section. 

Measurement Systems 

Several systems for the implementation of the basic LIF tech- 
nique (see table 3) were discussed at the workshop. These systems 
are described in this section along with pertinent details related 
to the rationale for the approach and its status. 

Lidar . - The basic lidar system concept consists of a laser source, 
transmitting optics, a receiver telescope, a detection system, and 
signal processing electronics. Use of the lidar configuration, in 
general, offers some advantages related to more efficient use of 
laser power by permitting a large measurement volume, and since the 
sample volume is located some distance away from the sample plat- 
form, uncertainties associated with platform contamination or sur- 
face losses are minimized. However, the main driving force for 
the lidar configuration is that it minimizes the laser-induced OH 
interference problem (e.g., see previous section) by permitting a 
large area sample beam thereby reducing the energy density, E, at 
the sample site. The principal disadvantages of the lidar config- 
uration are its much greater susceptibility to interference from 
scattered solar flux and calibration. 
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Direct calibration of lidar systems is difficult at best 
even when dealing with stable molecules. No satisfactory approach 
is currently available for calibration when dealing with a short- 
lived species such as OH. Therefore indirect methods are general- 
ly employed to extract quantitative OH concentrations from the 
return signal. 

Three groups are currently employing or proposing the lidar 
configuration for OH measurements. The group at the Ford Motor 
Company has pioneered the use of the lidar configuration with a 
number of flight tests on the NASA CV-990 aircraft. The group at 
Goddard is currently developing the lidar configuration for both 
OH and NO detection (National Aeronautics and Space Administration, 
1983) and the group at the Jet Propulsion Laboratory is developing 
a lidar system for OH detection. 

(a) Ford Motor Co.: The group at Ford Motor Co. was the first 

to recognize the advantages offered by the lidar configuration 
and, in fact, implemented the lidar approach for ambient OH mea- 
surements. Their early research and development efforts have led 
to the only aircraft flight tested (e.g., on NASA's CV-990) lidar 
system for OH measurements. A schematic of this lidar system is 
shown in figure 5, and is described in more detail by Wang et al. 
(1981). The salient features of the Ford lidar system are: (1) 
282-nm excitation with detection of OH fluorescence at 309 nm, (2) 
broadband detection (e.g., = 5 nm) of OH fluorescence, (3) termina- 
tion of receiver field of view (FOV) to reduce solar background, 

(4) normalization of the OH return signal to the nitrogen Raman 
scattering and (5) use of return signals from an "on" and "off 
wavelength to correct for NRF . 

The excitation laser beam is generated by a Nd-YAG laser 
pumped tunable dye laser system doubled into the ultra-violet 
(UV). A portion of the UV beam is split off to pass through a 
water vapor discharge cell where the OH fluorescence from the 
cell is monitored to establish that the laser is properly tuned. 

A Lyot depolarizer is used to depolarize the laser output. The 
sample volume is illuminated by an expanding beam during a pre- 
selected gate interval. The delay of the gate interval is typi- 
cally chosen so that fluorescence induced by the laser pulse is 
monitored only for distances greater than 8 m and less than 16 m. 

A Newtonian telescope pointing along the path of the expanding 
laser beam collects the fluorescence generated. The return sig- 
nal is imaged through a sequence of filters centered at 308 nm 
and having a spectral band pass of about 3 nm. For each laser 
pulse, a portion of the return signal is split from the tele- 
scope's output and directed through interference filters centered 
at 302 nm to a second photomultiplier. This provides a reference 
signal proportional to nitrogen Raman scattering for each OH fluor- 
escence return. The OH fluorescence is referenced to the Raman 
signal thereby providing a methodology for indirect calibration 
of the lidar system. The laser system is typically tuned on and 
off the OH fluorescence line with a temporal separation of about 
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Figure 5. Wayne State Univ./Ford 
Motor Co. lidar system. 


1 second. The on and off return signals are used to correct for 
NRF . To correct for the background component of the return sig- 
nals that is unrelated to the laser wavelength (e.g., solar back- 
ground), the detection electronics are gated on again several 
microseconds after each laser pulse. The resulting signal is due 
primarily to solar background. To minimize the solar background, 
the line of sight of the telescope is terminated by a backstop. 

Use of the backstop is one of the features that distinguishes it 
from the Goddard approach and increases its flexibility in that 
the solar background noise is not a strong function of solar 
zenith angle or platform altitude. On the CV-990 the backstop 
is provided by the wing tip of the aircraft. 

The OH concentration obtained from the Ford system is sub- 
ject to the uncertainties of the measurement of the OH fluores- 
cence relative to the N 2 Raman scattering. The OH fluorescence 
efficiency is only known to within + 35%. Nonethless, referencing 
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the OH fluorescence to nitrogen Raman scattering assures consis- 
tency from measurement to measurement, and any improvement in the 
accuracy of the value for OH fluorescence efficiency can be trans- 
lated directly to the OH measurements. 

(b) Goddard Space Flight Center: Two concepts for the lidar 

approach to implement the LIF technique for measurements of am- 
bient tropospheric OH are currently under consideration at the 
Goddard Space Flight Center. The two system concepts differ in 
the bandwidth used to spectrally resolve the OH fluorescence re- 
turn signal and the excitation wavelength. A broadband width 
(i.e., = 5 nm) approach is currently in a brassboard configuration 
for aircraft test flights. A narrow-bandwidth approach with 308-nm 
excitation is currently in the conceptual evaluation phase although 
the two concepts share many common lidar components. 

The similarities between the Ford and the Goddard broadband 
systems are 282-nm excitation with 308-nm detection and indirect 
calibration via N 2 Raman scattering. An aspect of the Goddard 
lidar concept which distinguishes it from the Ford Motor Company 
approach is the method used to reduce sky background. The Goddard 
system is pointed upward with the background filtered through an 
adjustable polarizer. Looking upward reduces the solar background 
because less atmosphere is present to scatter light toward the re- 
ceiver, and the polarizer takes advantage of the fact that Rayleigh 
scattered light (the source of the background) is strongly polar- 
ized while fluorescence from the hydroxyl is not. Insertion of a 
polarizer decreases the collection efficiency for the fluorescence 
by a factor of two while achieving a much greater reduction in the 
background thereby resulting in a net improvement in the signal- 
to-noise ratio. 

The basic lidar systems used for both the narrowband and broad- 
band filtering concepts share many common components. The common 
source of both instruments is a frequency-doubled tunable dye 
laser pumped by a doubled neodymium-YAG laser. The frequency is 
checked by splitting a small amount of the 282-nm radiation off 
and passing it through a fluorescence cell in which OH molecules 
are produced by the photolysis of water with a xenon continuum 
lamp. For the broadband system the collected radiation is direct- 
ed into a 5-nm bandwidth spectral selection system which isolates 
the Rayleigh scattered return at 282 nm and the OH fluorescence 
signal centered at 308 nm. Before processing, the 308-nm radia- 
tion is passed through a pair of adjustable multilayer dielectric 
Brewster angle polarizers to minimize the background of polarized 
scattered light and a 5-nm spectral filter centered at 308 nm. 

Both the 282-nm Rayleigh return and the 308-nm OH signal are de- 
tected using photomultiplier tubes and processed using either a 
gated photon-counting system or a gated integrator depending upon 
the signal levels. 

Scattered solar flux is by far the most serious noise source 
for this system. At low altitudes much of the stray light has 
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been multiply scattered which decreases the polarization of the 
scattered light thereby increasing the solar contribution. Maxi- 
mum polarization of scattered solar radiation occurs when viewing 
90° to the solar direction and with the platform at the maximum 
available altitude. Consequently a configuration having a fixed 
viewing direction can best take advantage of the polarization of 
skylight when the sun is low in the sky and limits the usefulness 
of such a system to mid to upper tropospheric altitudes and 
early morning or late afternoon operation. 

The Goddard lidar configuration tends to minimize the fluor- 
escence interference from other sources because there are no for- 
eign objects or sources of contamination near the measurement vol- 
ume (i.e.^ no flow tube walls). However, because the current sys- 
tem has no provisions for rapid tuning or double pulsing between 
on-line and off-line operation, measurements in heavily loaded air 
(below the boundary layer or in plumes) are precluded. 

The proposed narrow -bandwidth lidar system differs from the 
broadband system primarily in the choice of the excitation and 
detection wavelengths and the use of a multiple Fabry-Perot etalon 
for spectral discrimination. The proposed narrow-bandwidth system 
would use 308 nm excitation instead of 282 nm. Detection would 
be from the same manifold at 309 nm. Since the individual fluor- 
escence lines arising from the hydroxyl molecule are approximately 
0.001 nm wide, one could in principle build a single line detector 
with a bandpass 5000 times narrower than that of the broadband sys- 
tem (i.e., 5 nm/.OOl nm = 5000). 

The principal advantage for the lidar configuration gained 
from operating with a reduced bandpass is the reduction of the 
background. The 5000x reduction in bandpass would provide a 
comparable reduction in the scattered solar flux, and therefore 
in the high background limit the noise is reduced by approximately 
70. In addition the Qi ( 1 ) line (307.84 nm) of the (0,0) manifold 
lies in an absorption feature of the solar spectrum, further reduc- 
ing the intensity of the solar scatter by an altitude-dependent 
factor which can be as large as 10. Additional benefits arise 
from the fact that the absorption cross-section in the (0,0) mani- 
fold is about 3 times larger than for the (0,1) transitions. This 
results in a direct factor of three increase in signal. Finally, 
the ozone absorption cross-section at 308.6 nm is about 30x small- 
er than that at 282 nm, and the 0(1 d) quantum, yield is about 0.55. 
These two facts combine to reduce the laser-produced OH by a fac- 
tor of at least 50. 

The overall benefits of the narrowband detection scheme are 
offset to some degree by the lower number of fluorescence photons 
arising from a single emission line compared to the number in a 
whole band. Although no studies of the yield of individual lines 
of OH have been performed, it is estimated that the Qi ( 1 ) line 
represents about 4% of the total band strength. Accounting for 
this reduction in fluorescence photons results in an estimated 
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10-fold improvement in overall performance if solar scatter con- 
tinues to dominate the noise. 

(c) Jet Propulsion Laboratory; To date all the field experi- 
ments using LIF and lidar systems to detect OH radicals in the 
atmosphere have used frequency-doubled dye lasers operating near 
282 nm. The primary reason for using this excitation wavelength 
is that the resulting fluorescence is centered near 309 nm which 
is significantly red-shifted from the laser wavelength so that 
relatively simple filtering techniques can be used to discriminate 
between elastically scattered laser light and fluorescence. In 
this section a 308-nm excitation scheme as proposed by the JPL 
group ( Laudenslager et al., 1980, 1981, 1982, 1984; McDermid et 
al., 1983a, b, c, d) will be discussed in terms of its potential 
for reducing NRF and laser-induced OH interferences and increasing 
the OH fluorescence signal. The components required for implemen- 
tation of this scheme and some of their advantages are also 
considered. 

Figures 1 and 2 illustrate schematically the 282-nm and 308- 
nm excitation schemes. In the 282-nm scheme OH is excited to the 
v'=l level; different rotational transitions correlating with low 
rotational levels of the ground state can be used for excitation. 

In order to be detected at 309 nm the v'=l level must be relaxed 
by collisions with atmospheric molecules, principally N 2 and 02 r 
to the v'=0 level. Both V-T and R-T processes occur with almost 
unity collision efficiency and with multiquantum (R) transitions 
showing appreciable probability, ultimately producing close to a 
Boltzmann population distribution in the excited state. Compet- 
ing with relaxation is electronic quenching which drastically re- 
duces the fluorescence efficiency, especially in the troposphere. 

It has been shown that the rate constants for electronic quenching 
show a dependence on the rotational and vibrational state (McDermid 
and Laudenslager, 1980, 1982; Copeland and Crosley, 1984). Thus, 
if the fluorescence intensity at 309 nm is to be used to infer 
hydroxyl radical concentrations, the rate constants for all the 
state-to-state energy transfer processes must be known and incor- 
porated into a kinetic model of the excited state. This is an 
extremely complex problem and is directly impacted by the accuracy 
and completeness of the data set of state-to-state rate constants. 
For excitation at 308 nm directly into the v'=0 level the effects 
of vibrational relaxation can be eliminated from the analysis. 
Quenching and rotational relaxation must still be considered but 
far fewer levels are involved, allowing more detailed and precise 
calculations to be made (McDermid et al. , 1983b) . 

Using excitation at 308 nm immediately realizes an increase 
in fluorescence since the absorption cross-section is some 3 to 4 
times larger than at 282 nm. However, because of the proximity 
of the laser and emission wavelengths, a high-resolution, high-con- 
trast spectral filter is required. Suitable devices have already 
been developed for remote sensing experiments, particularly wind 
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speed measurements (Rees et al., 1981a, 1982a, b; Killeen et al., 
1982; Hays et al., 1981), and could readily be adapted to an OH- 
lidar system. 

Figure 6 shows a detailed energy level diagram pertinent to 
the 308-nm excitation scheme envisioned by the JPL group. Exci- 
tation is carried out in the Q 21 I rotational transition at 307.847 
nm to the Fo manifold of the excited state. For slow rota- 
tion (low J), Hund's case (a) is an appropriate description of the 
state; for higher rotational levels in the state and 
for all levels of the state, Hund's case (b) is appropriate. 
(See Appendix A.) The consequence of this is that the so-called 
O 21 I satellite line is essentially fully allowed and has a similar 
transition probability to the main 0^1 line. It has been observed 
(McDermid et al. , 1983c; Stepowski and Cottereau, 1981) that the 
F^ and F 9 manifolds of the state do not mix appreciably, even 
after collisional relaxation. This result is not surprising from 
consideration of the coupling between the angular momentum J and 
the spin S since the process Fx <-> F 2 essentially involves a re- 
versal of the angular momentum vector. This fact has very impor- 
tant consequences regarding the fraction of the total fluorescence 
which can be detected with a narrow spectral bandwidth filter. 
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The 308-nm excitation scheme 
prepared by JPL. 



There is a fortuitous overlap, within AX < 0.001 nm, of the rota- 
tional transitions O 22 , Q 122 , 023, and Qi23 at 308.986 nm, at which 
wavelength the JPL group proposes to set their detector. It can 
readily be shown from considerations of the relative rotational 
transition probabilities that almost 20% of the total fluorescence 
from a Boltzmann distribution in the F 2 of the state will 

be detectable in a 0.001-nm spectral bandpass at 308.986 nm. How- 
ever, if excitation is to the manifold then only 3% of the 
total fluorescence would be observed. This indicates the impor- 
tance of exciting in one of the satellite branches such as 021* 

A key element in the 308-nm scheme is the use of a triple 
Fabry-Perot etalon filter with an imaging photon detector (IPD) 
(Rees et al., 1980, 1981b). The spectral resolution of the 
filter is «0.001 nm and the imaging detector allows simultaneous 
monitoring of a number of different (relatively closely spaced) 
wavelengths. Thus it is possible to measure the OH fluorescence, 
a rotational Raman normalization signal, and the background level 
simultaneously on a single laser pulse. This eliminates the sta- 
tistical correlation problems discussed elsewhere in this report 
and also avoids the need for two laser systems as in the Georgia 
Tech system. Also, the interferences caused by broadband spec- 
tral sources, such as NRF and solar scatter, are greatly reduced. 
For example, detector bandwidths of approximately 5 nm have been 
used in the 282-nm scheme. By employing a 0.001-nm filter these 
interferences would be reduced by a factor of 5000. It is inter- 
esting to note that, due to the presence of OH in the Sun, all OH 
fluorescence lines correspond to absorption features in the solar 
spectrum. The absorption in the solar spectrum at 308.986 nm is 
particularly strong because of the overlapping of a number of 
lines and this leads to an additional reduction in the level of 
solar scatter by a factor of «10 thus giving a total reduction 
in the interference level from solar radiation of 5 x 10^. 

The artificial generation of OH through photolysis of ozone 
as noted for 282-nm excitation can occur in the 308-nm scheme but 
at a significantly reduced level. The absorption cross section 
for ozone at 308 nm is 1.28 x 10“^^ cm^ and at 282 nm is 3.7 x 
10“13 cm2. The quantum yield for production of 0(^0) is also 
reduced at 308 nm, (})[0(3 d)] = 0.7, compared with 282 nm, 4)[0(1 d)] 

= 0.9, at 295 K. At higher altitudes (lower temperatures) the 
quantum yield at 308 nm falls rapidly. Thus the net interfer- 
ence level is reduced by a factor of “40 in the 308-nm scheme and 
this eliminates the need for special techniques to control this 
interference. For the same level of artificial OH interference 
some 40 times more laser energy can be used, within the limits 
posed by saturation effects (see Appendix B) , with a proportional 
increase in laser-induced fluorescence. 

In the 308-nm scheme the energy gap between the excitation 
and detection wavelengths is «100 cm“^ and thus vibrational 
Raman scattering by any atmospheric component is not a potential 
interference. However, rotational Raman scattering by nitrogen 
and oxygen has been considered even though the scattering cross- 
sections are very small ( «10“29 ^,^2 sr“l). (There is some 
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uncertainty regarding the absolute magnitude of the Raman cross- 
section for N 2 at wavelengths shorter than 350 nm.) For excita- 
tion in low , Q 2 I' ^1 rotational lines there are no overlaps 
in the proposed detection region near the Q 2 bandhead. The use 
of nearby Raman lines for normalization is attractive, as noted 
above . 

A promising coherent source for excitation at 308 nm is the 
XeCl excimer laser. Excimer lasers represent a new class of laser 
which operates directly in the UV and has good efficiency (>2%) . 
Excimer lasers are high-pressure (typically 1 to 5 atmospheres) 
transversely excited electric discharge lasers somewhat similar 
to the TEA CO 2 laser. Transversely excited gas lasers are inher- 
ently very rugged and reliable devices if designed properly. To 
date, several excimer lasers have been operated for many pulses 
with little loss of output energy on a sealed gas system using 
fast gas recirculation. Some chemical and particulate cleanup 
may be necessary, and for the case of the rare-gas/f luoride lasers 
some gas makeup may also be required. Lasing occurs from bound 
excited molecules to the ground states of the rare-gas/halide 
molecules which are unbound or, as is the case for XeCl, very 
weakly bound. The bound-free transition gives rise to a broad 
spectral gain bandwidth, typically 1 to 2 nm wide for the excimer 
lasers. These lasers are scalable to high pulse energies in the 
ultraviolet (0.1 to 40 J) and with gas recirculation can be 
operated at high repetition rates (1 to 1000 Hz). The temporal 
pulse widths of excimer lasers can also be varied from picosec- 
onds to several hundred nanoseconds. 

Excimer lasers can be tuned directly over their gas bandwidth 
using intracavity dispersive optics to obtain narrow spectral out- 
put (Pacala et al., 1982, 1984). Table 5 lists the tuning ranges 
and operational characteristics of the most efficient excimer 
lasers. Although the various developed excimer lasers do not 
offer total wavelength coverage in the UV, the high energy output 
from the excimer laser allows for efficient stimulated Raman 
shifting of the tuned output to other UV, visible and VUV wave- 
lengths. At JPL the XeCl excimer laser has been operated with 
narrow spectral bandwidths (<0.001 nm) (Pacala et al., 1984) and 
has demonstrated direct tuning to several absorption lines of the 
OH, A-X, 0-0 transition including the Q 21 I line discussed above. 

A continuous tuning range from 307.5 to >308.5 nm has been demon- 
strated at JPL for XeCl. It is interesting to note that the XeBr 
laser operates in the spectral region suitable for 282 nm excita- 
tion of OH. The XeBr laser has a gain region from 281.3 to 282.3 
nm and pulse energies of 100 mJ pulse”^ have been obtained from 
such a laser (Sze and Scott, 1978). 

The excimer laser source should be customized to meet the 
needs of atmospheric measurements. In order to achieve long- 
lived operation on a self-contained gas mixture the laser body 
must be constructed from the proper materials and some means of 
gas cleanup should be incorporated into the flow system. To ob- 
tain high-repetition-rate operation the design of the electronic 


2 8 



TABLE 5. WAVELENGTH RANGE AND OPERATIONAL CHARACTERISTICS OF SELECTED EXCIMER LASERS 


TYPE OF 
LASER 

WAVELENGTH EFFICIENCY 
(nm) (%) 

ENERGY 

(J) 

PULSE 

WIDTH 

(ns) 

REPETITION 

RATE 

(Hz) 

COMMENTS 

EXCIMER LASERS 








DISCHARGE 








ArCl 


175 





MANY X's EACH TUNABLE OVER 

ArF 


193-194 





1.0 TO 2.0 nm. NARROW SPEC- 

KrCl 


221-223 





TRAL WIDTH. SINGLE MODE HAS 

KrF 


248-249 

0.1-4 

0.1-40 

1-250 

1-1000 

BEEN DEMONSTRATED, FROM OSCIL- 

XeBr 


282-282 





LATOR-AMPLIFIER, 1 kHz REPE- 

XeCl 


307-309 





TITION RATE REQUIRES GAS RE- 

XeF 


351-353 





CIRCULATION. 

EXCIMER PUMPED DYE 








(XeCl-308 nm) 








FUNDAMENTAL 


320-800 

0 . 1 - 0.2 

< 0.12 

1-100 

1-1000 

A LONGER VISIBLE WAVELENGTH 

FREQUENCY DOUBLED 


220-320 

0.01 

<0.08 

1-50 

1-1000 

EXCIMER WOULD BE MORE EFFI- 








Cl ENT AS A DYE PUMP. 

RAMAN SHIFTING OF 

STOKES SHIFTS 






(XeCl LASER 308 rm) 

1ST 

2ND 3RD 







353 

414 499 





TUNING OF EXCIMER LASER PRO- 

D 2 

339 

377 426 





DUCES TUNABLE STOKES OUTPUT. 

CH 4 

338 

375 422 






N 2 

337 

359 392 

<50 




ANTI-STOKES SHIFTS ALSO POS- 

T1 


405 





SIBLE AT LOWER EFFICIENCY. 

Pb 


459 






Ba 


475 





RAMAN OSCILLATOR-AMPLIFIER 

B1 


475 





CONFIGURATION IMPROVES EFFI- 








CIENCY. 


circuitry is critical. The concept of magnetic saturable reactor 
switches has been developed at JPL to solve this problem. In 
order to obtain narrow-bandwidth tunable operation, an excimer 
laser oscillator- amplifier arrangement is necessary. Several ex- 
cimer lasers with Raman shifting cells will probably be suitable 
for field measurements of multiple species. The design of these 
components is very much constrained by the parameters required 
for the particular atmospheric measurement. 

In summary, compared to the 282-nm scheme, the 308-nm scheme 
as proposed by the JPL group would yield an increased fluorescence 
signal and a reduction in the level of interferences by >10^ thus 
providing a realistic sensitivity at the [OH] = 1 x 10^ molecule 
cm"3 level. Also, the availability of simple, compact and effi- 
cient tunable XeCl excimer lasers operating near 308 nm offers ad- 
vantages for the field instrument. Detection of OH in laboratory 
experiments has been accomplished using these lasers. The ability 
to tailor the energy, repetition rate and pulse length of these 
lasers made them especially suitable for application to the OH 
measurement problem. 
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In Situ Systems at Atmospheric Pressure . Three systems or ap- 
proaches for implementing the in situ LIF technology at ambient 
pressures for measurements of OH were discussed as part of the 
workshop. These systems are now described. 

(a) Georgia Institute of Technology: 

In situ sampling at atmospheric pressures represents one of 
the most direct and earliest approaches for implementing the LIF 
technique for OH measurements. One of the first groups to initi- 
ate this approach from both ground based and airborne platforms 
was from the University of Maryland and is now at Georgia Tech. 
This section will describe the current state of technology of the 
Georgia Tech LIF in situ system. 

The basic excitation scheme of the Georgia Tech approach 
consists of 282-nm excitation with detection of OH fluorescence 
at 309 nm using an = 5-nm narrow-bandwidth filter. The LIF system 
at Georgia Tech has been designed to minimize two of the dominant 
noise sources plaguing in situ ambient pressure systems; i.e., 
laser-generated OH and broadband nonresonant fluorescence. The 
method used for reduction of OH interference has been that of re- 
ducing the laser pulse duration as was first suggested by Wang et 
si* (1976) and Hanabusa et al. (1977) and then again by Davis et 
al. (1981b) . The rationale for short pulses was discussed in the 
previous section. In particular, if the excitation pulse can be 
made sufficiently short, then the lifetime of the 0 (^ 0 ) (and con- 
sequently the OH) produced by photolysis of O 3 lags the laser 
pulse. The approach to minimize broadband nonresonant fluores- 
cence is, in part, based upon results obtained by Davis et al. 
(1981a, b) during the GAMETAG missions which clearly indicated 
that one of the dominating noise sources was temporally uncorre- 
lated background fluorescence. One obvious approach to correct 
for a broadband background noise source is to measure the signal 
+ background at Aon ~ 282 nm and the background only for Agff, 
slightly different from Aqj^. However, because of the uncorre— 
lated nature associated with this background noise, the Xoff 
signal must be obtained close enough in time to effectively freeze 
the atmospheric component contributing to the fluorescence. 

The optical hardware for the Georgia Tech system is illus- 
trated in top and side views in figure 7. The sample probe is 
illustrated in figure 8 . The salient features of this configura- 
tion are (1) the use of two Nd:YAG laser pumped, frequency-doubled 
dye laser systems for generation of Aon ^off wavelengths; 

( 2 ) pockel cells to achieve short pulses; (23) 282-nm excitation 
with 308-nm detection and (4) an OH calibration system. The two 
complete laser systems provide the mechanism for obtaining the Aon 
signal (e.g., OH signal plus background) and the Aoff signal 
(e.g., background only) with control of the temporal separation 
of the two signals down to less than 1 msec. With the use of the 
pockel cells, UV excitation pulses in the range of 0.15 ns to 4.5 
ns can be achieved starting from a basic laser pulse width of 9 
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DRP - Dye Reservoir and Pump 

DY - Frequency Doubled Nd:YAG Laser 

END - Energy Monitor Diode 

FCS - Frequency Calibration System 

FD - Fast Diode 

FSO - Fast Sampling Oscilloscope 

OCA - Optical Coupling Arm 

P - 90° Prism 

PDL - Pulsed Dye Laser 

SR - Seat Rail 

BC - Beam Combining Mirror 

COS - Complex Optical System 

(a) Top view. 
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AMP - Dye Laser Amplifier OSC - Dye Laser Oscillator 


BS 

- Beam Splitter 

OT 

- 2' X 6' Optical Table 

DY 

- Frequency Doubled NdrYAG Laser 

P 

- 90° Prism 

DRP 

- Dye Reservoir and Pump 

PC 

- Pockel Cell 

FCS 

- Frequency Calibration System 

PDL 

- Pulsed Dye Laser 


FSO - Fast Sampling Oscilloscope SHG - Second Harmonic Generation Crystal 

GLP - Gian Laser Polarizer SR - Seat Rail 

LTSG - Laser Triggered Spark Gap VI - Vibration Isolator 

H - Mirror 

(b) Side view. 

Figure 7. Georgia Tech OH system. 
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Figure 8 . Georgia Tech OH sample chamber. 


ns. The currently demonstrated reduction in the laser-generated 
OH O 3 /H 2 O interference utilizing a pulsewidth of 0.5 ns is a 
factor of 10 relative to that observed from a 9-ns pulse. 

An important component of any in situ sampling system is the 
sample chamber and associated inlet plumbing. This is particular- 
ly true for a reactive trace species such as OH. The Georgia 
Tech OH sampling manifold shown in figure 8 has been designed for 
aircraft operation and can be completely sealed off from the air- 
craft cabin. Full aperture gate valves on both ends of the OH 
fluorescence chamber permit the isolation of the chamber while 
the aircraft is on the ground, thus preventing excessive contami- 
nation of the chamber liner. These gate valves also serve to 
isolate the fluorescence chamber from rain and snow while operat- 
ing under severe flying conditions. In the latter regard when 
water droplets do appear on the laser injection window, a third 
gate valve can be closed and a pressure equilibration valve 
opened. At this point, the laser injection window can be taken 
off and cleaned without disrupting the aircraft cabin pressure. 

The intake manifold consists of a 3.5-in. i.d. black anodized 
pipe (Teflon coated) which is bent through an arc of ~ 50“ . For 
typical tropospheric conditions the lifetime of OH is on the 
order of 1 to 2 sec and the radial diffusion time at tropospheric 
pressures is on the order of 90 to 200 msec. Since the time of 
transit of an OH molecule to the point of sampling for an aircraft 
air speed of 130 knots is 30 ms, OH losses due to homogeneous 
chemistry are negligible. Note also that the fluorescence chamber 
shown in figure 8 has been designed to handle as many as 12 PM 
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tubes and associated optics thereby permitting efficient use of 
the available excitation energy. The use of multiple PMT's is, 
in effect, similar to the use of a multipass excitation chamber, 
but, based on the Georgia Tech group's experience, results in 
lower total noise levels than a multipass system. 

One of the more attractive features of the in situ sampling 
systems is the potential for direct absolute calibration in terms 
of known OH concentrations. Figure 8 also shows the position of 
a calibration chamber in which known amounts of OH can be gener- 
ated, although at levels well above ambient for calibration. A 
gas injection port (GIP) for injection of a gas to scavenge am- 
bient OH from the air stream is also available. 

Although several noise sources contribute to defining the 
operational signal-to-noise ratio for the aircraft system, fluor- 
escence from the walls of the sample chamber is an important one 
directly associated with the design of the chamber. This noise 
source is derived from scattered 282-nm radiation being absorbed 
by small amounts of contaminants on the walls of the chamber. The 
photon noise from N 2 and CH 4 Raman scattering and aerosol fluores- 
cence is normally much lower than that of wall fluorescence under 
clean tropospheric air conditions. However, when sampling in con- 
tinental boundary layer air over geographical regions where sig- 
nificant anthropogenic sources are present, the atmospheric aero- 
sol loading is at least 10-50 times higher than in remote tropo- 
spheric boundary layer air, and the corresponding noise signal 
approaches, and in some cases exceeds, that from wall fluorescence 
The solar flux noise source is typically three orders of magnitude 
lower than that resulting from wall fluorescence. This is a di- 
rect result of both the design characteristics of the collection 
optics and the use of time gating on the signal processing elec- 
tronics . 

(b) Harvard University; 

As with the other LIF systems discussed in this report, a 
central concern of the Harvard group is the need to minimize 
photochemical OH interference while providing a system that can 
be directly calibrated for response to known amounts of OH and 
for evaluating system and background noise in the absence of OH. 
The basic scheme utilizes 282-nm excitation with detection of the 
309-nm fluorescence in an in situ sampling configuration. The 
distinguishing features of the system proposed by the Harvard 
group are ( 1 ) the use of a high pulse rate, low pulse energy cop- 
per (Cu) vapor laser source; and (2) a system configuration for 
calibration and noise measurements. 

Figure 9 illustrates the system proposed by the Harvard group 
As noted earlier, laser-induced OH is, among other parameters, de- 
pendent upon the laser pulse width and the peak flux. In particu- 
lar, it is generally recognized that to minimize this OH interfer- 
ence, a short laser pulse and/or low peak flux is desirable. Con- 
versely, to maximize system sensitivity, high average powers and/ 



Cu vapor laser 



Figure 9. Schematic drawing of proposed Harvard Univ. optical system. 


or high peak flux are desirable. To accommodate these conflict- 
ing goals, the system approach suggested by the Harvard group is 
the use of a high-repet it ion -rate pulsed laser to provide high 
average powers with low peak flux. In particular an atomic cop- 
per vapor laser is proposed as the source to pump a dye laser 
which is in turn doubled into the ultraviolet spectral region. 

This is in fact the approach that has been developed by the 
Harvard group for balloon-borne stratospheric OH measurements. 

The proposed Cu-laser system would operate at rates in the 15 to 
20 kHz range at 510.6 nm. This system would provide an average 
output power of 20 W, corresponding to 1 mJ /pulse and a peak power 
of 50 kW with a 20-nsec pulse width. This would be used to pump 
a dye laser for a projected output at 282 nm with a spectral band- 
width of <5 GHz and an average power of 50 mW. 

The advantages expected from such a system are illustrated in 
figure 10 (Stimpfle and Anderson, 1982), where R, the ratio of the 
fluorescence from OH formed photochem.ically in the laser pulse to 
that from ambient OH, is plotted against the pulse repetition 
frequency at a fixed average power. The values of R shown here 
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Figure 10 . Influence of varying the pulse repetition at a fixed 
average energy. 


are calculated assuming that each pulse interrogates a different 
sample volume element. Note that for a fixed average power, the 
peak laser flux must decrease with increasing repetition rate. 

It is clear that under the assumed conditions, the fluorescence 
from laser-induced OH is drastically reduced at the higher repeti- 
tion rates. The dependence of the OH interference for the Harvard 
approach was discussed in the previous section and given by equa- 
tion ( 10 ). It was noted in the previous section that the OH in- 
terference could be reduced as much as desired (e.g., E approach- 
ing 0 ) but at the expense of the signal-to-noise ratio. Accord- 
ingly it is important to realize that the results illustrated in 
figure 10 are constrained by the requirement that each laser 
pulse interrogate a new sample volume element. To realize this 
the characteristic sample element length and the laser repetition 
rate must be consistent with the time required for the sample air 
flow through the sample chamber. 

The sample chamber proposed in the Harvard system is an aero- 
dynamic pod mounted exterior to the aircraft fuselage. The pro- 
posed design of the pod would include a multipass cell for increas- 
ing traversals of the UV beam through the scattering volume and 
PMT's with interference filters for wavelength discrimination. 

The proposed pod would permit injection of C3H5, O3 and H2O at 
the entrance of the pod during flight. Injection of the C3H6 
would be designed to scavenge ambient OH from the airstream there- 
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by permitting inflight background noise measurement in the absence 
of ambient OH. Injection of O 3 and H 2 O would provide the possi- 
bility of evaluating the susceptibility of the system to laser-in- 
duced OH. And finally, a Lyman a lamp coupled with a light trap- 
ping cell would be located upstream of the detection axis as a 
means to photodissociate H 2 O to produce OH. This would provide a 
mechanism to demonstrate that the OH detection system was operat- 
ing and give a rough check on the OH sensitivity. 

In Situ Systems at Low Pressure . One of the approaches that is 
currently being pursued at Portland State University and Institut 
for Chemie, KFA, Julich, for minimizing background noise is re- 
ferred to as baric filtering. The group at Portland State Univer- 
sity refers to this technique as fluorescence assay with gas ex- 
pansion (FACE). (See, e.g.. Hard et al., 1979, 1980.) In this 
mode of operation the basic OH excitation scheme can utilize 
either the 282-nm or the 308-nm excitation scheme described above. 
The salient feature of the baric filtering is that the excitation 
of OH occurs in a partial vacuum. Because of this, one of the 
more interesting aspects of this approach is that there is typi- 
cally a time delay between the excitation and the fluorescence 
event. This offers the potential for temporal discrimination be- 
tween the OH fluorescence and background noise. Much of the back- 
ground interference occurs simultaneously with the laser pulse 
while the temporal behavior of the OH fluorescence is controlled 
by radiative lifetime and collisional quenching. At low pressure, 
the OH fluorescence is delayed with respect to the excitation 
pulse. Moreover for most molecules nonradiative deactivation pro- 
cesses dominate at atmospheric pressures, and under these condi- 
tions the fluorescence efficiency is low. If the gas sample is 
excited at low pressures, the increase in fluorescence efficiency 
compensates for the decrease in density of the target species 
thereby leaving the fluorescence per unit excited volume nearly 
constant. For OH the radiative lifetime of the A^z state is 
approximately 700 nsec. In the pressure range of a few torr, net 
lifetimes of 100 to 300 nsec are achieved, permitting the use of 
time-delayed detection to suppress the early background due to 
Rayleigh, Mie and Raman scattering and rapidly decaying NRF. This 
time delay reduces the need for spectral filtering, especially 
with 308-nm excitation, where the required spectral resolution 
would otherwise be high but where the excitation cross section is 
largest. 

Figure 11 illustrates a typical timing sequence for the baric 
filtering approach. The PMT is gated off during the laser excita- 
tion pulse to avoid saturating effects. After the laser pulse de- 
parts from the excitation cell, the PMT is switched on. Still 
later, after a total delay that is adjustable, the gated charge 
integrator or gated photon counter is switched on for an adjust- 
able duration. The optimal gate time is dependent upon the sig- 
nal strength for the background and the OH fluorescence. The 
best SNR that can be achieved by gating out a single fluorescence 
background is 
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/2kBAs - 4 \ 1/2 

SNRmax = s(— — 

\2kB/ks - 3 / 

where the k's are the observed fluorescence decay rates for the 
background (B) and signal (S), and S is the OH signal strength. 
This reduces to the signal-limited case (SNR = for kB>>kg. 

It is significant that this limit can always be theoretically 
reached, whenever kB>>2ks. Furthermore, the maximum SNR rises 
very rapidly with kB/ks, as shown in figure 12, reaching 80% of 
the limiting value when the rate constant ratio R has increased 
from 2 to 3. 

Another advantage of baric filtering is the reduction of pho- 
tolytic interference, caused by laser-induced photolysis of other 
molecules to produce OH. While laser-induced OH interference is 
still present, transient kinetic effects (e.g., excitation photons 
leaving the detection zone before an appreciable quantity of OH is 
produced) reduce its dominating characteristics. The Portland 
State University group has experimentally verified the absence of 
photolytic interference in their system by sampling room air at 
45% RH, enriched with ozone to a concentration of 3.9 ppm. After 
expansion to 4 torr, no OH signal was observed at the 10® cm“3 
level, with or without delayed gating. When the delay was removed 
and the final pressure was raised to 150 torr, photolytic inter- 
ference equivalent to an ambient OH concentration of 8 x 10^ cm“3 
was observed. From the latter observation, a worst-case interfer- 
ence at 4 torr is calculated to be equivalent to 1 x lo5 OH cm“3 
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Figure 12. Effectiveness of optimally delayed gating 

vs. ratio of first-order decay rate constants 
of a fluorescent background B and signal S. 

in ambient air containing 5 x lol2 03 cm“3 and 2.5 x 10^^ H 2 O cm“3. 
Tests have also been conducted using H 2 O 2 and HOMO; these gases 
are weaker absorbers at 282 nm, but they offer less temporal sup- 
pression of this interference (Hard et al., 1980). H 2 O 2 at ~ 100 
ppm and HONO at ~ 0.1 atm, with subsequent dilution by a factor 
not greater than 1000 , give detectable interference equivalent to 
10' OH cm"3. This translates to negligible interference for the 
baric filtering technique at ambient concentrations of these gases. 
Finally, since there are other possible photolytic parents, air 
from a propene/NOx smog chamber reaction was sampled in the NO, 

NO 2 and O 3 regimes, and in each case the OH signal disappeared 
after the lights were turned off. 

A particularly attractive feature of the baric filtering 
technique is the possibility of calibration of the overall system. 
Because the inlet flow rate in the baric filtering approach is 
relatively low, absolute calibration can be obtained by sampling 
from a chamber containing a known amount of OH. In a closed cham- 
ber irradiated by UV fluorescent lamps or sunlight, the OH concen- 
tration can be measured by the rate of disappearance of a more 
abundant species with which only OH reacts, provided the reaction 
rate constant is available from other work. One approach that 
has been used by the Portland State University group involves the 
use of an irradiated 150-liter Teflon bag containing a propene/NO/ 
H 2 O starting mixture. The propene concentration is measured re- 
petitively by gas chromatography with flame ionization detection. 
Because propene is also removed by reaction with ozone, the cali- 
bration is restricted to the NO-rich regime. The calibration 
procedure suggested by the group from Institut fUr Chemie consists 


38 



of generating known amounts of OH in the vicinity of the inlet 
orifice through photodissociation of H 2 O by Lyman a radiation. 
Careful design of a calibration flow system would permit a direct 
calibration of a low-pressure sampling system. 

The main disadvantage of the expansion method is the loss of 
OH by wall reactions. About 50% of the measured species is cur- 
rently estimated to be lost in the sampling process. Uncertain- 
ties in the many factors entering this calculation illustrate the 
importance of an absolute calibration as well as proper aerodynamic 
design. 


Comparison of LIF Measurement Techniques 

A significant portion of the time at the workshop for those 
participants cognizant of the LIF techniques was devoted to two 
activities; (1) the development of a central equation describing 
the LIF process and the various systems for implementation of the 
LIF technique and (2) the estimation of the performance of each 
of the systems described above. The results of these activities 
are presented in this section. 

Central LIF Equation . - The variety of LIF techniques (as well as 
others) that are currently under development or have been proposed 
is a testimony to the importance as well as the difficulty of ob- 
taining valid measurements of tropospheric OH. While the math- 
ematical basis for LIF detection has been reasonably well under- 
stood for many years, the participants of the workshop felt that 
the development of a "central equation" which encompassed the dif- 
ferent nomenclatures/forms of equations used to describe different 
OH LIF sensors would prove valuable both to those groups directly 
involved in such measurements as well as to the broader atmospheric 
chemistry scientific community. Moreover, the development of a 
"central equation" from which the working expression for each LIF 
instrument could be extracted as a specialized case would provide 
a common basis for comparisons of current and future LIF systems. 

The approach used was first to provide a spectroscopic des- 
cription of the energy levels and transitions employed in the LIF 
methods for the structural details of the OH radical. Second, 
the so-called central equation was derived by relating the expected 
count rate in a LIF system to the absolute OH concentration. Fi- 
nally, the assumptions and/or approximations necessary to extract 
a working expression for each LIF approach from the central equa- 
tion were explicitly stated. This then related the approach and 
observed quantities of an LIF technique described in the previous 
section to the central LIF equation. This effort resulted in 
agreement among the participants on the values of all parameters 
to be used relative to the spectroscopic details of- !the OH radical 
and provided a common language for the discussion of system de- 
tails. A detailed discussion of the development of the fundamen- 
tal equation resulting from this effort is presented in Appendix 
A. The central equation and the working expression for each LIF 
technique are given in figure A7 and will not be repeated here. 
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It was also recognized that these fundamental equations did 
not address the issue of saturation phenomena in LIF measurements. 
Since saturation effects are very complex and mathematical analy- 
sis quite difficult, it was decided that a brief discussion should 
be included with a particular emphasis on the status of work in 
this area. This discussion is presented in Appendix B. 

Performance Estimates . As a result of the workshop activities 
related to agreement upon the central LIF equation (discussed in 
the previous section) and the development and comparison of the 
working expressions for OH interference (discussed in an earlier 
section), the anticipated and in some cases the achieved perfor- 
mances of several LIF techniques were compared for common sample 
atmospheric conditions. Since the implementation approaches were 
in different stages of development, one estimate was made for 
current capability and another was based upon anticipated tech- 
nology improvements or proposed systems. Tables 6 and 7 present 
these performance estimates in the form of signal-to-noise ratios. 
Two test atmospheres were postulated, each of which contained a 
concentration of 5 x 10^ OH cm“3 . a tropical marine boundary 
layer was assumed with 1 atm total pressure, 25 ppbv ozone, and 
25 torr water vapor. A middle tropospheric condition was assumed 
with 0.5 atm total pressure, 60 ppbv ozone, and 0.3 torr water 
vapor. 

It has been noted earlier that the OH interference and OH 
fluorescence signals are both a function of the laser energy. 

Since the OH interference is a corrective term, the SNR can in 
principle be made as large as desired by increasing the laser 
energy. In practice however the laser-generated OH becomes too 
large to be realistically correct. Accordingly, a ground rule for 
this comparative evaluation was established such that the level of 
laser-induced OH must be less than 20 percent of the true OH con- 
centration. The level of interference calculated for each LIF ap- 
proach is given at the head of each column along the row denoted 
"Laser-Generated OH". The basis of comparison is the expected sig- 
nal-to-noise ratio for the assumed conditions. Integration times 
of 1 and 20 minutes (e.g., column 1) were assumed for this evalua- 
tion. Two levels of background NRF noise (e.g., column 2) were 
assumed. While the nature of its origin and its meteorology have 
not been established, its magnitude was assumed to be proportional 
to detector bandwidth, and values of 10^ to 10^ equivalent hydroxyl 
radicals cm“3 for a 2.7-nm bandwidth are typical for both boundary 
layer and midtropospheric atmospheres, with the smaller figure 
applying to over-ocean midtropospheric observations. The variance 
of NRF background exceeds Poissonian by a "Degradation Factor" 
(column 3) which depends upon the design of the instrument, speci- 
fically the time interval between "on-line" and "background only" 
measurements. This degradation factor is expected to be unity 
when this interval is short, and at 1 sec the Ford group has 
observed a value of 5 (unpublished results). 
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TABLE 6. SIGNAL-TO-NOISE RATIOS ([OH] = 5 x 10^ cm“3 ) 

BASED UPON CURRENT TECHNOLOGY 



However, based upon information from GAMETAG flights as to the back- 
ground levels, the Georgia Tech group feels a more conservative MDC 
for a S/N = 2 ranges from 8 x 10^ cm“3 to 1.5 x 10® cm“3 for 20- 
minute integrations. 

Poisson statistics are applicable due to short on/off resonance time; 
therefore no degradation due to NRF is expected. 

^ Actual performance level for ground-based operation with integration 
time of 30 minutes. 

For the background levels given, the minimum detectable concentration 
(MDC) for a S/N = 2 ranges from 1.8 x 10= cm“3 to 5.8 x 10= cm =. 
However, based again on information from the GAMETAG flights as to the 
background levels, the Georgia Tech group feels a more conservative 
MDC for a S/N = 2 ranges from 3 x io= cm“3 to 6 x 10= cm = for 20- 
minute integrations. 

6 Projected performance level based on actual ground-based performance 
with integration time of 30 minutes. 



TABLE 7. SIGNAL-TO-NOISE RATIOS ([OH] = 5 x 10^ cm-3 ) 

BASED UPON ANTICIPATED TECHNOLOGY IMPROVEMENTS 
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MARINE BOUNDARY LAYER (1 atm, 25 ppb O^ , 25 torr H^O) 


1 25. Oa b b 12.0 

10/ 5 ebb 6.1 

20 c b b 3.1 

1 8.5 0.1 4.8 4.4 

10° 5 c d d 2.0 

20 c d d 1.0 

_ 1 110. Oa b b 5.6 

10' 5 16 c b b 2.7 

20 c b b 1.4 

1 38. oa 0.5 22.0 1.9 

10° 5 10 c d d 8.8 

2^ c d d 4 . 4 

MIDDLE TROPOSPHERE (0.5 atm, 60 ppb O^ , 0.3 torr H 2 O) 


0.6 

0.1 

15.5 

3.1 

0.8 

13.4 

2.7 

0.7 


Laser-Generated 
OH (m.olec cm 


0 
1 — 1 

X 

I — 1 

0 
1— 1 

X 

1 — 1 

0 

i-H 

X 

in 

rH 

2.6 X 10'^ 

7.5 

1 

— 

36.0® 2.7 8.1 

35.0 

3.5 

107 5 

— 

c d d 

7.0 

0.7 

20 

— 

c d d 

1.8 

0.2 

1 

— 

12.0® b b 

12.1 

3.0 

108 5 

— 

c b b 

2.4 

0.6 

20 

— 

c b b 

0.6 

0.1 

1 

— 

154.0® 12.2 36.0 

157.0 

15.5 

107 5 

32 

c d d 

31.4 

3.1 

20 

— 

c d d 

7.9 

0.8 

1 

— 

54.0® b b 

55.0 

13.4 

108 5 

20 

c b b 

11.0 

2.7 

20 

— 

c b b 

2.8 

0.7 


a Based upon GAMETAG flight information relative to the background 
levels, the Georgia Tech group believes a more realistic estimate of 
the minimum detectable concentration (SNR =2) to be in the range of 
3.5 X 105 cm“3 to 7 x 10 cm~3 for a 20-minute integration time. 

° Background dominated by solar noise, results for 10' and lOO back- 
ground are the same. 

^ Poisson statistics are applicable due to short on/off resonance time; 
therefore no degradation due to NRF is expected. 

Background dominated by solar noise, Poisson statistics are applicable, 
® For the same reason noted in footnote (a) the minimum detectable con- 
centration (SRN = 2) is estimated to be in the range of 1.4 x 10^ cm~3 
to 3 X 10® cm“3 for a 20-minute integration time. 
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with the exception of the entries in table 6 by the Ford 
group, the SNR's are computed from the working signal expressions 
given for each group in Appendix A. Appendices C - G were provided 
by the respective groups to describe the numerical calculations 
leading to the entries in tables 6 and 7, 

Among the groups providing entries for table 6, only Ford 
and Georgia Tech have systems that are aircraft ready and have in 
fact been flight tested. The LIF techniques not specifically 
listed are not sufficiently developed to warrant comparison in 
the "Current Technology" table. The performance levels given for 
the Ford system are based upon measurements obtained from a ground- 
based platform and for a 30 minute integration time. Also, the 
Ford entries are given only for a degradation factor of 5, corres- 
ponding to the actual degradation noted by the Ford group with a 
1-second delay between their Aon ^off wavelengths. It is 

important to note that the SNR's given for the Ford system have 
not been achieved under flight conditions. Entries are given for 
the Georgia Tech system only for a degradation factor equal to 1. 
Recall that the Georgia Tech system utilizes two laser systems, 
which permits and ^Qff signals with temporal separations of 

less than 1 msec. The NRF between Xqh ^off expected to 

be correlated, resulting in no degradation in SNR with integration 
time (e.g., see footnote b to table 6). Note also that the SNR 
calculated for the test conditions established by the workshop is 
considerably larger than that actually expected by members of the 
Georgia Tech Group (e.g., see footnote a to table 6). The per- 
formance given in footnote a for the Georgia Tech system is based 
upon flight experience and in particular reflects the fact that 
the Georgia Tech Group considers the background noise adopted for 
this comparison as too optimistic. It is of interest to note 
that the predicted laser-generated OH interference is comparable 
for the two atmospheric LIF systems but is about an order of mag- 
nitude lower for the Portland State system which employs low- 
pressure excitation. 

The anticipated improvements in the performance of those 
systems listed in table 6 as well as the anticipated performance 
of proposed systems are given in table 7. The improvements in 
performance for the three systems also listed in table 6 are 
typically less than a factor of 2. In general the new technology 
improvements represent higher laser pulse frequency , providing more 
signal per unit integration interval, shorter pulse width, and im- 
proved spectral resolution for better background discrimination. 

The Ford group also included improvements in optical configuration 
which wili provide higher signal throughput. The additional three 
systems included in table 7 represent approaches that are antici- 
pated with in the next 2 to 5 years. Note, that the Goddard ap- 
proaches included in table 7 are uniquely limited by solar scat- 
tering rather than NRF background. 

Summary Comments. Although the literature is sizable and rapidly 
growing, the LIF technique is still young as a means of making 


43 



hydroxyl measurements under ambient conditions. This fact is re- 
flected most strongly in the uncertainty which exists on the ques- 
tion of whether or not the best variant of the technique has been 
identified. Remote (lidar) and in situ (atmospheric and low pres- 
sure) variants which are being investigated have been described. 

The agreement on the mathematical formalism developed at this work- 
shop should, through its usage by workshop participants and other 
investigators in their future work, aid in the presentation and 
understanding of current and new ideas for the advancement of the 
LIF techniques. From the performance intercomparisons some maturi- 
ty seems evident in the nearness with which several variants are 
approaching a common standard of performance. The temptation is 
to conclude that perhaps each of them has been nearly optimized 
by adjustment of instrument parameters to offer results as good 
as that variant of technique permits. The temptation to reach 
this conclusion should probably be resisted a little longer. 

Indirect Methods 

Indirect methods for determination of ambient OH concentra- 
tions are based upon the reaction of OH with a given tracer gas. 

The abundance of the tracer gas must be known before and after 
reaction with OH and the reaction rates must be known. Two indi- 
rect approaches were discussed during the workshop. The first, 
referred to here as the global average approach, relies princi- 
pally on the global mass balance of certain trace constituents 
and accordingly provides a globally averaged estimate of tropo- 
spheric OH. The second approach referred to here as the radio- 
chemical approach is an in situ sampling technique that relies 
on measuring the production of ^“^002 via reaction of ^^CO and OH. 
Details of both of these approaches are discussed below. 

Global Average . Effective utilization of this approach for OH 
determination requires that the tracer species satisfy a number 
of conditions: 

« Its atmospheric abundance in time and space must be ac- 
curately defined. 

• It must be known that removal occurs predominantly by re- 
action with OH and the rate constants must be accurately 
determined. 

• The sources and source distributions in time and space 
must be accurately known. 

Within this framework the uncertainties are further reduced if 
reaction with OH has modest or low activation energy and the 
molecule is reasonably long lived. 

_ Over the last decade attempts have been made with tracer 
chemicals that are either predominantly natural in origin or pre- 
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dominantly anthropogenic in origin. and are two of the 

naturally occurring tracers that have been utilized to estimate 
OH levels. The sources however are so complex and varied 

that such an approach is at best unreliable. The problem of com- 
plex sources is minimized, although not eliminated, when is 

utilized. The dominant source of is cosmic rays and this 

source is independent of season and latitude. Compared to the 
cosmic ray production rate ( «20 x 10^5 molec yr“^) anthropo- 
genic sources are negligible. is also produced in signifi- 

cant but uncertain quantities (2 to 10 x lO^S molec yr“l) by the 
oxidation of organic matter. Current atmospheric ^^CO data are 
all taken at ground level, but seasonal (figure 13) and latitude 
gradients have been noted. These limited data have been fitted 
with the help of a 2-D model (Volz et al., 1981) to suggest an 
average tropospheric OH abundance of 7 ± 3 x 10^ molec cm“3. 

These values are a factor of 3 lower than those estimated by 
Weinstock and Niki (1972) almost a decade ago, largely because of 
their underestimation of ^'^CO abundance. It must be realized, 
however, that the lifetime of l^CO is short (5 ± 2 months), and 
a finer resolution of the organic source term and vertical gra- 
dients is needed to further refine these estimates. 



14 

Figure 13. Seasonal variation of CO at 

15°N. (From Volz et al., 1981) 


Recently Singh et al. (1983) have suggested t)\at simpler and 
more easily measured molecules such as propane can be used to 
study the seasonal behavior of OH, provided the source term 
is adequately characterized. Indeed, a poorly characterized 
source function provides a major element of uncertainty when OH 
estimates are derived from mass balances of naturally occurring 
chemicals . 

In recent decades, synthetic organic chemicals have been in- 
jected into the air environment in such large amounts that measur 
able background levels are present. Because of the exclusive man 
made source of such chemicals, the uncertainties associated with 



the source term could, in principal, be eliminated. A preliminary 
screening led to the selection of methyl chloroform (Singh, 1977a, 
b; Lovelock, 1977) as the most suitable molecule for the follow- 
ing reasons: 

• Its atmospheric concentration was high enough and could be 
accurately measured 

• The usage was such that emissions were virtually identical 
to production 

• An agreed-upon OH rate constant was available 

I 

Synthetic organic chemicals thus could be monitored for their re- 
moval rate in laboratory— like conditions. 

Although a good deal of methyl chloroform data is now avail- 
able (figures 14 and 15), the earlier results (Singh 1977a, b) 
have not changed significantly: 

• The tropospheric lifetime of methyl chloroform is 8 to 10 
years, corresponding to a mean OH concentration of 6 x iq 5 
molec cm”3 (table 8) 

• The north/south gradient of methyl chloroform is too large 
to be accounted for by the distribution of its sources and 
therefore an assymetry of OH distribution must be invoked 

To best interpret methyl chloroform data, one must use a 2-D model, 
even though the basic features of OH distribution can be adequately 
obtained with the help of a multibox model. 



Figure 14. Latitudinal distribution of methyl 

chloroform. (From Singh et al., 1983) 
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Figure 15. Growth of methyl chloroform at 40°N. 
(From Singh et al., 1983) 


TABLE 8. ESTIMATION OF RESIDENCE TIMES, INTERHEMISPHERIC EXCHANGE 
RATE AND HYDROXYL RADICAL ABUNDANCE FROM COMPARISONS OF 
MEASURED AND PREDICTED RESULTS 



Measured 

Parameters - Dec. 1981 
Weighted 

Global average 


Predicted best values 
from emissions data 
and a 2-box model 

Estimated mean 
OH concentration 
(T=265 K) 


Compounds 

average 

concen- 

tration 

X (ppt) 

NH/SH con- 
centration 
ratio 

R 

Ta* 

(year! 

Tc* 

(year) 

(ppt) 

R * 
P 

X 10^^+ CH, 

Uri 

(cm^ molec"T s ^) (rtolec 

10^ 
an ^) 

Fluorocarbon 12 

294 

- 

69 

1.2 

294 

1.07+0.01 

— 

— 

Fluorocarbon 11 

179 

1.08 

55 

1.2 

179 

1.08+0.01 

— 

— 

Methylene Chloride 

29 

1.81 

0.9 

1.0 

29 

1.96+0.12 

4.27exp(-1094A) 

5.1 

1,2 Dichloroethane 

25 

2.64 

0.6 

1.2 

24 

2.63+0.25 

8.86e«p(-1094A) 

3.7 

1,1,1 Trichloroethane 136 
(Methyl Chloroform) 

1.34 

9.3 

1.2 

136 

1.21+0.02 

5.40exp(-1810A) 

5.8 

Te t r ach lor oe tha ne 

17 

5.80 

0.6 

1.2 

17 

2.71+0.27 

9.44exp(-1199A) 

5.2 


* Ta and respectively are mean atmospheric residence time and mean interhemispheric exchange rate. The 
errors in Rp are based on a value of a = 0.96 0.04. 

t Reconmended rate constants from Hampson (1980); 1,2 dichloroethane rate constant is frctn Howard and Evenson 

(1976) . The tenperature dependence is assuned to be that of itethylene chloride (R. AtJcinscn, private ccmnunicaticn) . 
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^^though methyl chloroform was the best initial choice of an- 
thropogenic traces to study atmospheric OH removal, other choices 
can now be made. Molecules such as dichloromethane , 1, 2 dichloro- 
ethane and tetrachloroethane are three other anthropogenic mole- 
cules suited for this purpose {Singh et al., 1983). A two-box 
model (table 8 ) shows that removal rates of these molecules are 
consistent with average OH concentrations of 4 to 5 x 10 ^ molec 
cm -3. Because of their relatively short lifetimes compared to 
methyl chloroform (table 8 ), these chemicals can show OH latitu- 
dinal gradients with greater sensitivity. There is no doubt that 
atmospheric removal rates of OH levels which are deduced from 
these additional chemicals can provide strong support for the 
reliability of indirect techniques outlined above. 

Uncertainties in such methods are not inherent to the tech- 
nique but are often due to differences in measurements associated 
with calibration problems. Although these differences 
have been resolved substantially over the years, further work is 
necessary. Another source of uncertainty is in the emissions 
data. ^|caus| 2 °f often complex sources of naturally occurring 
gases ( CO, CO, C 3 H 3 etc.), synthetic organic molecules offer 
a superior opportunity. To date, the synthetic molecule most 
commonly employed for OH determination has been methyl chloroform. 
Other similar molecules should be employed to gain further confi- 
dence in the predictive capabilities of the global average. Mole- 
cules such as dichloromethane , 1 , 2 dichloroethane and tetrachloro- 
ethane are strong candidates for this purpose. These somewhat 
more reactive molecules may provide additional information on the 
seasonal and latitudinal behavior of OH concentrations . Two- 
dimensional models would be ideally suited for the interpretation 
of these data. 

Radiochemical Method . The radiochemical technique developed at 
Washington State University consists of a measurement of the rate 
of carbon monoxide oxidation. This rate is an important atmo- 
spheric chemistry parameter in its own right, but current models 
indicate that it is also a direct measure of hydroxyl concen- 
tration. Species other than hydroxyl are known to oxidize carbon 
monoxide (for example, CH 3 O) , but are not expected in sufficient 
concentration in the atmosphere to cause more than 5 % of the total 
oxidation. In what follows we will tacitly neglect these other 
oxidants . 

Two rather different approaches have, for technical reasons, 
been used in a _ ground-based (Campbell et al., 1979) and airborne 
radiochemical instrument. This discussion centers on the airborne 
instrument, but a brief discussion of the ground-based instrument 
is included for completeness. 

The ground-based instrument illustrated schematically in fig- 
ure 16 maintains the normal atmospheric chemistry with as little 
perturbation as possible. A known small increment (10%) is made 
in the ambient carbon monoxide concentration, using carbon-14 mo- 
noxide, and the resulting rate of production of carbon-14 dioxide 
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Figure 16. WSU flow reactor for measurements of hydroxyl 
radical concentrations in tropospheric air. 


is measured over an interval of about 7 seconds. By assuming 
steady-state unperturbed chemistry in the instrument the hydroxyl 
concentration can be immediately deduced, given the value of the 
rate coefficient for the OH + CO reaction. 

The crux of the technique lies in ensuring that the chemistry 
is unperturbed. The sampled air is pulled into a quartz-walled 
flow reactor, where the turbulence in the flow external to the 
instrument is allowed to partially decay, and then carbon-14 mo- 
noxide of very high purity is quickly mixed into the stream from 
an array of nozzles. The flow, by then almost laminar, proceeds 
down the reactor tube for 7 seconds. The central reg'ion of the 
flow, which has not been in diffusive contact with the walls of 
the reactor, is then diverted for extraction of the carbon dioxide. 
After separation from the unreacted carbon-14 monoxide (which is 
in considerable excess), the carbon dioxide is transferred to a 
beta counting system to determine its carbon-14 content. 

Perturbations to the ambient chemistry can arise from changes 
in photolytic insolation (minimized by the use of the quartz re- 
actor, which is transparent to all photolytic radiation), from 
loss of reactive species from the flow to the walls of the reactor 
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(a subject studied extensively both by numerical flow modeling 
and by wind tunnel tests), from the increment in carbon monoxide 
(model studies indicate a somewhat less than first-order depen- 
dence), and from the effects of the beta particles emitted by the 
carbon-14 monoxide (which can be simulated by use of an external 
source of beta radiation). In practice the cumulative error due 
to all of these effects can be kept small (<20%) provided the at- 
mospheric turbulence intensity is not too high. 

A different configuration is needed for an instrument to be 
mounted on an aircraft, especially a high-speed aircraft, and 
this has required a change in the principle of operation. Super- 
ficially the airborne instrument resembles the ground-based one. 
The incoming air is decelerated to a relatively low velocity, the 
carbon-14 monoxide is mixed into it, and after a delay to permit 
reaction, the carbon dioxide is trapped out and separated for 
counting. There is, however, a major difference in approach 
caused by the impossibility of admitting sunlight to the deceler- 
ation and reaction regions of the airborne instrument; the opera- 
tion no longer assumes steady state to be obtained in the chemis- 
try. In place of the assumption of steady state we substitute a 
general objective of converting each hydroxyl radical in the col- 
lected air stream into a molecule of carbon-14 dioxide employed. 
Actual conversions range from 5% to 30%. Nevertheless the ideali- 
zation remains a useful way of viewing the instrument's operation. 

In practice the operating conditions are chosen less with the 
idea of maximizing the conversion of hydroxyl to carbon-14 dioxide 
than with that of maximizing the signal-to-noise ratio and minim- 
izing those errors whose magnitude cannot be determined indepen- 
dently. We include in these errors uncertainties in the details 
of the chemistry which occurs between the time the air enters the 
instrument and the time it leaves the reaction chamber, where 
steps are deliberately taken to quench any remaining radicals and 
bring all chemical changes to a halt. 


Figure 17 shows the major flows in the instrument and the 
methods of measuring them. The air is twice decelerated before 
it fully enters the inlet, once by passage close to the stagnation 
point on an axisymmetric airfoil mounted externally to the air- 
craft's boundary layer, and again by passage through a second 
quasi-stagnant region just inside the mouth of the inlet. At 
this second deceleration point the central part of the flow pro- 
ceeds down the inlet, while the air entering further from the 
axis is reflected by the pressure gradient and forms a reverse 
flow region immediately adjacent to the inside wall of the inlet. 


For both of these deceleration steps, which reduce the velocity 
from as much as 250 m s“^ to a nominal 27.4 m s“^, the associated 


boundary layers are very thin (because of the short development 
times). In particular, they are much thinner than the 'insulat- 
ing' air layer between the airstream actually being sampled and 
potential radical sinks, such as the metal face of the decelerat- 


ing airfoil or the metal walls of the end of the inlet tube. 


Consequently, negligible radical losses (<1%) are expected up to 
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Figure 17. WSU hydroxyl instrument MK III. Principal 
flow path schematic. 


this point, provided only that the flow not be grossly turbulent. 
Note that upper limits for these losses can be measured by the 
same techniques employed for measuring losses in the other parts 
of the apparatus (see below), and such measurements will be made 
when an opportunity offers. 

As the sample traverses the length of the inlet to the point 
of injection of the carbon-14 monoxide, a distance of 26.7 cm, the 
flow is still further decelerated by graded suction through the 
inlets 's porous wall. The final velocity at the injector is re- 
duced to 11.1 m s“^. The loss of hydroxyl by chemical decay dur- 
ing transit through this section, which requires 0.0149 s, is 
small, even if not entirely negligible. Chemical lifetimes in 
the incoming air in reasonably unpolluted regions can be ex- 
pected to be in the range from about 0.67 s at sea level (150 
ppbv carbon monoxide content) to about 3.4 s at 13 km altitude. 
The corresponding losses due to decay during transit through the 
inlet range from 2.2% to 0.4%. 


Wall losses are minimized by the suction. Air which has 
come close enough to the wall to have suffered loss of radicals 
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by diffusion and surface recombination also has a high probabil- 
ity of being removed from the flow by being sucked through the 
porous wall. Laboratory measurements have put upper bounds on 
these losses of about 4%, provided the flow is reasonably stable. 
(Early measurements on the CV-990 aircraft in August 1983 dis- 
closed much higher losses, on the order of 15%, associated with 
grossly turbulent flow in the inlet. Consultations with Pratt 
and Whitney indicated that the source of the problem was mis- 
alignment of the major axis of the external deflecting airfoil 
from the direction of the airflow prevailing in its location.) 

The admixing of the carbon-14 monoxide into the airstream is 
the next step in the process. Mixing occurs at the throat of the 
diffuser section with a diametrical expansion ratio of 1.91. The 
walls of this diffuser are also porous and subjected to graded 
suction which removes 17% of the incoming flow, again minimizing 
the losses of reactive species from the main part of the flow. 

The injector itself consists of four nozzles of 0.1 mm tubing. 
Measurements of the turbulence intensity and scale length in the 
region immediately downstream indicate prompt and thorough disper- 
sion of the carbon-14 monoxide, which is typically supplied at a 
mixing ratio of several hundred ppm in helium carrier gas. The 
exit velocity of the carrier gas from the individual nozzles 
(about 30 m s“^ ) is chosen to give a well-developed, rapidly 
broadening jet behind each nozzle, without adding significantly 
to the momentum flux in the main flow. 

Losses of radicals to the injector nozzles themselves may be 
calculated from the extensive data compiled for heat transfer 
from small cylinders in a flowing gas. These calculations show a 
maximum loss at low aircraft altitudes of about 2%, falling to 
less than 1% at 13 km. 

From the diffuser the air passes into a reactor section, 
where sufficient time is allowed for the OH + CO reaction to pro- 
ceed part way to completion. Again the walls are porous, with a 
graded suction which reduces the mean axial velocity from 2.54 m 
s“l at the inlet end to 0.51 ro s~^ at the exhaust. At the exhaust 
end a screen is interposed to quench any remaining radicals and 
more sharply define the reaction time, which has a nominal value 
of 0.218 s. Laboratory studies of the wall losses, which are 
greatest in this section because of the lowered axial velocities, 
put upper bounds of about 5% under conditions of reasonably 
stable flow. 


The subsequent operations focus on efficiently collecting all 
the carbon dioxide in the sampled air and separating from it the 


unreacted carbon-14 monoxide, which is in excess of the desired 
carbon-14 dioxide by factors ranging from 10^ to 10^. Stringent 
precautions must be taken to prevent spurious (especially surface 


catalytic) conversion of 
Separation ratios of 10^ 


carbon-14 monoxide to the labeled dioxide, 
are attained in practice by means of a 


three-stage freezing/sublimation process. A corresponding require- 


ment is that the carbon-14 monoxide injected be free of contami- 
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nant carbon-14 dioxide; proportions of less than 10“^ can now be 
consistently obtained for periods of a few seconds after final 
purification. 

In addition to the collection of the carbon dioxide, an ali- 
quot of the reactor effluent is collected in such a way as to 
represent an integrated sample over the duration of the measure- 
ment, which is typically 100 s. This aliquot is both measured in 
quantity (moles) and counted to determine the mixing ratio of the 
carbon-14 monoxide in the reactor. The flow rate out of the re- 
actor (and also out of each suction chamber) is measured, together 
with the pressure and temperature of the airstream in the reactor. 

The sensitivity of this technique (see Appendix H for a de- 
tailed discussion) is limited by two sources of noise; counting 
errors and separation inefficiencies. Counting errors arise from 
the finite background of the carbon-14 beta counting system, and 
require that sufficient activity be collected (as carbon-14 di- 
oxide) to provide at least 0.25 counts per minute. In general 
this can be achieved using practical run durations of about 2 
minutes, even for very low ambient hydroxyl concentrations . 

Separation inefficiencies in practice are the limiting fac- 
tor in the attainable sensitivity. Equations (H8) and (H9) re- 
late the ratio of carbon-14 dioxide to carbon-14 monoxide to the 
ambient hydroxyl concentration; values of this ratio less than 2 
X 10“^ may not be presently measurable. The corresponding ambient 
hydroxyl concentrations are 5 to 10 x 10'^ cm”^ . 

Long-Path UV Absorption 

Absorption spectroscopy has, in general, been a versatile 
approach for direct detection of many atmospheric species. This 
is a technique that is implemented by measuring the change in the 
radiance from a source that is due to absorption occurring as the 
radiance is directed over a known path length. The major advan- 
tages of this approach have been its simplicity of implementation 
and its versatility for detection of a wide range of species. 

The disadvantages are associated with interference from the ab- 
sorption features of other atmospheric molecules and the diffi- 
culty of extracting small signal changes embedded in the source 
radiance. In addition, for trace species, the optical path must 
be on the order of 1 to 10 km, thereby presenting difficult imple- 
mentation problems from other than ground-based platforms. 

Traditionally, the IR spectral region has been characterized 
as the fingerprint region for atmospheric measurements. This 
characterization stems from the fact that most polyatomic mole- 
cules have vibrational-rotational bands in this region, and in- 
deed the OH radical exhibits absorption features in the near IR 
as well as in the UV. However, the IR absorption has not proven 
suitable for detection of ambient OH, due primarily to the exces- 
sive interference from ambient water vapor. Implementation of the 
long-path absorption technique for measurement of OH utilizes the 
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strong, narrow, and well-resolved rotational absorption lines in 
the near UV. The intense Qi ( 2 ) line in the (0-0) band of the 
transition at 307.9951 nm has been selected as 
most suitable for the detection of OH because (a) it is one of 
the strongest lines, and (b) there is a second weaker line, the 
021^2), at 308.0006 nm in its immediate vicinity which can also 
be recorded to give a uniquely characteristic absorption spectrum 
for OH. The presence of atmospheric OH would therefore be identi- 
fied by two spectral features, which lends much higher confidence 
to its observation. An instrument utilized by Hiibler et al. (1982) 
is described schematically in figure 18. A laser system, con- 
sisting of a mode-locked argon ion laser which pumps a frequency- 
doubled dye laser, serves as the light source. The emitted UV is 
pulsed with a pulse duration of less than 50 psec and a repetition 
rate of 82 MHz with an output power of 2 mW. The spectral half- 
width of the UV output is about 0.1 nm. For comparison, the 
pressure-broadened width of the Qi ( 2 ) line is 0.002 nm. The re- 
solution of the monochromator used here is 0.002 nm. To reduce 
both the UV flux density and the divergence, the laser beam is 
expanded in a Cassegrain type telescope to a diameter of 0.3 m. 

The light path length was determined by a mirror 1.5 to 5 km away. 
The reflected light is focused on the entrance slit of a Spex 
double monochromator with a 0.05-nm-wide portion of the spectrum 
scanned repetitively by 1092 slits (10 pm wide, 10 mm high) 
placed radially around the periphery of a rotating metal disc 
(200 mm diameter). A stationary mask with a rectangular opening 
insures that light passes through only one slit at a time. The 
spectrum is scanned at a repetition rate of 6 kHz and the wave- 
length-dependent intensity distribution ' is monitored by a photo- 
multiplier. This signal is amplified and converted to digital 
signals corresponding to wavelength intervals of 0.0003 nm. The 
180 segments of each scan are stored and averaged separately in a 
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Figure 18. Principle of apparatus for 

long-path optical absorption 
using laser light. (From 
Hiibler et al . , 1982) 







fast signal averager, which is synchronized by an infrared light 
barrier. The scans are further integrated over 5-minute intervals 
and then stored. 

A calibration for the position of the OH lines in the spec- 
trum is repeated every 5 minutes. This is achieved by recording 
the strong absorption from OH radicals in an acetylene-air flame 
placed in the unexpanded laser beam. Whenever the spectral posi- 
tion of the Qi(2) line has changed more than 0.0006 nm the data 
of the previous 5 minutes are discarded to avoid deterioration of 
the spectral resolution. 

To detect an OH concentration of 10^ molecules cm”3, a mea- 
surement accuracy sufficient to resolve a differential absorption 
on the order of 10“^ is required. Because of noise introduced 
mainly by light intensity fluctuations, a noise level as low as 
10“^ has not as yet been achieved with a 5-minute measuring in- 
terval. Therefore a number of integrated scans are superimposed 
with respect to the Qi(2) reference line from the flame. In this 
way compensation for any long-term wavelength drift of the spec- 
trograph is achieved. Absorption features with an optical density 
of about 10“4 have been extracted from the integrated data via 
correction for the baseline profile. In general, the resulting 
spectrum is a superposition of spectra from several constituents 
such as SO 2 , CH 2 O, etc. Deconvolution has yielded characteristic 
OH spectra; an example is shown in figure 19 for an optical path 
of 9.7 km. The average OH concentration inferred from these data 
is 3 X 1Q6 OH molecules cm~3. 


The spectra from long-path UV absorption spectroscopy can 
provide convincing evidence of the presence of atmospheric OH, 
even after correction for the contribution from other atmospheric 
constituents. The minimum detectable OH level is a function of 
the optical path length that can be used, and this in turn is 
strongly influenced by the atmospheric conditions. 


Demonstrated 


sensitivities on the order of 1 x 10^ OH molecules cm 3 for path 
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Figure 19. Atmospheric OH UV absorption 
measurement at Deuselbach on 
Sept. 24, 1980. (From Hiibler 
et al. , 1982) 
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lengths of 3 km have been achieved. It should be noted that di- 
rect scaling with the path length is often not realizable because 
of the increased turbulent noise or transmission losses experienced 
with increases in the optical path. 

Measurement Validation 

Validation of measurements of trace atmospheric species is 
clearly an important phase in the development of any given instru- 
ment. Validation of atmospheric measurements, however, has, in 
general, an extra dimension of difficulty as compared to the nor- 
mal laboratory measurements encountered in many scientific fields. 
This new dimension is the inability to control the environmental 
parameters surrounding a given measurement. This leads to an in- 
herent lack of repeatability of the measurement conditions and 
hence a loss of one traditional method for validation (i.e., re- 
peatability). Under such circumstances certification and/or vali- 
dation of a measurement technique (and results) must rely upon 
comparison with model results, calibration standards, and/or in- 
tercomparisons of measurements made by instruments using different 
principles of detection. Of these three, the latter two are more 
desirable since the purpose of many atmospheric measurements is 
model input or model verification. 

Calibration . The use of standards in the field calibration of a 
given instrument provides a measure of the basic response and lin- 
earity of the instrument. Such a calibration generally entails 
the introduction of known amounts of the desired species into the 
sample chamber or volume used by the instrument and preferably at 
concentration levels expected to be encountered in the atmosphere. 
For OH, the establishment of a known calibration source is not a 
trivial task, particularly at ambient atmospheric levels. More- 
over, the range of sampling geometries that have been employed 
for OH instrumentation (e.g., remote vs. in situ) and the varia- 
tions with a given geometry (e.g., ambient vs. low pressure) in- 
dicate the range of difficulties that might be expected in estab- 
lishing calibration procedures for OH. It is clear that no 
universally accepted calibration procedure currently exists for 
OH, and each instrument approach has in most cases been "cali- 
brated" using procedures uniquely applicable to the given instru- 
ment. These various "calibration" approaches have been discussed 
in the previous sections. 

Notwithstanding these individual approaches, it was the con- 
sensus of the workshop participants that a more unified approach 
to OH calibration should be considered. Such an approach should 
be able to accommodate the variety of sample geometries and con- 
ditions, and in particular the lidar and in situ sampling geome- 
tries. One possible approach to a "universal" calibration station 
is being considered by the Ford group for their lidar system. 

Their approach entails the use of a number of medical UV lamps 
arranged around the circumference of a 20- to 30-cm-diameter , 10- 
meter-long tank. A multipass absorption cell configuration would 
be included in the design to provide means to measure the OH con- 


56 



centration within the cell via UV absorption. Direct calibration 
of the lidar approach appears to present the more difficult of 
the calibration tasks because of the large sampling volume. 

The possibility has been raised of modifying the design of 
this device so that other instruments could also be calibrated. 
While there appear to be no problems in calibrating these instru- 
ments with such a device, the feasibility of modifying the device 
to supply the relatively large flows required by the in situ tech- 
niques requires further study. It was agreed that the availabil- 
ity of such a device at Ford may materially assist in obtaining a 
uniform calibration of all instruments in use. 

Beyond the calibration of individual instruments, questions 
were raised about the stability of any and all calibrations. De- 
termination of the stability of the individual instruments must 
be relegated to the investigators and little information is avail- 
able at the present time. Should an instrument have problems in 
this area there are a number of steps which can be taken to rec- 
tify the lack of stability, to the extent justified by the current 
accuracy requirements. These steps include, for example, laser 
frequency stabilization, and do not add greatly to the complexity 
or cost of the instrument. 

Instrument Intercomparisons . No two hydroxyl instruments have 
ever made measurements at the same place at the same time, and so 
no intercomparison data can be presented. Notwithstanding the 
short history of hydroxyl measurements and the paucity of instru- 
ments, it may still be contended that this unusual situation 
should not persist longer than necessary. The techniques are 
exacting and often novel, so a certain scepticism of the data 
produced by any one instrument by the scientific community is 
only to be expected (and even, indeed, to be hoped for). This 
scepticism is unlikely to abate until the extent of agreement of 
measurements made by different instruments on the same samples of 
the real atmosphere is established. It is a reasonable, although 
not necessarily correct, presumption that initial instrument in- 
tercomparisons will disclose disagreements attributable to stabil- 
ity, interference or other problems. The existence of these pro- 
blems should be ascertained at an early time by instrument inter- 
comparisons, which may also serve as a convincing validation of 
the instrument calibrations. 

Any intercomparisons should be carried out in real atmo- 
spheres representative of the atmospheres to be encountered dur- 
ing global measurement programs. There is no reason to antici- 
pate that species present in heavily contaminated atmospheres 
will affect hydroxyl concentration measurements (i.e., act as 
interfering species in the broad sense of the word) , but data 
showing divergences between different types of instruments in 
unrepresentative atmospheres will not be discounted until such 
time as they are shown to occur also in more representative at- 
mospheres. Within this restriction, a reasonable variety of at- 
mospheres should be employed during instrument validation. 
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The calibration stability of instruments should not escape 
scrutiny during intercomparisons. For this reason it may be de- 
sirable to operate instruments for longer periods between recali- 
brations than would normally be the case in actual field measure- 
ment programs, in order to disclose drifts in calibration more 
clearly . 

Collection of data from intercomparisons, in the face of 
time-varying hydroxyl concentrations and differing integration 
times/sensing times from one instrument to another, will require' 
reasonably good signal-to-noise ratios in the measurements. Most 
instruments have a signal-to-noise ratio which is proportional to 
the first or higher power of the hydroxyl concentration over the 
range of concentations to be found in the atmosphere. This dic- 
tates a minimum desirable OH concentration for intercomparisons, 
such as a concentration which gives an expected standard deviation 
of no more than 15% in the ratio of the concentrations measured 
by any two instruments. Such an intercoraparison will require a 
minimum concentration, under the atmospheric conditions of the 
first section of table 6 , of 5 x 10^ cm”3 , and more if conditions 
are rendered less favorable by higher concentrations of ozone or 
NRF (presumably aerosol) species. 

Intercomparisons must be carried out between instruments 
mounted on the same aircraft. However, it also appears highly de- 
sirable to intercompare instruments on the ground, especially early 
in their development. It may furthermore prove desirable to make 
intercomparisons between an instrument on the ground (perhaps on a 
tower or a mountain top) and an instrument on an aircraft flying 
nearby. This last type of intercomparison may provide more robust 
calibration of aircraft instruments than can otherwise be obtained, 
especially for in situ instruments where inlet losses need consi- 
deration, and may provide a check on the stability of instruments 
when installed on aircraft. 

Ground intercomparisons could also profit from the inclusion 
of systems which are inherently incapable of or not ready for air- 
borne operation. This could include the long-path UV absorption 
methods and the Washington State University ground-based indirect 
OH instrument. 

Sites for any intercomparison will need to be carefully con- 
sidered from the point of view of spatial homogeneity of hydroxyl 
concentrations (including the gradient near the ground and varia- 
tions in the UV albedo of the underlying surfaces) and temporal 
variability (which will make comparisons of instruments with dif- 
fering sensitive times more difficult). The minimum concentration 
for adequate statistics must also be considered. 

MEASUREMENT OF PEROXIDES 

The major emphasis of this workshop has been on the evalua- 
tion of techniques for detection of OH. This of course is a 
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consequence of its importance in atmospheric chemistry and the 
lack of reliable measurements. The following sections, however, 
describe results of the workshop activities concerning other H^Oy 
related molecules of atmospheric interest. The deliberation 
given to these constituents is admittedly brief and by no means 
as complete as the preceding sections on OH. This should not be 
construed as a lack of importance of any of the remaining H^Oy 
species but more a comment on the development status of many of 
the techniques. 

The group was aware of three H02“related measurement tech- 
niques, of which two (Stedman, University of Denver; Cambell, 
Washington State University) measure all species which react with 
NO to produce OH, and give a value which is the equivalent HO 2 
concentration. In fact the limited rate constant data available 
show all peroxyl radicals to react with NO at similar rates, so 
these measurements would produce something close to the sum of 
concentrations of peroxyl species. 

The one technique which measures HO 2 explicitly is a frozen 
radical technique developed at KFA in Julich (Mihelcic et al. , 
1978). This is a working instrument which is relatively portable 
but has limited sensitivity (about 10^ cm“3) which may limit its 
use. 


The Stedman technique uses "radical multiplication" in the 
presence of large concentrations of NO and CO; the detected spe- 
cies is NO 2 . It has high sensitivity and is an operating, compact 
instrument. It requires calibration and its calibration stability 
has not been assessed. 

The Washington State University instrument is a modification 
of the radiochemical hydroxyl instrument. An additional injec- 
tion of NO is made, converting RO 2 to OH, which is then measured. 
This is a simple modification of any radiochemical hydroxyl instru- 
ment. Sensitivities of less than 10° cm”^ with accuracies of 30% 
are expected. 

The group was addressed by George Dawson of the University 
of Arizona, who reported upon a promising condensation wet-chemical 
technique, and Mike Rodgers of the Georgia Tech group, who reported 
upon a unique differential photofragmentation method. The NCAR 
group (Kok and Lazrus) was also reported to be making good pro- 
gress towards an improved wet-chemical method. A brief discus- 
sion of the wet-chemical and photofragmentation techniques will 
be given below. 


Wet Chemical 

Hydrogen peroxide is a highly soluble gas and therefore 
should be collected with high efficiency by condensation sampling. 
And indeed, readily detectable concentrations of H 2 O 2 have been 
found in condensate (Farmer and Dawson, 1982). Quantitation of 
the method, however, requires further work. The problems in the 
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measurement of gas phase hydrogen peroxide (H2O2) are not with 
the analytical techniques for detecting H2O2 but in the use of 
aqueous solutions for trapping gas phase H2O2. The problem with 
trapping gaseous samples of H2O2 in the aqueous solution is the 
reactions of ozone (O3) with the solution to both produce and 
destroy H2O2. For example, 

O3 + HO H2O2 + O2 

O3 + H2O2 ^ H2O + 2 O2 

The predominance of one reaction over the other depends on 
the concentrations of the reactants. Recent papers by Heikes et 
al. ( 1982 ) and Zika et al. ( 1982 ) detail these reactions more 
completely. These reactions make the accurate measurement of 
gas phase H2O2 in ambient air questionable at this point in time. 
While the first reaction has specifically been identified as pro- 
ducing H2O2 in aqueous solution, other reactions involving ambient 
air species are also possible. 

The most sensitive analytical techniques currently available 
for H2O2 determination, luminol chemiluminescence (Kok et al., 
1978 ) and peroxidase enzyme fluorescence (A. L. Lazrus, private 
communication, 1982 ), both use aqueous solutions as the medium 
for H2O2 analysis. These techniques can quantitate less than 1 
X 10“o M of H2O2 in 1 mL of aqueous solution. 

To accurately collect and analyze ambient gas phase H2O2 
using trapping techniques, methods must be found to remove the 
ozone from the air stream and/or reduce the activity of the 
collected H2O2 and H2O. This is a very difficult task since H2O2 
is very reactive and can be decomposed easily. Preliminary in- 
vestigations are under way to examine various techniques to 
accomplish this. Concurrently the development of a gas phase 
H2O2 source is also under development. 

Photofragmentation/Differential Fluorescence 

Detection of H 2 P 2 « Given that the LIF approach to monitoring hy- 
droxyl radical concentrations under atmospheric conditions proves 
successful, then LIF offers the potential for the development of 
sensitive detectors for other HxOy species. Many HxOy species 
(such as H2O2) produce OH radicals during UV photolysis; accord- 
ingly a detection scheme might consider utilizing the OH produced 
from photofragmentation of a parent species to establish its con- 
centration. As an example of this approach, we consider the de- 
tection of H2O2: 

H2O2 + hvj ->• 2 OH 
OH + hv2 ^ OH* 

OH* ->■ OH + hV3 
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since H2O2 is present in concentrations far in excess of 
naturally occurring OH, even a small fractional photolysis of H2O2 
would produce OH concentrations much greater than the naturally 
occurring OH concentrations. Consequently this approach for H2O2 
detection would suffer negligible interference from the natural 
concentration of OH. Large interferences could potentially occur, 
however, from the photolysis of other H^Oy species. For H2O2 de- 
tection the following interferences must be considered: 


(a) 

HNO3 + hv 


OH + NO2 


(b) 

H2O + hv 


OH + H (X 

< 242 ) 

(c) 

HCOOH + hv 

-V 

HCO + OH 


(d) 

CH3OOH + hv 


CH3O + OH 


(e) 

O3 + hv 

-v 

0 (^D) + 

O2 


0(1d) + M 


0 ( 3 p) + 

M 


0 (^D) + H2O 


2 OH 



Examination of the absorption cross section data along with 
typical concentration ratios indicates that there is no spectral 
region where H2O2 may be detected without interference by one or 
more of the above processes. Thus, in order to establish the 
concentration of H2O2 / some means of correcting for this interfer- 
ence must be developed. The proposed method, differential 
fluorescence, utilizes the fact that near the peak of the ozone 
absorption cross section, process (e) dominates over all other 
means of producing OH radicals. By comparing the OH produced by 
photolysis in this region (for example, = 266 nm) with that 
produced in some other wavelength region where process (e) and 
the primary process are the major processes (X2 = 222 nm) one 
can correct for the interference from process (e) by the follow- 
ing relation: 

1 F2 02 

[H2O2] = ( [OH] 2 - [OHlx) 

2oh 0 ^2 f*!®! 

2 2 


where ox 02 are the absorption cross sections for ozone 

at Xx = 266 nm and X2 = 222 nm, respectively; ^^^02 
absorption cross section for H2O2 at X2 = 222 nm; Fx and F2 are 
the laser fluences at 266 and 222 nm; and [OHlx and [OH] 2 are the 
OH concentration levels resulting from the photolysis of atmo- 
spheric air samples at 266 and 222 nm. 

The limit of our ability to determine the H2O2 concentration 
is set by the precision to which this interference can be corrected 
for. It is in this regard that the extreme sensitivity of exist- 
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ing LIF OH sensors can be utilized to the fullest extent. The OH 
concentration produced by the photolysis of these species is fully 
three orders of magnitude greater than the detection limit of cur- 
rent OH sensors. 

Detection of CH-^OOH. Unlike H2O2 described above, CH3OOH can be 
detected with little or no interference from other atmospheric 
gases. The detection of CH3OOH would begin with the primary pro- 
cess 

CH3OOH + hvx CH3O + OH 

with the subsequent detection of CH3O by; 

CH3O + hv2 ^ CH3O* 

CH3O* + M > CH3O + M 

CH3O* > CH3O + hv3 

The laser-induced fluorescence detection of CH3O has been 
reported by several investigations (Inoue et al., 1979 , 1980 ; 
Ohbayashi et al., 1977 ). The principal interferences to the 
detection of CH3O are believed to be; 

(a) CH3ONO + hv CH3O + NO 

(b) CH3ONO2 + hv CH3O + NO2 

ASSESSMENT OF TECHNIQUES 

In this section the findings of the workshop participants 
relative to the measurement capabilities of the techniques consid- 
ered will be summarized. As noted above, the largest percentage 
of the workshop's discussion was focused on techniques to measure 
OH in general and on the various LIF techniques in particular. 
Accordingly, this section will address the status of the OH tech- 
niques first and then techniques for the remaining species as a 
group. 

Table 9 reviews the OH techniques that were considered by 
the workshop participants. Under the "Techniques" heading the 
spectrum of instruments is identified by key words or acronyms 
that are intended to highlight the measurement concept and a par- 
ticularly distinguishing feature in the implementation of each 
concept. The group developing the instrument is also identified 
primarily as an additional method to characterize the particular 
technique. It is recognized that other laboratories may also be 
pursuing one or more of the concepts listed here. The comments 
that are given are directed primarily at the operational status 
of the technique at the time the workshop was convened. Foot- 
notes on these comments attempt to update this status. 
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TABLE 9. STATUS OF OH TECHNIQUES CONSIDERED AT WORKSHOP 


Techniques 

Group 

Comments 

• 

LIF/Lidar 

Ford Motor Co. 

A/C tested with 
with measurements^ 

o 

LIF/In Situ 

Georgia Tech. 

A/C tested with 
measurements® 

• 

LIF/Low Pressure 

Portland State 
University 

Ground-based 

measurements 

• 

Radiocarbon Tracer 

Washington State 

Ground-based 


(photostationary & 
flow reactor) 

University 

measurements 

• 

Long-Path UV 

Institut fiir Chemie 

Ground-based 


Absorption 

KFA, Julich 

measurements 

• 

Global Average 

SRI 

Global Average 

• 

LIF/Lidar 

NASA Goddard 

A/C test planned 

• 

Radiocarbon Tracer/ 
flow reactor 

Washington State 
University 

A/C tests planned^ 

• 

LIF/Low-Pressure 

Institut fur Chemie 
KFA, Julich 

Lab tests under way 

• 

LIF/High-rep-rate 

laser 

Harvard University 

Lab tests under way 

• 

LIF/308-nm excitation 

Jet Propulsion 
Laboratory 

Lab tests under way 


a Additional flight tests in October 1984 and April 1985 as 
part of NASA's Global Tropospheric Experiment (GTE), 
b Additional flight tests in April 1985 as part of NASA's GTE. 
c Flight test in April 1985 as part of NASA's GTE. 


The grouping of the techniques listed in table 9 is intended 
to highlight those techniques that had been developed to the point 
that tropospheric OH measurements were reported to the workshop 
(e.g., first six) or flight tests were planned (e.g., middle 
two) or were in an early development phase (e.g., last three). 

Of the six systems reporting OH measurements, only the LIF/lidar 
system developed by the Ford group and the LIF/in situ system 
developed by the Georgia Tech group had undergone some level of 
flight testing. The long-path UV absorption developed by the KFA 
group (and also the ERL Aeronomy Laboratory group) , the LIF/low- 
pressure system developed by the Portland State group (develop- 
ment by the KFA group also under way) , and the radiocarbon tracer 
systems developed by the WSU group had also reported measurements 
of OH from ground-based instruments. Of these three instruments 
only the LIF/low-pressure system appears to be compatible with 
future aircraft measurements. The global average approach repre- 
sents the last of the six concepts that reported tropospheric OH 
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measurements to the workshop. The specific question of aircraft 
compatability for this method is moot since no OH instrument as 
such is involved, although it clearly must rely upon global mea- 
surements of selected tracers from both ground and aircraft plat- 
forms. Of the remaining five techniques, two had definite plans 
for aircraft tests and the last three were deemed compatible with 
future aircraft missions. 

The majority of the techniques listed in table 9 employ some 
variation of the LIF concept. This in itself reflects a general 
consensus of the participants relative to the techniques avail- 
able for OH measurements, namely, that some variant of LIF offers 
the potential for providing sufficient sensitivity and specificity 
in an aircraft-compatible instrument to meet the needs for measur- 
ing tropospheric OH. At this point, however, the participants 
felt that it was premature to attempt to provide a more in-depth 
assessment or ranking of the OH techniques. Each of the aircraft- 
compatible techniques offered the potential for providing measure- 
ments under some conditions that would be useful to atmospheric 
models. This assessment is due in part to a lack of OH measure- 
ments in general, which tends to elevate any valid measurement 
to a status that might not otherwise be offered. Notwithstanding 
the current state of OH measurements, the various measurement 
techniques do in fact possess sufficient merit for continued con- 
sideration. As is often the case, those that are less well de- 
veloped (i.e., low-pressure sampling and 308-nm excitation coupled 
with narrow-bandwidth detection) seem to offer advantages to cir- 
cumvent the shortcomings of the more developed techniques. To 
realize those advantages, however, will require application of 
new technologies (e.g., excimer laser and high -resolution filters) 
and adaptation of instrument system to aircraft environment. 

While there appears to be no insurmountable technological obsta- 
cle, implementation of these new techniques in an aircraft envi- 
ronment will be no mean feat. 

The participants reviewing the LIF techniques recognized that 
each variant of the techniques must address interference from 
laser induced OH and nonresonant fluorescence. It was further 
agreed that techniques utilizing 282 nm excitation had a poten- 
tially larger interference problem than the 308 nm schemes. How- 
ever, it was recognized that this interference was understood 
sufficiently well to permit correction for a significant level of 
the interference. It was also recognized that such corrections 
are at best a last resort effort. A more desirable approach is 
to use an instrument configuration (e.g., lidar, low-pressure 
sampling) and/or implementation (e.g., short pulse, low pulse 
energy, 308-nm excitation, narrow-bandwidth detector) to minimize 
the noise source. The NRF background represents a noise source 
that is less well understood although it appears to be related to 
aerosol fluorescence and is assumed to be proportional to detector 
bandwidth. This results in an order of magnitude difference be- 
tween boundary layer and free troposphere noise levels. At the 
present NRF cannot be predicted and LIF systems must rely solely 
on instrument design and/or operational methods to minimize its 
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impact. Because of the apparent dependence of NRF on aerosol 
loading some of the LIF techniques listed in table 9 may ulti- 
mately be restricted to operation in a relatively clear environ- 
ment . 


An aircraft-compatible radiocarbon tracer instrument is 
scheduled for flight test as part of NASA’s GTE. The sensitivity 
of this technique will depend upon the use of materials with low 
catalytic activity. The long-path UV absorption instruments show 
no promise of being adaptable to aircraft use. Nonetheless, 
these instruments should remain as important ground-based tech- 
niques which can act as secondary calibration standards. 

The workshop discussions related to techniques for measure- 
ment of HxOy species other than OH were, at best, brief. Because 
of this, specific recommendations cannot be presented. As a 
general comment, however, it was clear that there are few if any 
instruments available for field application. There is a clear 
need for additional funding to encourage new approaches as well 
as to bring the few existing techniques to a point that field 
experience can be gained. 

CONCLUDING REMARKS AND RECOMMENDATIONS 

Before listing recommendations that have emerged from the 
workshop discussion, it is worth noting that while the workshop 
was charged with rather broad objectives, in actual fact the 
major focus of its deliberations was an assessment of a single 
measurement concept (i.e., LIF) for a single species (i.e., OH). 
This single concept currently includes a wide range of variants 
which do, in fact, exhibit fundamental differences even though 
they stem from the same basic phenomena. For the most part, the 
LIF variants that have emerged have done so because there has 
been no clear demonstration of the superiority of any one tech- 
nique or LIF variant for measuring OH. This is an important re- 
sult which emerges as the basis for the more specific recommen- 
dation noted below. Another important result that emerges from 
the workshop is the realization that the LIF concept is perhaps 
the only direct detection scheme that offers the potential for 
combining sufficient sensitivity and selectivity for tropospheric 
OH measurements from an airborne platform. This realization is 
manifested by the continued interest in developing one or more 
of the LIF variants. 

The discussions of the various LIF concepts have resulted 
in a better understanding of factors that have affected their 
successful implementation. Laser-induced OH, for example, is a 
noise source that plagued the credibility of early LIF OH measure- 
ments. It is now believed that this is one of the noise sources 
that is understood and can be predicted to an extent that correc- 
tion for a significant fraction of this interference should be 
feasible. The non-resonant fluorescence is a noise source that 
has only recently been recognized. It is one that can totally 
dominate LIF instruments and its origin and meteorology are not 
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well understood although methods for minimizing or eliminating 
its impact now appear to be in hand. A major accomplishment of 
the workshop has been the development of a common format for de- 
scribing the relationship between observed signal and hydroxyl 
concentration. Such a formalized format should provide a mecha- 
nism for a common repository of fundamental spectroscopic con- 
stants for OH, as well as a common format for evaluating current 
and future LIF systems. 

SPECIFIC RECOMMENDATIONS 

1. Alternate techniques for OH measurement should continue to 
be explored and supported if they are found to offer significant 
advantages. These include the following: 

a) LIF measurement with sample expansion to low pressures 

b) LIF measurement using higher repetition rates 

c) LIF measurement using narrower band detectors 

d) LIF measurement using 308-nm excitation in the 0-0 band 

e) LIF measurement using short laser pulse widths 

f) Measurement techniques analogous to the radiochemical 
method but using other, possibly non-radioactive, reagent 
species . 

2. A better understanding, at least at the empirical level, is 
needed of the dynamics of the OH radical, in order to assess 
the improvements to be expected in reducing interference in 
the LIF method from ozone photolysis. 

3. NASA should encourage the assembly of carefully reviewed com- 
pilations of data, such as quenching constants for OH* and 
NO 2 *, of interest to experimentalists. This effort could par- 
allel similar efforts on behalf of modellers. 

4. There should be coordinated studies, and dissemination of re- 
sults of studies, of the nature of the nonresonant fluorescent 
background to LIF measurements. Features of interest include 
its meteorology, on which some data may already be available, 
its spectrum, and its two-point autocorrelation function. 

5. Development of OH calibration techniques should be expedited, 
including exploration of common calibration systems such as 
that proposed by the Ford group, and improvements in the open- 
air long-path absorption technique. 

6. Intercomparison of OH instruments at an early time should be 
encouraged. NASA should endeavor to see that the best cali- 
bration feasible is utilized in each field expedition, and 
both the relative calibration factors and calibration stabil- 
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ities of all instruments should be determined and treated as 
part of the record of the field expedition. 

Through this report and otherwise, NASA should request the 
community to consider possible techniques for aircraft mea- 
surement of HO2 f H2O2 and CH3OOH. 
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APPENDIX A 


FUNDAMENTAL EQUATIONS OF THE LIF HYDROXYL MEASUREMENTS 


INTRODUCTION: OBJECTIVES AND ORGANIZATION 

The primary objective of this appendix is to define the equation 
relating OH concentration to the observed detector count rate for a 
given laser flux. That equation, which will be referred to henceforth 
as the "central equation" for absolute sensitivity, will be cast in 
terms of standard spectroscopic quantities such as oscillator strengths, 
Honl-London factors, rotational and vibrational quantum numbers, etc., 
such that it may be tested directly under laboratory conditions. 

A second objective will be to relate the working expressions used 
by each of the OH LIF groups to this central equation, thereby estab- 
lishing a common language for the comparison of experimental results 
and the definition of LIF system performance. 

Three sections comprise this appendix. The first details the rele- 
vant structural details of the OH radical via a spectroscopic descrip- 
tion of the energy levels and transitions employed in the LIF method. 
Second, the central equation relating observed count rate to absolute 
concentration will be derived, and third, that equation will be related 
to the observed quantities cited by each of the experimental groups 
engaged in the measurement of OH in the Earth's atmosphere. 

SPECTROSCOPY OF THE OH RADICAL 

2 2 . . 

This section describes the structure of the A E - X n transition 
of the OH molecule. The electronic spin of the hydroxyl radical is 
1/2, which produces a doublet structure. The selection rule, AJ = 0 
± 1 for the total angular momentum quantum, number, means that each 
component of the doublet has three main branches. Two vector coupling 
schemes are of importance here. The first (Herzberg, 1950), conven- 
tionally called Hund's case (a), assumes that the interaction between 
nuclear rotation and electronic motion is very weak. The electronic 
motion is coupled to the internuclear axis, as shown in figure Al. The 
electronic angular momentum, n (angular momentum will be referred to in 
units of h) , is defined as the sum of A and E , the electron orbital 
angular momentum and the electron spin, respectively. The total 
angular momentum, J, then assumes the values J=f^, J^ + 1, fi + 2,..., 
and need not be an integer. 

If, however, A = 0, as in the A state of OH, the spin vector, S, 
is not coupled to the internuclear axis and is not defined. Hund's 
case (a) then becomes meaningless, and is replaced by Hund's case (b) . 

In this case, the angular momenta. A, and N form a resultant K, which 
is shown in figure A2 . The quantum num.ber, K, is the total angular 
momentum, apart from spin, and may assume the values A, A + 1, A + 2, 
etc. It must be an integer. The two components of the doublet struc- 
ture thus correspond to J = K ± 1/2. The coupling cases are simply 
conceptual devices and represent limiting cases; for instance, as J 
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(or K) increases, a molecule that 
approach case (b) . 



Figure Al. Vector diagram 
of Hund's case (a). 


was originally close to case (a) may 



Figure A2 . Vector diagram of 
Hund ' s case (b) . 


There exist, then, six main branches: R]^ , R 2 , , Q 2 r ?!/ and P 2 , 

for AJ = ±1, 0, -1, respectively. There are also six satellite branches 
for which AK ^ AJ. These are shown in figure A3, along with the main 
branches. The superscript indicates the AK of the transition, and the 
subscript indicates the upper to lower doublet transition. Figure A3 
assumes that the state conforms to case (a) for which we can effec- 
tively divide each of the bands into two subbands: ~ 

2n3/2 - 2 e. 

The rotational term values for the x2n state have been calculated 
by Hill and Van Vleck (1928) : 

F^(J) = B^[(J + 1/2)^ - - l/2v/4(j + 1/2)^ + Y(Y - 4)X^] - 

F2(J) = B^[(J + 1/2)^ - X^ - 1 / 2/4 (J + 1/2) ^ + Y(Y - 4)X^] 

- D (J + 1)"^ 

V 

where Y A/B-^, with A e coupling constants which measure the 
strength of the coupling between the spin and the 
projection of the orbital angular momentum on the 
internuclear axis 

F^ refers to the components with J = K + 1/2 

F 2 refers to those with J = K - 1/2, as summarized here: 
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Molecular constants for OH have been reviewed recently by Huber 
and Herzberg (1979) . That review is included here in table A1 to 
provide the spectroscopic data base for this and subsequent sections. 


State 

T 

e 

CO 

e 

-1 

cm 

CO X 
e e 

-1 

cm 

B 

e 

-1 

cm 

a 

e 

-1 

cm 

D 

e 

10 ^cm ^ 

r 

e 

0 

A 

Observed 

Transitions 


89459. 

1232.9 

19.1 

9.247 

0.078 

2 

2.046 

C A 

55820.7 









C X 

(88223) 


82130 

2954 


15.2179 


16.16 

1.0809 

Di-^X 

81759.7 


69774 

660 


5.086 


9.29 

1.869 

B -> A 

35965.5 


32684.1 

3178.8 

92.91 

17.358 

0.786 

20.39 

1.0121 

A-i-'-X 

32402.3 

x^n 

0 

3737.76 

84.881 

18.910 

0.7242 

19.38 

0.96960 

1/2 ■*- 3/2 

126.23 


Summary with complete literature search can be found in Huber and 
Herzberg (1979) . 


Table Al. Spectroscopic constants for O H. 


Figure A3 (Anderson, 1971) shows that both Fn and F 2 are double. 
The rotational term values given previously do not reflect this, as 
the interaction between the rotation of the nuclei and the projection 
of X on the internuclear axis L has been neglected. This interaction 
increases with increasing rotation. The components correspond to the 
opposite symmetry with respect to inversion at the origin of coordi- 
nates and are usually designated and where i can be 1 or 2, 
depending on the doublet in question. 
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K' 



Figure A3. The transition array 
including satellite branches for 
A2z-x2n transition of the hydroxyl 
radical. 


The relative transition probabilities or line strengths for the 
^ transition have been expressed for an arbitrary coupling con- 

stant by Earls (1935) . His expressions represent specialized cases 
of formulas derived by Hill and Van Vleck (1928) . Earls' expressions 
are given by the equations below. The upper sign in the expression 
always refers to the main branch, and the lower sigh represents the 
satellite branch. 


S = 


S = 


^ 2J + i j ^ (2J +1) ± U[ (2J + 1)^ - 2A] j 


^ 2J + 1 \ /(2J + 1)^ - 2 ± U[(2J + 1)^ - 8J + 2 (A - 4) ]' 


('■ 


32 


0 , 


12 


Q 


12 


S = 


^2J + I j ^ (2J + 1) ± U[ (2J + 1) + 2(A - 4)] j 


Qt 


"12 


S = 


^ 2J + i j ^(2J + 1 


) ± U[ (2J + 1)^ + 2 (A - 4) ]' 
32 


Q. 


21 


S = 


^ 2J + i j ^(2J + 1)^ - 2 ± U[ (2J + 1) ^ - 8J -2A) ] 


32 


) 


Q, 


'Q 


21 
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/2J + l\/(2J + 1) ± U[ (2J + 1)^ - 2A) ] 
V J + 1/ V 32 



21 


where U — — 

/(2J + 1)2 + A(A - 4) 

The values of J are those of the state. 

The basic parameters in the structure of the transition discussed 
in this section will be used in the sections that follow to calculate 
the observed count rate for a given laser flux, pumping transition, 
scattering geometry, pressure, and OH concentration. 
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THE CENTRAL EQUATION FOR 
LASER- INDUCED FLUORESCENCE ABSOLUTE CALIBRATION 

In the spirit of writing an equation relating observed signal-to- 
OH concentration, we begin by assuming that the radiation field used to 
induce fluorescence is capable of pumping any combination of rovibronic 
transitions. Given the capability of modern lasers to select a single 
rovibronic transition, however, the resulting summation over pumped 
transitions will collapse to a single term in the equation defining the 
kinetics of mixed radiation and collisional energy transfer appropriate 
to this problem. 

A fully detailed treatment of this problem would require the 
solution of a coupled set of time-dependent equations of the form; 


dN (J' ,v' ,E) 
dt 


laser 


v' ,J' 


'Production rate to 
quantum level J',v' 
in electronic state' 
E' from laser 
.excitation 

^v' ,J' 


Production rate from' 
v,J energy transfer 
into v',J' from all 
other levels 


x^n 


v' ,J' 


[Removal by v,J energy) 
transfer within A^E / 
(state ) 


[Collisional deactivation ) 
to X^n state from v'J' | 
(level of A^E ) 


^rad 

I Removal by radiation 
to lower rovibronic 
state 


Examples of the functional form of these terms are indicated below 
following standard kinetic and spectroscopic notation. (See Steinfeld, 
1978, and Herzberg, 1950.) 


The state-to-state energy transfer terms into v' ,J' from all 
other levels Vj_,Jj in the A^E state will be of the form: 

P , = V k^^(v' ,J' ;v. ,J.)[M][OH(v. ,J. ,A^E)] 

V ,J Z— / 13 1 J 

i/ j 
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where the (v' , J ' ; vi , ) ' s are the bimolecular rate constants repre- 
senting the rates of collisional energy transfer and [M] is the bath 
gas. Because of the distinct difference in rates for O 2 and N 2 , a 
full solution will require two such sums, one for M = N 2 and one for 
M E O 2 . 

The rate of laser-induced excitation will take the form: 




N (v 

X 




F (v) g (v) dv 


[OH(X^n) ] 


where all quantities have been previously defined. _ We will adopt the 
convention for quantum levels such that for laser- Induced excitation 
followed by collisional redistribution and radiation, the sequence in 
quantum numbers is : 


(v ,J ,X n) 

X X 


laser-induced , 
excitation 




collision , 
transfer 


(V ,J 




radiation , 
quenching 


(v" ,J”,X^n) 


iphe two major candidates for laser— induced fluorescence in the 
atmosphere are indicated in figure A4 ; specifically. 


Option li 

2 , 282 nm , 

<''x = “'■^x'’^ "> absorbed’ <''a = ^''^A'* 

followed by 

(V = o.j.Ah) (V = 0,j",x^n) 


or 


Option 2 


30 9 nm 


(v = 0,J ,X'^n) Zll > (v. 


X 


X 


absorbed 


0,J^,A I) 


followed by 
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(V 


= 0 , J' ,A'Z) 


309 nm , 
emitted ' 


(v" = 0,J",X n) 


R-R collisional 
relaxation 


2 


V = 1 


Laser-induced 
excitation 
282 nm 


X^n 


V = 0 



Fluores- \ 
cence 
314 nm 


v-v collisional 
relaxation 



R-R collisional 
relaxation 


Fluorescence 
309 nm 


Laser-induced excitation 
Option 1 


fl uorescence 




V = 0 


Laser-induced 
excitation 
309 nm 


X^n 


0 



R-R collisional 
relaxation 


Fluorescence 
309 nm 


F2(J) 


Laser-induced excitation 
Option 2 


fluorescence 


Figure A4 . Two possible options for laser-induced 
fluorescence of hydroxyl radical. 


Given that option 2 is much simpler to treat because of the 
absence of v - v' transfer in the upper electronic state, we will 
write our "central equation" for the more complicated case represented 
by option 1. 


For tropospheric OH measurements (and in fact for stratospheric 
measurements down to a total pressure of n^lO torr) , rotational and 
vibrational relaxation rates are so fast in the time scale of fluores- 
cence that, to an excellent degree of approximation, we can represent 
the process of laser-induced fluorescence as : 


1 . 


2 

A rate of excitation from the n) level pumped by the 

laser to (v^,J^,A^E) 


2 . 


Fraction of molecules which (a) collisionally relax from 
(v^,Ja/A2z) -> (v' = 0,J',A^Z) to (b) those which leave the 
(va,Jax-^^^) level by collisional relaxation but remain in 
the A^I state, leave the (v^,Jj^,a2e) level by radiation to 
the X^n state, and those which are rovibronically quenched 
by collision 


3. Fraction of molecules in (v' = 0 ,J',a2z) which jadiate in 
competition with collisional quenching of the A^T. state 
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That is, to an adequate level of accuracy, the number of photons 
observed per pulse, Sqh/ iri a specific (v',J',a2e) -> (v",J",X^n) rovi- 
bronic transition is a product: 


^OH = 


Rate of excitation J 

from (Vj^, Jx/X^n) to / 

(vA/Jj^fA^Z) integrated) 
over a laser pulse \ 

{ Fraction of molecules 
in (v' = 0,J',A^Z) which 
radiate in competition 
with collisional 
quenching of the A state. 


{ Fraction of molecules 
collisionally coupled from 
-> (v' ,J' ,a2e) 
relative to all other removed 
terms from (v^,Jj^,a2z) 

( Fraction of the total 
number of photons emitted 
- \ within the detection vol- 
I ume which are collected 
\by detection optics 


{ Volume of scattering 
region defined by the 
overlap of the laser 
beam and the detector's 
field of view 


We proceed to write the explicit expressions for each of these terms 
in order of their appearance above: 


Step 1 : 


Rate of excitation 
fro«, (V J X2n) to , . 

xntegrated I x x 

over a laser pulse 


V V . 

X A 


2J, + 1 
A 


ue 

me 


At “ 
o 


t)g(v)dv 


[0H(X n)i 


where At is sufficiently long to encompass the full width of the 
laser pulse. 

Step 2 : 

Fraction of molecules 
collisionally coupled 
from Cva.Ja’^^^) 

Cv',J',a2e) 


V 


+ 

Q R 
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where ^ is the first-order rate constant for vibrational 

A relaxation from to v' 

is the first-order rate constant for collisional 
^ deactivation from the state 

k^ = ^ A(v' = 0,J';v" = 0,J") is the first-order rate for 

J" radiation from v' , J' to all allowed lower levels 
in the X^n state 


Step 3 : 


Fraction of molecules 
in Cv' = 0,J' ,A^Z) which 
'radiate in competition 
I with collisional quenching 
of the A state 


Y, A(v' = 0,J' ;v" = 0,J") 
allowed 
J" 


S A(V 
allowed 
J" 


= 0,J’ ;v" = Q,J") + k;! 


where the A(v' = 0,J';v" = 0,J")'s are the Einstein A coefficients 
for emission from J' to all available J"s. 

Summarizing, we have: 


, TTe^ 


At “ 
L_ o 


F(v,t)g(v)dv 


X [OH(X n)] 


V . V 

A 


. + k^ + 

^ v’ Q R 


A(v' = 0,J';v" = 0,J") 

allowed 

J" 

E A(v’ = 0,J’;v" = 0,J") + k^ 
J Lai lowed 
J" 


/ 


edV 


scattering 

volume 


[Tn] 


We proceed to explicitly define the quantities required to deter- 
mine the rate of laser excitation from the X^n state to the A^E state. 


Given the X doubling of the state and the symmetry rule for 
dipole transitions, the detailed definition of N(Vj^,Jj^,6) is important. 
We will use the convention that the normalized population refers to the 
fractional population which can absorb laser radiation at a single 
frequency. 
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Table A2 specifies the term values for the Fx(K) and F 2 components 
of the X^n state. 


K 


J 


J 

1 

0.0 cm ”^ 

1 - 1/2 

127 

1/2 

2 

84 

2 - 1/2 

188 

1 - 1/2 

3 

202 

3 - 1/2 

290 

2 - 1/2 

4 

355 

4 - 1/2 ■ 

430 

3 - 1/2 

5 

544 

5 - 1/2 

608 

4 - 1/2 

6 

768 

6 - 1/2 

824 

5 - 1/2 

7 

1028 

7 - 1/2 

1077 

6 - 1/2 

8 

1322 

8 - 1/2 

1366 

7 - 1/2 

9 

1651 

9 - 1/2 

1692 

8 - 1/2 

10 

2016 

10 - 1/2 

2053 

9 - 1/2 


Table A2 . Rotational term values for the 
doublet components of the X^n state. 


The normalized population in each of F^ (K) and F 2 (K) components is 
tabulated in table A3, with explicit account taken of the X doublet 
structure . 


Syscea 

K 

J 



'" 3/2 

1 

1-172 

.099 

.099 


2 

2-1/2 

.099 

.099 


3 

3-1/2 

.075 

.075 

F^(K) 

4 

4-1/2 

.044 

.044 


5 

5-1/2 

.021 

.021 


6 

6-1/2 

.009 

.009 

J - K + 1/2 

7 

7-1/2 

.003 

.003 

2„ 

" 1/2 

1 

1/2 

.027 

.027 

2 

1-1/2 

.040 

.040 


3 

2-1/2 

.037 

.037 

FjW 

4 

3-1/2 

.025 

.025 


5 

4-1/2 

.013 

.013 


6 

5-1/2 

.006 

.006 

J - K - 1/2 

7 

6-1/2 

.002 

.002 


Table A3. Normalized population in each rotational 
level up to K = 7 in both doublet components. 

Table A4a summarizes the literature values of the band oscillator 
strengths for the v^ - transitions of interest. At the present 
time, the only vibrational bands of interest to this subject are those 
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of -^the 0 -> 0 and 0 1 bands; for completeness, the 1 1 band is given 

for subsequent steps in this document. 


Band 

Oscillator 

Strength 

Value 

References 

^0 0 

9.5 X 10"‘* 

Huber and Herzberg (1979). 
See table Al. 

^0 1 

2.7 X lO"'* 

Huber and Herzberg (1979). 
See table Al. 

^ 1 

4.5 X lo"^ 

References in table Al. 


Table A4a. Summary of band oscillator strengths. 


Table A4b presents the dependence of the oscillator strength in the 
rotational level from the work of Wang and Killinger (1979). Results 


j " 

17(0, 0, </', J 
P Q 

r») 

R 

P 

Tj . j ; 

Q 

R 

1 

0.9145 

0.9136 

0.9118 

1.0002 

0.9982 

0.9943 

2 

0.9150 

0.9132 

0.9105 

0.9986 

0.9947 

0.9889 

3- 

0.9154 

0.9127 

0.9091 

0.9953 

0.9894 

0.9817 

4 

0.9156 

0.9120 

0.9074 

0.9902 

0.9825 

0.9728 

5 

0.9156 

0.9111 

0.9056 

0.9834 

0.9738 

0.9623 

6 

0.9155 

0.9101 

0.9036 

0.9749 

0.9634 

0.9500 

7 

0.9151 

0.9088 

0.9013 

0.9647 

0.9513 

0.9362 

8 

0.9147 

0.9074 

0.8989 

0.9528 

0.9377 

0.9208 

9 

0.9140 

0.9058 

0.8962 

0.9393 

0.9224 

0.9038 

10 

0.9131 

0.9039 

0.8933 

0.9242 

0.9056 

0.8853 

11 

0.9121 

0.9109 

0.8902 

0.9076 

0.8873 

0.8654 

12 

0.9109 

0.8997 

0.8868 

0.8894 

0.8675 

0.8441 

13 

0.9095 

0.8972 

0.8832 

0.8698 

0.8464 

0.8215 

14 

0.9079 

0.8946 

0.8794 

0.8488 

0.8239 

0.7976 

15 

0.9061 

0.8917 

0.8753 

0.8265 

0.8001 

0.7725 

16 

0.9041 

0.8886 

0.8709 

0.8028 

0.7751 

0.7463 

17 

0.9020 

0.8852 

0.8660 

0.7780 

0.7491 

0.7190 

18 

0.8996 

0.8816 

0.8613 

0.7520 

0.7219 

0.6908 

19 

0.8970 

0.8778 

0.8561 

0.7250 

0.6938 

0.6617 

20 

0.8941 

0.8737 

0.8506 

0.6970 

0.6648 

0.6319 

21 

0.8911 

0.8693 

0.8447 

0.6681 

0.6351 

0.6014 

22 

0.8878 

0.8646 

0.8386 

0.6384 

0.6046 

0.5703 

23 

0.8842 

0.8597 

0.8321 

0.6080 

0.5736 

0.5387 

24 

0.8805 

0.8545 

0.8254 

0.5770 

0.5421 

0.5068 

25 

0.8764 

0.8490 

0.8182 

0.54 5 6 

0.5102 

0.4747 


Table A4b. Values of the Franck-Condon factors q(0,0,J',J") and the 
correction factor Tji ^jn for the P, Q, and R branch transitions of 
OH. The following spectroscopic data are used; Di/Bq = 1.20 x 10~'^, 
D 2 /B 2 = 1-01 10"4, riQ =_1.012iR, r20 = 0.9706^, = 0.897 x 10 I 8 

cm , 02 - 1.056 X 10 3 cm~^. (From Wang and Killinger, 1979) 
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are suitunarized in terms of a correction factor Tji ^jn for the P, Q, and 
R branch transitions of OH. This should not be confused with the term 
values in the previous development. 


The Honl-London factors defining the rotational line strength 
have been tabulated by Dieke and Crosswhite (1961). Although more 
recent treatments are available, what we seek is a self-consistent 
treatment relating absolute OH concentration, measured^ absorption, band 
oscillator strength, and normalized line strength. This question has 
been addressed directly in the OH band oscillator strength measurement 
by chemical technique (Golden et al., 1964). 


Table A5 presents the tabulated Honl— London factors vjhich are 
correctly normalized for: 


Z S^, = 4(2J" + 1) 

j, J 



2 2 

Table A5. Rotational-transition probabilities for a S II transition 
with the coupling constant a = -7.55. (The K values are those of the 
final state.) (From Dieke and Crosswhite, 1961) 
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We turn next to the question of the overlap integral: 



(v) g (v) dv 


In order to standardize the nomenclature, we review briefly the 
relationship between peak absorption cross section, integrated absorp- 
tion cross section, integrated flux, and the above integral. Toward 
this end, we consider the case of two Gaussian envelopes summarized 
in figure A5 . The Doppler width of OH absorption, Avj-,, is somewhat 

less than the FWHM of the laser emission envelope, Av^. . 

L 



Figure A5 . Illustrative Gaussian envelopes. 

We define the Gaussian absorption line shape (see, for example, 
Killinger, 1978) : 


r r 'i -1 

a(v) = o p(v) = a Av^ exp 

O 0\ 7T / D 


-4(ln 2) 


(V - 


Av, 


noting that: 


4 In -1 


-4 (In 2) 


(V - Vj^)‘ 


Av. 


dv = 1 


Thus , 


r f 

ja(v)dv = p(v)dv = 


Note that: 


= Ffi f / ^J\J" 
% me v',v"\4(J" + 1) 
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The frequency dependence of the laser flux is given by: 

2 


F(v) = 


4 In -1 


j Av^ exp 


-4(ln 2) 


(V - v^)‘ 


Av, 


If the positions of the Doppler profile and the laser line 
width v;^ are coincident (which is typically the case) , the integral 
is straightforward: 


I 


F(v)o(v)dv = F a 


A In 2\ 1 \ 1 


O o\ TT 


exp 


-4(ln 2) 


^(v - (v - 


dv 


Av! 


Av 


L '_l 


= F a 2 


O O \ TT 


1„ L_ 

\ 2 2 
\AVl + AVj^, 


It is frequently of interest to relate this expression to the 
absorption at the peak of the absorption line: 


a(v) 


V 0\ TT 

peak 


4 In 2 


1/2 


I = o 

Av„ / peak 


So , 


1/2 

\4 In 2 / ^ D peak 


Thus , 


* 

I 

J 


F(v)o(v)dv = 


F o 

o peak 


Av 


2 . 1/2 

D \ 


2 2 
,AVl + AVj^; 


This provides a simple means of comparison given the laser 
line width, spectral character of the laser, and total flux. 
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RELATIONSHIP BETWEEN THE CENTRAL EQUATION 
AND THE WORKING EXPRESSIONS FOR EACH RESEARCH GROUP 

The objective of this section is first to link the central 
equation to each of the published working expressions used by each 
group, and second to link each of the equations used by the various 
groups to each other to eliminate ambiguities. Schematically, both 
objectives are shown in figure A6 . 



a (om 

Ford — S = A[OH] ■ n (An/n) 


Georgia Tech — D = x x x V x [ OH ] 


harvard 


— = jF(X)a(X)dA 


F(X)a(X)dA [OH] [cTn]VQ = Cq^EQ[OH] 


_ r2 ^ 

Goddard — S = Inane — ~2 


Portland State g = (P^/P ^)N GLa (OH) (n, /n) YO, 0 QT T nT 

OH c f amb p o 1 Imamgd 


Figure A6 • Schematic representation of relationship 
between the central equation and the working expres- 
sions for each research group. 
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Ford Research Group (See Wang et al. , 1981.) 

The working expression used by the Ford group is expressed in 
terms of the observed count rate: 

a (OH) 

S = a[ 0H] — n(An/n) 


instrumental parameter determined by overall 
excitation and collection efficiency 


[OH] 


hydroxyl concentration, cm 


a (OH) integrated absorption cross section for OH; 
specifically, 

a (OH) = r cf 
o o V -»■ V, 

X A 

2 2 -13 

where r^ = e /(m.c) = 2 . 8 x 10 cm 

and f is the band oscillator 

V -5- V, 

X A 

Av effective line width of the transition which takes 

into account the ratio of the laser line "envelope" 
to the Doppler width of the transition; for all the 
laser systems currently in use, the envelope of 
longitudinal modes from the doubled output of the 
dye laser is equal to or somewhat broader than the 
Doppler line width of the absorption, as noted in 
figure A5 

n fluorescence efficiency under ambient conditions 

An/n fraction of OH population in rotational level from 
which excitation originates 


We proceed next to the specifications of the correlation between 
each of the terms in the central equation and the quantities cited 
above. The identification is straightforward: 


(An/n) = N(v ,J ,0) = 

^ ^ 2 ^ 2 . 

V J 


(2J + 1) exp,[-(hc/ke)T(v ,J )J 

X XX 


(2J + 1) expI-(hc/k0)T(v ,J )] 

X XX 


^ V' ^ = 0,j';v" = 0,J") 

A 

. + E A(v’ = 0,J';v" = 0,J") + k^ 


a (OH) = TTr cf 
O O V -»■ V, 

X A 
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Georgia Tech Reserch Group (See Bradshaw and Davis, 1982.) 

The LIF signal, D, in counts/pulse, is expressed in terms of 
efficiency terms such that 

D=E, xE^xE, xE xVx [OH] 

X , f d e 



with p representing the number of laser photons per shot 

'^1 at the piimping wavelength 

a effective absorption cross section for OH at X 

Af 1 

a cross sectional area of laser beam 

^1 

f. fraction of total OH population in quantum states 
^ that can be pumped at 282 nm with a laser line 
width, A A 

S saturation parameter which is unity under con- 

^1 ditions in which the period of time required 

to deplete the level pumped is much longer than 

the repopulation time 

g /fluorescenceX 

^f \efficiency / 

with reciprocal of fluorescence lifetime 
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k, first-order rate constant for dissociation 
d 

k^ bimolecular electronic quenching rate constant 
[m] concentration of quenching species 


E 


d 


( optical \ 
detection )=y, 

efficiency/ 12 2 2 


with 



fraction of total fluorescence falling within the 
optical transmission window 

optical collection efficiency 
optical filter transmission factor 
quantum efficiency of photomultiplier tube 


E 

e 


V 


electronicX 
detection j 
efficiency / 


( volume of 
sampling 
region 


PMT signal pulses counted 
PMT signal pulses emitted 


a X 2. 
^1 


with previously defined and £ equal to the length 

of the fluorescence volume 

The relationship between the above quantities and the corresponding 
terms in the central equation is as follows. First, note that since 
the units of D in the above expression are counts/pulse, we can multi- 
ply D by the laser repetition rate Np to extract the observed count 
rate S, since 

S = DN 

P 

Considering the efficiency factors in order, we consider E^ , the 
optical pumping efficiency at X^, first. The quantity 1 
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in the exponent is the probability that a given pulse will induce a 
transition from v,J to the upper state. Since this quantity is much 
less than 1, we can expand the exponential such that 



Thus , 



With reference to the central equation, we miay equate 



The fraction of total OH population in the quantum state pumped by the 
laser, f^, is 

(2J^ + 1) exp[-Chc/k0)T(v^,J^) 

" V V r (2J + 1) expI-(hc/kBn(^ ,J )] 

V J Br ^ ^ 


and the saturation parameter is 

- (1 + Sg^,p) 

The fluorescence efficiency, E^, in the Georgia Tech expression is 
equal to 


I k- + k, + k Iin] 
\ f d q 


F 

V 

V 



I 


, I 

, I 

k , 

+ 

k„ 

+ k 

V, V 


Q 

R 


, Y, A(v' = 0,J';v" = 0,J") 

_J^ 

, + k^ + k^ E A(v' = 0,J';v" = 0,J") + kj 


Georgia Tech 


Central equation 


There is an important technical point contained in the difference 
between these two expressions because there are two distinct processes 
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2 . . . 

involved. The first is v-v transfer in the A Z state in competition 
with collisional deactivation and radiation from the vibrational level 
initially pumped by the laser. The second is fluorescence from the 
ground vibrational level in competition with collisional deactivation. 

The third term, the optical detector efficiency, is best taken 

in product with the scattering volume: 




Y, 


X Z, 


X V = 


edV 


Sea tter ing 
Volume 


[Tn] 


with a direct correspondence 


A(v' ,v")S(J",Br) 

Det Transition . 

E E ^ A(v',v")S(J",Br) 

V’ J” Br 


^2 


Z. = T 
^2 


= n 

^2 


The last term, the electronic detection efficiency, depends on the 
gating technique employed and, for tests of the absolute quantities 
discussed to this point, is best set equal to 1 by opening the detec- 
tion gate to include the entire period from laser pulse to approximate 
three folding periods of the emission. 

Harvard Research Group (See Anderson et al., 1980.) 

The observed fluorescence signal, Sqjj, in units of counts per 
second, is expressed in terms of the product: 


S 


OH 


r 

I 

J 


F(X)o(X)dX 


[OHHeTpIV Q 


The quantities on the right hand side of the above expression are 
related to the central equation in the following way: 
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F(X)o(A)dA = N(v , J , e)f 


S(J ,JJ 


X A \ Tie 


XX V v^\2J^(Rr) + 1/mc 


F(^)s('^)dv 


Harvard 


Central equation 


1 

< M 

< 

1 

A(v’ = 0,J’ ;v" = 0 , J") 

A 

J" 

^ V Q R 

L A(v' = 0,J';v" = 0, J") + kj 
. J" Q 


Harvard Central equation 

where all quantities have been previously defined and 


[eTn]V 


I 


cdV 


Scattering 

Volume 


[Tn] 


Harvard Central equation 

Goddard Research Group 


An expression for the remote laser-induced fluorescence signal 
which is valid for a lidar system similar to those described earlier 
for the measurement o£ tropospheric hydroxyl radical or nitric oxide 
is the following: 


S 


z2 


Inane A 

t 


2^ 

4tt 


T T 
12 


zl 



z 


where S 
I 


number of photons counted in a time interval (t]^,t 2 ) 
following a single laser shot and tj = 2zj/c 

number of photons emitted in the laser shot 


n molecular concentration 

a effective absorption cross section for the transition 

being pumped; this includes the cross correlation of 
the laser bandwidth with the absorption profile and 
the Boltzmann distribution of population for energy 
levels in the molecule 
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n fluorescence efficiency of the molecule; this includes 

a cross correlation between the OH emission spectrum 
and the detector bandpass 

e efficiency of the detection system; i.e. , the number 

of photons counted divided by the number of photons 
entering the system 

A area of receiver telescope; (A-(-/4Trz2) is the geometri- 

^ cal collection efficiency of the arrangement 

T, ,T„ atmospheric transmissions at the transmitted and 
^ ^ received wavelengths, respectively; in general, 

is of the form exp (-at) where a is a (hopefully 
small) constant 

0 overlap between transmitted laser beam and field of 

view of the receiver telescope; O is usually zero at 
very short ranges, increases to 1 at 0.1 - 1.0 km 
range, and remains 1 thereafter 

z range 

No term to correct for saturation has been included here, but it 
may be necessary under some operating conditions. 

The quantities on the right hand side are related to the central 
equation in the following way: 


oI = Jl 


= F(X)a(A)dA 


r 

V 

A 

v' 


I 

, I 

, I 

k 

+ k^ + 

k„ 

V. -*■ V 

A 
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E A(^ 

J" 


= 0,J';v" = 0,J") 


= 0,J’;v" = 0,J") + k 


1 / dz 

— T.T„ 0 — 

t 4-n 12 2 

z 


Portland State Group 

The observed signal, in units of counts/second, is 
Sqh = [OH]T^(Pj/P^^b)NpGL<,^(0H) (Hj^/n) YO^O^«T^T^nT^T_j 
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where T chemical transmission of inlet nozzle and tube 

c 

Pf/Pamb ^3.tio of pressure in fluorescence region to ambient 
pressure 


Np number of laser photons per second 

G White cell gain 


L 

a (OH) 
o 


n^/n 


Y 

a 

T 


m 


M 

T 


detected length per pass 

2 2 

absorption cross section for the ( n,v=0 to n,v=l) 
transition 

2 

fractional population of the detected n rotational 
level 

fluorescence yield 

overlap between absorption line and laser line 


overlap between monochromator bandpass and (0,0 + 
1,1) emission 

collector solid angle 

transmission of absorbing filter 

transmission of monochromator 


quantum efficiency of photomultiplier 
transmission of OH fluorescence by delayed gate 

efficiency of charge discriminator in detection of 
single photons 


The terms of the above working expression are related to those of 
the central equation as follows. The terms T^ and (Pf/^amb^ unique 

to the FACE sampling process and have no equivalent terms in the central 
equation. However, the meaning of the term Pf/Pamb implicit in the 
integrated fluorescence yield, Y, at lower pressure: 


N GL0„n = 
p 


F(^)g('^)dv cdV 


O (OH) = f 

° '"x '"a\2J. + 1 


S(J ,JJ\ 2 
X A \ ire 


me 


(n^/n) = N(v^,J^_,0) 
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0,J';v" = O.J") + kj 


In order to get a feel for the magnitudes involved, one may 
evaluate the expected count rate for typical values of the various 
parameters. The excitation term 


N(v ,J ,0)f , 

X X v.v 


S(J,Br) \ 
2J' (Br) + l) 


ne 

me 


F(v)g(v)dv 


J 


is essentially the product of the number of laser photons per unit area 
emitted in a shot with the probability that a hydroxyl radical will 
absorb a photon. The form ranges from 6 x 10l3 cm~2 for the Harvard 
high-repetition-rate system to perhaps 5 x 10^'^ photons pulse“^ cm“^ 
for the lidar systems. The effective cross section is approximately 
10“!^ cm.2 . Thus, the excitation term ranges from 6 x 10“5 to 5 x 10"^. 

The next term is the fluorescence efficiency of the hydroxyl 
radical. This is dependent upon pressure and temperature. Calculated 
values based upon measured values of quenching and transfer cross 
sections for N 2 , O 2 , and H 2 O are presented in figure A7. Typical 
values are 1 to 2 x 10“^ throughout the troposphere. 



Altitude, km 

Figure A7. Variation of quenching 
param.eter Q with altitude using 
constants cited in this report. 
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The next term is the product of the scattering volume and the 
geometrical collection efficiency. This ranges greatly for various 
systems. For lidar systems that have very large excitation volumes 
but low geometric collection efficiencies, this number ranges from 1 
to 20 cm^. For the contained flow systems using filter detections, 
values on the order of 0.1 to 0.6 cm^ per photomultiplier tube are 
usual. Systems using monochromators achieve 1 to 4 x 10~3 cm3. 

The [Tn] term represents the efficiency with which a photon 
entering a collection system is detected. Typical values for this 
quantity are 0.01 to 0.05. 


Finally, the OH concentration for midday midlatitude conditions 


is on the order of 5 x lo* 


.-3 


cm The low pressure sampling systems 
experience a concentration decrease on the order of 400X, but the 
fluorescence efficiency term is increased by approximately 200X, so 
the effective OH concentration could be considered to be on the order 
of 2.5 X lo6. 


Taking the product of all these terms yields numbers on the order 
of 0.001 to 1. Specific estimates of counting rate for individual 
instruments are included in the section on direct intercomparison. 
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APPENDIX B 


SATURATION PHENOMENA IN LIE MEASUREMENTS 

Optical saturation is the phenomenon in which the laser-induced 
rates of absorption and spontaneous emission between two levels, induced 
by a laser, become comparable to or greater than the spontaneous 
emission and collision rates connecting these levels. This results in 
the excited state population N^ acquiring a value of similar magnitude 
to that of the ground state, Ng. Under these conditions the observed 
fluorescence signal, which is proportional to N^/ is no longer linearly 
dependent on laser intensity I, but increases at a slower rate, and in 
principle ultimately becomes independent of I. 

Below is described a conceptual picture of optical saturation 
using a two-level picture. This is, however, inadequate for the de- 
scription of a real experiment involving a molecule, such as OH, for 
several reasons, which will be explained briefly; these are the multi- 
level nature of the electronic states and energy transfer among them 
(Kotlar et al., 1980; Crosley, 1981) and effects due to spatial (Daily, 
1978a, b) , spectral (Killinger and Wang, 1979) , and temporal (Van Dijk 
et al., 1979; Lucht et al. , 1980) fluctuations in the laser pulse. 

There exists insufficient knowledge concerning these processes and 
a lack of the needed fundamental parameters to permit one to calculate 
from first principles the effects of saturation at a given intensity. 

The same general concerns (collisional rates, collisional population 
redistribution among excited state levels, overlap of the laser and 
absorption lines) exist in somewhat different form in the linear regime. 
For those experiments in which an OH calibration cell is continually 
monitoring fluorescence from a split portion of the laser beam, oper- 
ation in the saturation regime should pose no problems with the data 
analysis. This will be true if the spectral power density beam profile 
and collisional environment within the calibration cell are made the 
same (through focussing, filtering and choice of buffer atmosphere) as 
in the test section. For experiments in which continual reference to 
such calibration cell is not made, however, optical saturation can 
complicate the data analysis, and operation in a linear regime is 
recommended. 


A TWO-LEVEL MODEL AND HOMOGENEOUS LASER 

Here, to obtain a conceptual picture of optical saturation, we 
consider a two-level treatment (fig. Bl) which ignores the multilevel 



u 


5,-<-eU 


Figure Bl 
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vibrational and rotational nature of OH, the temporal evolution of the 
populations, and temporal, spatial, and spectral fluctuations of the 
laser. 

The upper state u is pumped by the laser at a rate , and returns 
to the lower state £ by stimulated emission, spontaneous emission, and 
collisional quenching at the respective rates Bgl, A, and Q (all in 
sec"l units). The B's are Einstein absorption/stimulated emission 
coefficients related by the level degeneracies Bg = (gji/gy)Ba 
calculable knowing A. Applying a steady state approximation toN^ 
and noting that the total population Nq = is constant yields 

dN„ 

_ , 0 = - (B^I + Q + A) 


N = 
u 


B IN 

a o 

B I (1 + (g /g ) + Q + a) 
a. Jo u 


In the linear regime, where B^I << Q + A, the signal is proportion- 
al to B^I N„/(Q + A) , since << 1 and Ng - N^; (Q + A) is the 
fluorescence efficiency. As I increases to the saturation regime, the 
signal increases less rapidly than a linear dependence; as I becomes so 
large that B^I >> Q + A, then 


g N 

N 

^ ^u 

Were the situation this simple for OH, it would suggest that one should 
operate in a saturation regime not considering other factors such as 
laser-produced OH, since the signal becomes independent of I and, in 
particular, independent of collisional effects. This was the basis or 
the suggestion (Daily, 1977) and initial flurry of activity (Baronavski 
and McDonald, 1977a, b; Allen et al. , 1977; Daily and Chan, 1978; 
Boncyzle and Shirley, 1979; Pasternack et al., 1978; Allen et al., 1979, 
Van Dikj et al. , 1979) applying saturated LIE as a diagnostic tool for 
radicals in flames; of special interest was its potential use in turbu 
lent flames where single shot data must be used and the collisional 
environment thus cannot be characterized. However, a recognition of 
the quantitative problems described in the section below has demanded 
some parametric treatments and the bulk of the present effort has 
reduced to only one research group (Lucht et al., 1980). 

The two-level model nonetheless permits an approximation of the 
laser intensities for the onset of optical saturation in OH at 1 atm 
of air. We calculate I for the case in which Bgl = Q + A. At 1 atm, 

0 'V 1.5 X 1Q9 sec“l >> A. For the Qi(2) transition in the (1,0) band, 

Ba 8 X 109 cm^ erg so that B^I = Q corresponds to I - 2000 erg/cm . 
For a pulse energy of 1 mJ, a pulse length of 10 nsec, and _ a O.OIA 
bandwidth in the ultraviolet, this corresponds to a beam diameter of 
^2 mm (i.e., a power of 4 MW/cm2 at this bandwidth). 
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INADEQUACIES OF THE TWO-LEVEL MODEL 


Problems with the simple picture presented above are two- fold; 
laser beam inhomogeneities and energy transfer among the many vibra- 
tional and/or rotational levels in a real molecule. 

For a 10-nsec laser, one finds that none of the spatial, spectral, 
or temporal profiles are reproducible from shot to shot, and that, for 
example, there can be a time dependence of the mode structure (chirping) 
during a single pulse. Thus the degree of saturation will vary with 
space and time. If the laser beam characteristics were well known, 
these effects could be quantified. It has been shown (Daily, 1978a) 
that a Gaussian beam never reaches the saturation asymptote due to 
residual edge intensity remaining in a linear regime when the center is 
highly saturated; this can be accounted for if the laser beam spatial 
profile is known. Effects of mode structure for a variety of cases 
of relative mode spacing and width, and absorption line breadth, have 
also been treated (Killinger and Wang, 1978; Killinger, 1978) and can 
be calculated for defined conditions. Similar conclusions apply to a 
known temporal profile of the laser pulse (Van Dijk et al., 1979). 

None of these remains constant enough for the corrections to be 
applied on a prior basis; measurements of the variation for an indi- 
vidual pulse are virtually impossible at the present time. Fortunately, 
experience in saturation experiments (Baronavski and McDonald, 1977a, 
b; Allen et al., 1977, 1979; Boncyzle and Shirley, 1979; Pasternack 
et al. , 1978) indicate that the effects are well behaved on a pulse- 
averaged basis although empirical correction factors may be needed. 

This can be accounted for, with cancelling corrections, by using a 
continually referenced calibration cell having the same beam profile 
and average power density as in the test cell. 

The existence of rotational and/or vibrational levels in each 
electronic state can potentially lead to much greater errors. The two- 
level model would apply if Q >> R, V, denoting rotational and vibration- 
al relaxation rates (Crosley, 1981) . If however, R, V >> Q, the 
effective saturation parameter becomes Q times the upper state partition 
function (Kotlar et al. , 1980) . For OH the situation is intermediate 
(Crosley and Smith, 1980) and is complicated by the rotational level 
dependence of both the ratios Q/R in v' = 0 and Q/V for v' = 1 v' = 0 
transfer. 

Parametric approaches to these effects exist (Berg and Shackleford, 
1979; Lucht et al., 1980) but have yet to receive detailed experimental 
verification. More details are considered in a time-dependent model 
(Kotlar et al . , 1980), although the critical collisional rate constants 
(vibrational transfer with X^n^^) and the final state dependence of 
quenching collisions can only be guessed. An equipartition of energy 
scheme was used to suggest that a quenching collision from v' =0 
returns the OH to v" = 0 about half the time, and to higher v" other- 
wise. A vibrational transfer rate of ”^^8 x 10“^^ cm^ sec“^ within 
was assumed. Under these conditions the model indicated significant 
depletion of the v" = 0 population, hole burning in the rotational 
distribution within v" = 0, and a high ('\^10^) buildup of the total OH 
population into v" > 1 levels during and shortly after the laser pulse. 
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Because of the hole burning, the effective Boltzmann factor An/n at 
high I, averaged over the pulse, could be less than one-fifth the normal 
calculated value. 

The need to estimate key parameters precludes attaching quantita- 
tive significance to the results of this model, but the values used 
should be reasonable enough to draw operative conclusions. First, the 
use of a reference cell having the same collisional environment and 
laser intensity (that is, the same ratios V/Q, R/Q, and B^I/Q) in the 
test section should permit operation in a saturation regime with no 
greater errors than for the linear case. However, if a reference cell 
is not used, significant errors can result from saturation, due to the 
reduction in the effective Boltzmann factor, the calculation of whose ^ 
actual value depends crucially on unknown rates. Although the intensity 
at which saturation should become noticeable can be estimated from 
currently known rate constants, it is best at present to ascertain 
empirically (measuring the signal as a function of I) that one is not 
appreciably saturating the transition. 
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APPENDIX C 


DETAILED PERFORMANCE CALCULATIONS : 
WAYNE STATE UNIVERSITY/FORD MOTOR CO. 


The laser-generated OH through ozone dissociation was defined in 
equations (5) and ( 6 ) of the main text. Repeated here, they are: 


[OH] 


laser 


A' 


c 


F(X) [OH] 


laser 

total 


A 


c 


[OH] 


laser 

total 


(Cl) 


[OH] 


laser 

total 


k[H20] 

ao4>E[03] k[H^O] + kj^[N2] ^^[ 02 ! 


(C2) 


where F(X) = 1 
X = AAt 

A = k[H20] + k^[N 2 ] + ^^[ 03 ] 

Expressing [H 2 O] in torr, [ 03 ^] in ppbv, and E in J cm~^ , and using 
rate constants k = 2.3 x 10“^^ citw sec“l and kj^ = kQ = 3 x io~ll cm^ 
sec~l (from Streit et al. , 1976) , a value of 0.88 for the quantum yield 
of 0(^D) , and Oq = 3.45 x 10“18 cm^ (from Calvert and Pitts, 1966), 
equation (Cl) becomes 

[0H]JJ^®J = (11 ± 5) X 10® A^[H20][03]E (C3) 

where A^ is a correction factor introduced by Wang et al. (1981) to 
account for deviations from a thermalized rotational distribution for 
the laser-generated OH and for other possible processes not included 
in the rate equation analysis. Wang et al. (1981) have also determined 
experimentally that the overall numerical coefficient including A^ is 
(1.1 ± 0.6) X 10^ under ambient conditions, so that Ag = 0.1 ±0.05. 
Equation (C2) with A^, = 0.1 is at least semiquantitatively valid in 
assessing the ozone interference level in any OH measurements under 
ambient conditions. Using these equations, the ozone interference 
level corresponding to various humidity and ozone concentrations can 
be calculated readily. For the results presented in tables 6 and 7 of 
the text, the laser operated at 2.5 mJ per pulse with a beam expanded 
to 25 cm, or approximately 0.005 cm”2. 

The signal equation used for the LIF-lidar system by the Ford 
group is 


1 - exp(-X) 
X 
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n (An/n) 


(C4) 


a (OH) 

(OH signal) = A[0H]-2 

‘"eff 

where [OH] is the concentration of OH, (An/n) is the fraction of the 
OH population residing in the rotational level from which the exciting 
transition originates, aQ(OH) is the integrated absorption cross 
section for the particular rotational-electronic transition, Arigff is 
the effective line width for the transition, n is the overall fluores- 
cence efficiency under ambient conditions, and A is a constant deter- 
mined by the overall excitation and collection efficiencies. This 
constant can be evaluated most conveniently by noting that the Raman 
scattered signal owing to nitrogen molecules occurs around 3025A, for 
which the excitation configuration and collection efficiencies are the 
same for all practical purposes. This N 2 Raman scattered signal is 
given by 

(N 2 signal) = A[N 2 ]aj^(N 2 ) (C5) 

where [N 2 I is the nitrogen concentration in air, and a^(N 2 ) is the 
cross section for Raman scattering of nitrogen. Elimination of A in 
these two equations gives 


[OH] 


• T\ [N«]a,,(N„) An 
(OH signal ) ^ 2-* R^ 2' eff 

(N^ signal) a (OH) n(An/n) 

" o 


(C6) 


The appropriate values for the parameters that occur in the above 
equations are summarized in table Cl. 


TABLE Cl. VALUES FOR THE PARAMETERS OF EQUATION (C6) 


(SEA LEVEL CONDITIONS) 


Parameter 

Value 


“r'N2> 

7.4 X 10 "^^ 

cm^/molecule 

aQ(OH); Q^(2) 

0.82 X 10 “^^ 

2 -1 
cm cm 

An/n (computed value) 

0.196 


n 

0.7 X 10 "^ 


Doppler width (computed value) 

0.11 cm ^ 


Laser line width 

0.1 cm ^ 


Collision-broadened width 

0.17 cm ^ 


effective line width 

0 . 3 cm ^ 
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In the OH measurements, the exciting radiation is tuned periodi- 
cally on and off the OH resonance, and the detection electronics are 
gated for a fixed duration both during and after the laser excitation. 
For ground-based measurements at Niwot Ridge, data are recorded from 
four photon-counting channels: the Spex photomultiplier count for on- 

line and off-line intervals, called A and B respectively, and the Raman 
monitor count for the same corresponding intervals, called R;^ and Rg. 
For the OH signal levels encountered, the Spex count is essentially in 
the linear regime and no correction due to saturation need be con- 
sidered. In this case, the Spex A count is composed of three parts: 

A = al^ + bl^ + S (C7) 

where I^ represents the on-line laser intensity, al;^ the OH fluores- 
cence portion of the count, bl^ the background fluorescence, and S the 
solar background. Similarly, 

B = b'l^ + S' (C8) 

where Ig represents the off-line laser intensity, b'lg the background 
fluorescence portion of this count and S' the solar _ background. _ For 
the Raman monitor, the solar background is comparatively negligible, 
but the effects of saturation must be considered. The Raman monitor 
counts can be represented implicitly by 

"a = - V’'s> 


- V"S> <“> 


where Rg is the saturation count and cl;^ and c ' Ig are the Raman counts 
without saturation. When the spectrometer is tuned to the nitrogen 
Raman line, the A and B counts can be represented as 


A 


R 




Br = X'l 


B 


(CIO) 


The simplified form of these equations is chosen to illustrate the 
intensity dependence and presence of background counts in the data. 

The "constants" a, b, b', c, c', x, x' include, of course, many factors 
such as collection efficiencies, concentrations of species, etc. _ It is 
the ratio, a/x, of OH signal to nitrogen Raman signal which we wish to 
use in equation (C6) to obtain the OH concentration. To extract the 
constant, a, from the data, we first compute the quantity 

DEL = A/R^ - B/Rg (Cll) 


and then perform a linear least-squares regression over the measurement 
time period of the values of this quantity versus the quantity Ra/^B ~ 
1. DEL is the difference between on-line and off-line signals normal- 
ized to equal intensity of excitation. The expression (R^/Rg -_1) is 
a measure of the difference between on-line and off-line intensities. 
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-1 . 


by 


The intercept of the regression a/c(l - Ra/^S^ divided 


^r/^A = 


-1 


(C12) 


to obtain the ratio of OH signal to nitrogen Raman signal, a/x, used 
in equation (C6) . The uncertainty in the value obtained for the inter- 
cept of the regression is given by 


a 


°RES 



<V^B - 

1 


(C13) 


where n is the number of data points in a measurement time period, i 
indexes the particular data points, x represents the mean value of x, 
and aj^Es is the standard deviation of the residuals of the data from 
the computed linear fit. 


Values for a, when expressed in equivalent OH concentration, g 
ranged from 0.4 x lo^ to 1.4 x 10° OH/cm^ with an average of 0.9 x 10 
OH/cm^. The lowest value closely approached the "shot noise" limit of 
the photon statistics for the measurement. For a typical half-hour 
measurement period with 12 data points, the average value corresponds 
to a detection limit of 0.5 x 10° OH/cm° for 95 percent confidence and 
0.7 X 10^ for 99 percent confidence. 


For airborne operations, the higher levels of solar background 
from the wingtip preclude the use of photon-counting detection. Charge 
digitizing apparatus is used instead, and although it is somewhat more 
complicated, saturation of the detection system is less of a problem 
and the analysis of the data can be more straightforward. The data 
acquisition system records fluorescence signal, solar background, Raman 
monitor, UV intensity, and water vapor discharge signal for each laser 
shot. The acquisition program computes and displays results on the 
basis of a 2000 laser shot run using a straightforward analysis 
routine similar to that described above. However, since all the data 
are stored on a per-shot basis, they may be analyzed later by further 
refined analysis programs. In fact, the analysis described below is 
what we currently prefer, as it tends to give somewhat better results. 

For simplicity, we have opted to analyze our data using statis- 
tical techniques appropriated for stationary processes. Recognizing 
that the background levels and other quantitatives measured vary with 
time, we have chosen time intervals for acquisition and analysis of 
data sufficiently short that the assumption of stationary processes is 
reasonable. For situations where the conditions of measurement are 
changing rapidly, the analysis may give misleading results. 

The simplest approach uses just the on-line and off-line fluores- 
cence signals and their respective UV intensities. Subtraction of 
equation (C8) from equation (C7) gives 
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(C14) 


DIFF = al^ + hl^ 


where 


DIFF = A - B 


and we have assumed S = S' and b - b' . 

Equation (C14) describes a linear equation whose two parameters, 
a and b, can be determined through a least-squares fit according to 
standard procedures. Accordingly, the param.eter a is given by 

a = {[DIFF*I^][I^2] _ [DIFF*I^]}/A (C15) 

A = [I^2][i^2] _ (C16) 

and the corresponding variance is given by 

V(a) - (*217) 

Here the square brackets denote the sum of the products^of the quanti- 
ties therein over the entire set of observations , and s is the sample 
variance. Since the laser operates for ten shots between frequency 
switching and there are a total of 2000 laser shots per run, one may 
form from each run 100 groups, each of which contains ten pairs of on- 
off observations, and apply the above equations to each group. The ^ 
resulting 100 OH values are then processed in the usual manner to yield 
a mean value and standard deviation. (There are two "usual" manners for 
the variance, adding up the individual variances or computing it from 
the 100 OH values, and they give essentially the same results.) Note 
that the separate solar background measurement is not used in this 
analysis, nor is the fluorescence explicitly normalized to UV intensity; 
that is accomplished in effect by the two-parameter fit. 
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APPENDIX D 


DETAILED PERFORMANCE CALCULATIONS: 
GODDARD GROUP 


Two systems were described. The first used polarization to reduce 
scattered solar flux, the principal source of noise. The second used 
a narrowband detector and a different excitation wavelength. The 
results, which consider these systems, were presented in the table 
emphasizing anticipated technology improvements (table 7 in text) 
because the analysis for each system required estimates of critical 
parameter even though each system would use existing technology. For 
the first, the information about the degree of polarization of the 
scattered solar flux was not measured but was estimated as 60 percent 
for viewing straight up with a solar zenith angle of 75 percent at 
sea level. The expected performance of this system therefore depends 
on the reliability of this estimate. For the second system, two 
estimates were required. First, the expected throughput of a Fabry- 
Perot interferometer operating at 309 nm, based on the performance 
exhibited by similar instruments working at other wavelengths, was 
estimated to be 3 percent. Second, the fluorescence efficiency of the 
proposed OH lines was estimated as 4 x 10“^, since these quantities 
have not been measured to date. As before, the projections of system 
performance depend upon the accuracy of the estimates of these 
quantities . 

The laser-generated OH was calculated from the expression 


[OH]^^®®^ = 0.92 X 10® [H20][02]E 

where [H 2 O] is in torr, [O 3 ] is in ppbv, and E is in J cm“^ . For the 
marine boundary layer case [H 2 O] = 25 torr and [O 3 ] = 25 ppbv. Using 
1 mJ per laser shot with a beam diameter of about 3.5 in. resulted in 
[OH]l^s^^ = 1 X 10^ molec cm“3. For the middle troposphere case the 
expression used was 


[OH] 


laser 


3.6 X 



E 


where E is in J cm . 
resulting in [OH]^^^®^ 


For this case, 

~ 6 X 10 ^ molec 


10 mJ per laser shot was assumed, 
-3 

cm . 


The expression for the remote laser-induced fluorescence signal 
which is valid for a lidar system similar to those described earlier 
for the measurement of tropospheric hydroxyl radical is the following: 


S = In OTie A 


1 - f 

t 47T j 


2 iji ip 0 — 
^ 12 ^ ^2 
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where 


S number of photons counted in a time interval (tQ^,t 2 ) 

following a single laser shot and tj = 2zj/c 

I number of photons emitted in the laser shot 

n molecular concentration 

a effective absorption cross section for the transition 

being pumped; this includes the cross correlation of 
the laser bandwidth with the absorption profile and the 
Boltzmann distribution of population for energy levels 
in the molecule 

n fluorescence efficiency of the molecule; this includes 

a cross correlation between the OH emission spectrum 
and the detector bandpass 

e efficiency of the detection system; i.e., the number of 

photons counted divided by the number of photons 
entering the system 

A area of the receiver telescope and (A^/4irz'^) is the 

^ geometrical collection efficiency of the arrangement 

atmospheric transmissions at the transmitted and 
received wavelengths respectively; in general, Ti is of 
the form exp (-at) where a is a (hopefully small) 
constant 

0 overlap between the transmitted laser beam and the field 

of view of the receiver telescope; 0 is usually zero 
at very short ranges, increases to 1 at 0.1 to 1.0 km 
range, and remains 1 thereafter 



Z range 

No term to correct for saturation has been included here but may 
be necessary under some operating conditions. 


The correspondence between the central equation and the Goddard 
signal expression is shown by the following: 



n = [OH(X^n) ] 
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n = 


1 

X 

> 

+ 

< 

M > 

1 


r ^ 

allowed J" 

A(v' =0 , J' ;v"=0 ,J" ) ~ 

kl + kl + kl 

L Va , Vx 0 R J 

r 

E 

L allowed J" 

M 0 
+ 

rj 

0 

II 

’> 

rj 

0 

II 

> 

< 


EA 

4tt 


r^2 

J '^ 1 - 


, 0(Z) 


/ 


EdV 


scattering 

volume 


Tn 


In order to get a feel for the magnitudes involved, one may 
evaluate S given some "reasonable" values of the various parameters for 
an OH measurement: 


I 

'^'lO per 

n 

7 - 7 

'^lO ' cm 

a 

"v^lO cm 

n 

'x-0.01 

E 

1 — 1 

0 

• 

0 

> 

ft 

2 

'^'500 cm 

Ti,T2 

^1 


'\'10 per millijoule of laser energy 


can vary a great deal with various choices of 
system parameters; 9 x 10“4 is a realizable value 
for a 10-mrad field of view and Z 2 = 150 m (corres- 
ponding to a time interval of 1 ysec) 


S = 10^^ X lo”^ X 10 X 10 ^ X 10 “^ X (5 X 10 ^) X (9 X 10 = 45 

This says nothing about the ability to measure OH because we know 
nothing about interferences, but it does say that one will not be signal 
limited in the absence of background. 

The principal source of noise for the lidar measurements is 
statistical noise from the large background count arising from scattered 
solar flux. So-called white fluorescence has been of a much lower 
magnitude because of the absence of any physical objects in the 
scattering volume. All the other systems do have some sort of flow 


/ 


^20 (Z) dZ 
■ „2 


Then 
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tube or beam stop which provides a potential source for contaminants 
as well as a scattering (and possibly fluorescing) surface. 

The procedure consists of four measurements made in almost simul- 
taneous pairs (time between measurements < 1 msec) . The system is 
tuned onto the hydroxyl line with the laser, and a signal and back- 
ground are measured. Then it is tuned off the OH line and a signal and 
background are measured. The measured quantities are Sqn/ 
and B^pp,- and Sqff broken down as shown: 

^ON ^ON FLUORESCENCE ^WHITE FLUORESCENCE ^BACKGROUND 
^OFF ^ ^WHITE FLUORESCENCE ^BACKGROUND 


Assuming Bqn is the same as the in Sqn and BQpp is the same 

for SqfF/ 


F + W 


^ON 


®0N 


^ ^OFF ®OFF 


Scaling F and W for laser power, we can solve for F: 

I, 


= <=ON - BoN> - I 


ON 


OFF 


^^OFF ®OFF^ 


The uncertainty of F, op, can be calculated using the statistical 
uncertainties of the four measured quantities (assuming Iqn Iqff 

are very well determined) . Then 

"on )' 

^OFF/ 

In practice, and crsnFF largest terms in the expres- 
sion, and (Ion/Ioff^ Then 


= o. 


+ o 


’on 


B 


2 + 


ON 


2 + 0 


OFF 


B 


OFF 


= /s 


ON 


ON 


and 


= /S 


OFF 


OFF 


These values are dominated by solar background photons and are there- 
fore almost equal. 


113 



Solar background fluxes have been estimated by scaling results 
from the Goddard balloon lidar flight of October 1980. These values 
vary depending on parameters such as viewing direction, solar zenith 
angle, telescope field of view, and overall detection efficiency. 

These estimates are included in our calculated signal to noise ratios, 
tabulated in the text. 
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APPENDIX E 


DETAILED PERFORMANCE CALCULATIONS 
GEORGIA TECH GROUP 

Operational Characteristics for Two Tropospheric Environments 

(1) CONDITIONS FOR SYSTEM INTERCOMPARISON 
Case A Marine Boundary Layer 

1 atm total pressure 
25 ppbv O3 
25 torr H2O 

5 X 10® ambient OH/cm^ 

Total background of OH equivalents: 

(i) lx lo’^ OH equivalents 

(ii) 1 X 10® OH equivalents 

probed OH produced via O3/H2O interference 

given by [OH] = 5.75 x 10^ E F 

Fio ns 

where E is in mJ 
Case B Mid Troposphere 

0.5 atm total pressure 
60 ppbv O3 
0.3 torr H2O 

probed OH produced via O3/H2O 

Interference given by: 

[OH] = 3.6 X 10® E F 

Fio ns 

fluorescence efficiency = 2 x 10“® 

(2) SIGNAL STRENGTH EVALUATION 

S (OH) = E\^ X Ef X E<3 X E 0 X V X [OH] 
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= optical pumping efficiency 


where 

Ef = fluorescence efficiency 

“ optical detection efficiency 
Eg = electronic efficiency 
V = effective sampling volume 
For current field system consisting of 

1. 85 J/pulse at 500 ps 1 cm^ beam 

2. two lasers at 20 pps each 
(one on line and one off line) 

3 . Eg = 1 .0 

4. V = 1.5 cm^/PMT 

5. 12 PMT's 

6. fluorescence sampled within a 50 A bandwidth 
filter = 0.75 

7. 0.3 quantum yield of PMT 

8. optical collection and transmission 
efficiency = 3 x 10 “^ 

F = 0.1 

^10 ns 

Case A Signal Strength 

S (OH) = 1.0 X 10^ [OH] per laser pulse 
at [OH] = 5 X 106/cm3 
then S(OH) = 600 photons/min 
Case B Signal Strength 

S(OH) = 2.0 X 10^ [OH] per laser pulse 
for 5 X 10^ OH/cm^ then 

S(OH) = 1200 photons/min 



(3) BACKGROUND (B) TOTAL SOURCE STRENGTH 

B = 1 X 10^ OH equivalents 
and B = 1 X lO^ OH equivalents 
Case A Background B 

B = 12000 photons/min 
and B = 1200 photons/min 
Case B Background B 

B = 24000 photons/min 
and B = 2400 photons/min 

(4) SIGNAL-TO-NOISE RATIO (S/N) 

S/N = S(OH) 

/ S(OH) + 2B 

(5) PERFORMANCE OF PRESENT GA TECH SYSTEM 

Equivalent Integration [OH] % [OH] S/N 

Background Time Interference Interference D=1 D-5 D-20, 

Case A - MARINE BOUNDARY LAYER [OH] = 5 x 10^ cm“3 


1 X 108 

20 min 

5 X 

105 

10% 

17. 

NA 

NA 

1 X 108 

1 min 

5 X 

105 

10% 

3.8 

NA 

NA 

1 X 107 

20 min 

5 X 

io5 

10% 

49. 

NA 

NA 

1 X 107 

1 min 

5 X 

105 

10% 

11. 

NA 

NA 

Case B - 

MID-TROPOSPHERE 

[OH] = 

= 5 X 

lO^/cm^ 




1 X 108 

20 min 

3 X 

10^ 

0.6% 

24. 

NA 

NA 

1 X 108 

1 min 

3 X 

io4 

0.6% 

5.3 

NA 

NA 

1 X 107 

20 min 

3 X 

10^ 

0.6% 

69 . 

NA 

NA 

1 X 107 

1 min 

3 X 

io4 

0.6% 

16. 

NA 

NA 


NA - Integration time de-rating factor is given as D=1 for 
GA Tech system using <100 ys time-separated lasers 
(of/off line) which will "freeze" the atmosphere under 
all conditions. 
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(6) SPECIAL COMMENTS ON CURRENT GA TECH SYSTEM 


Case A - For the background levels given, the irinimum detec- 
table concentration (MDC) for a S/N = 2 is seen to 
range from 2.6 x lO^/cm^ to 8.3 x 10^/cm^. However, 
based upon information from GAMETAG flights (as to 
the background levels) we feel a more conservative 
MDC for a S/N = 2 ranges from 8 x lO^/cm^ to 1.5 
X 106/cm3. These values are all given for 20 minute 
integrations . 

Case B - For the background levels given, the minimum detec- 
table concentration (MDC) for a S/N = 2 is seen to 
range from 1.8 x 105/cm3 to 5.8 x 105/cm3 . However, 
based again on information from the GAMETAG flights 
(as to the background levels) we feel a more 
conservative MDC for a S/N = 2 ranges from 3 x lO^/ 
cm3 to 6 X 103/cm3. These values are all given for 
20 minute integration times. 

(7) POTENTIAL IMPROVEMENT IN ANTICIPATED GA TECH SYSTEM 

1. 85 yJ/pulse 300 ps 1 cm^ 

2. Two lasers at 100 pps each 

3. All other factors the same as section (2) 

4. F/Fio ns = 0.05 

Case A 

S(OH) = 3000 photons/min 

B = 6000 photons/min and 60000 photons/min 

Case B 


S(OH) = 6000 photons/min 

B = 12000 photons/min and 120000 photons/min 


Equivalent 

Background 

I ntegrat ion 
Time 

[OH] 

Interference 

% [OH] 

Interference 

D=1 

S/N 

D=5 

D=1 

Case A - MARINE BOUNDARY 

LAYER 

[OH] 5 

X 10^ cm" 3 




1 X 108 

20 min 

2.5 

X 

105 

5% 

38. 

NA 

NA 

I—* 

X 

o 

00 

1 min 

2.5 

X 

105 

5% 

8.5 

NA 

NA 

1 X 107 

20 min 

2.5 

X 

105 

5% 

110. 

NA 

NA 

X 

1— ' 
o 

1 min 

2.5 

X 

105 

5% 

25. 

NA 

NA 
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Equivalent 

Background 

Integration 

Time 

[OH] 

Interference 

% [OH] 

Interference 

D=1 

S/N 

D=5 

o 

rH 

II 

Q 

Case B - MID 

-TROPOSPHERE 

[OH] = 5 X 10 

6 cm“3 




1 X 108 

20 min 

1.5 X 104 

0.3% 

54. 

NA 

NA 

1 X 108 

1 min 

1.5 X 104 

0.3% 

12. 

NA 

NA 

1 X 107 

20 min 

1.5 X 104 

0.3% 

154. 

NA 

NA 

1 X 107 

1 min 

1.5 X 104 

0.3% 

36. 

NA 

NA 


NA - integration time de-rating factor given as D-1 for GA Tech 
system using <100 ps on/off line cycling time. 

Case A - As noted earlier we believe that the workshop agreed- 
upon noise level is too optimistic for an in situ LIF 
system; thus we expect the MDC (S/N =2) to be in the 
range of 3.5 x lO^/cm^ to 7 x 105/cm3 for a 20 minute 
integration time. 

Case B - See Case A comment; MDC (S/N = 2) expected to be in the 
range of 1.4 x lO^/cm^ to 3 x lO^/cm^ for a 20 minute 
integration time. 
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APPENDIX F 


DETAILED PERFORMANCE CALCULATIONS: 
PORTLAND STATE GROUP 


For the Portland State University FACE instrument, laser- 
induced OH is calculated from 

laser = ^ 1 [^3 ^ amb [^ 2 ^^ amb 

where [OH] laser units of signal— equivalent ambient number 

density, and [ 03 ] 3 j„j 3 and [H^Oj^j^b are also ambient densities, all 
in molec cm“ . The coefficient C^^ is proportional to the square 
of the laser flux and also to the expression 


k 4 ( 1-exp (-k 2 [M] At ) ) 


k 2 (l-exp(-k 4 [M]At) ) 


F = Y [M] I 1 - 


k 2 [M] At 


k 4 [M] At 


k4 - k2 


which is strongly dependent upon pulsewidth and pressure. [M] is 
the total number density, k 2 is the total 0(^D) removal rate con- 
stant, and k 4 _^g th| groundstate thermalization rate constant. 

^1 = V.5 X 10 cm molec has been measured experimentally wi 
the FACE instrument sampling ambient air enriched with ozone. 


th 


The FACE signal function is 


^OH ^caltO^lamb 

where Ccal is measured in a large transparent chamber in which 
FACE measures Sq^ while [OH] 3^,5 is simultaneously measured by hy- 
drocarbon decay using gas chromatography with flame-ionization 
detection. Cq^I varies with excitation overlap and with alignment 
and cleanliness of the White-cell and detection optics; at [OH] = 

5 X 10° cm we have observed = 7 counts/min. 

The dominant noise sources in FACE are the photon-counting 
fluctuations associated with the OH signal and the fluorescence 
background as well as concentration fluctuations in OH and in the 
species causing the background. Baric filtering and temporal fil- 
tering suppress the background by factors of 10 each. After sup- 
pression, the remaining fluorescence background in urban air is 
in the 10' ambient OH cm“3 equivalent range. 

In tables 6 and 7 of the text, the signal-to-noise ratio is 
calculated as SNR=S/D S + 2B. Here S is the number of counts of 
OH fluorescence observed at [OH] = 5 x 10° cm“^ with the cali- 
brated response Ccal “ photons/min in the stated integration 
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times. B is the stated atmospheric-pressure background, reduced 
by a factor of 100 via baric and temporal filtering. The degrada- 
tion factor D = 1 is for photon-noise-limited condition. D = 5 
and D = 20 allow for possible large concentration fluctuations 
between the signal and background measurements. 

In table 7 of the text, Portland State assumes a reduction in 
laser linewidth, improving the overlap with the pumped transition 
by a factor of 2.5. 
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APPENDIX G 


DETAILED PERFORMANCE CALCULATIONS: 


HARVARD GROUP 

The technique proposed by the Harvard group for the measurement 
of OH density in the troposphere has many elements in common with a 
system that is currently used for in situ measurement of OH in the 
stratosphere. The proposed experiment utilized laser-induced fluores- 
cence detection of OH in the A-X band. The sampling configuration is 
the enclosed flow technique, whereby the air samples flow at high 
velocity through a 5-in.-i.d. flow tube that houses the detection optics. 

The theory of LIF detection of ambient OH is well documented else- 
where in this publication. When applied to the proposed technique, a 
signal equation of the following form is obtained: 

Sqjj (cts/sec) = 

J*F (photons/cm^sec) a (cm^) dl [oh] (molec/cm^) V(cm^) Q eTn (Gl) 

where, in addition to the self-evident terms, we have: 

V sample volume element defined by overlap of projected area of 
laser beam and projected area of detector field of view 

2 t 

Q fluorescence efficiency; that is, fraction of A I state OH 
radicals that fluoresce rather than collisionally deactivate 

e geometric collection efficiency of detector 

T optical transmission of detector optics 

n efficiency of converting a photon incident on PMT photocathode 
into an observed count 


In our working expression the nonvariable terms in equation (Gl) 


are collected into a single product given by C 


OH' 


^OH ^OH 




cts/sec 


Q • molec/cm3 


E (mW) Q[OH] 


(G2) 


mW ^ 


where 


'OH 


= / 


Ulaser ^ 


and g^ (X) represents the laser line width function. 

Deriving values of [OH] from the observed count rate due to OH 
fluorescence requires knowledge of Cqh/ E, and Q. Cqh is determined 
in the laboratory by direct measurement of Sqh with a known OH density 
in the sample volume element at a known circulating power, E, and a 
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known quenching environment, Q. E is measured by monitoring Rayleigh 
scattering at the laser pumping wavelength, and Q is derived from the 
local conditions from pressure and temperature measurements. 

During the calibration phase, E is recorded as the ratio of 
to [M] , the buffer gas density. The signal equation for the SR^y 
is analogous to the Sqh signal equation, the only difference being the 
nature of the scattering events observed, and we have: 

^Ray (cts/sec) = E[M] (G3) 

The value of C^ay be determined easily by experiment; however, 

its value is not used explicitly in data reduction since it cancels 
out in the final form of the working expression: 


'OH 


= C 


^Ray/[M] 

calibration 


Q[OH] 


(G4) 


The value of C^ay does, of course, provide useful diagnostic infor- 
mation regarding the actual number of passes achieved by the light 
tapping cell when it is coupled with an independent measure of the power 
in the input beam. The latter is provided by photodiodes which are 
used to monitor laser performance. 

The current stratospheric instrument has the following directly 
measured sensitivity: 


C 


OH 




5 X 10 


cts/sec 
Q • LohJ 


at 


Ray 


/[M] 


10 


-15 cts/sec 
molec/cm' 


and an input beam energy of 10 mW. 


The proposed tropospheric instrument with narrowed laser perfor- 
mance and increased power has a projected sensitivity of: 


C 


OH 


5 X 10 


„ 0.4i „ 20 mW 

0. 02 A 10 mW 


2 X 10 


cts/sec 
Q • LOHJ 


(G5) 


This is the sensitivity used for the predicted performance levels of 
our anticipated system in table 8 of the text. For the results given 
in table 8 an average power of 3 mW with Q = 8.7 x 10“'* is assumed for 
the marine boundary layer cases, and an average power of 20 mW with 
Q = 1.5 X 10“^ is assumed for the middle troposphere. The current 
detector mode with (1,0) pumping at 2820& and wideband (0,0) observa- 
tion centered around 3090R is assumed for these estimates. 


123 



The strengths of the various background signals present on the 
OH and Rayleigh detectors are relatively low. The estimates for the 
anticipated systems in table G1 are based upon the performance of the 
current stratospheric system. 

The ambient OH fluorescence signal is separated from the background 
signals by tuning the laser off the OH absorption line and directly 
measuring the background signal. This approach is sufficient for the 
majority of background sources, which are nonresonant. The notable 
exception is, of course, the laser-induced OH term. The magnitude of 
that source and its measurement are discussed below. 


The projected laser-induced [OH] values in table G1 are based on 
the following equation (see section in text on photolytic interference) : 


[03]tH,0] 

tO“haser = 3.8 . f (x) 


where we start with the nominal operating parameters : 


(G 6 ) 


E 2 0 mW 

N 17,000 

P 

2 2 2 

A Trr = IT (0.5 cm) = 0.78 cm 

f(x) 0.56 at 1 atm T = 296 K 

0.46 at 0.5 atm T = 256 K 

Equation (G 6 ) predicts the density of laser-induced OH that is 
formed and detected during a single pulse in a laser beam of peak flux 
equal to E/NpA. However, the light trapping cell multiply reflects the 
laser beam through the sample volume element at least 40 times. The 
net flux through the volume element is increased since the 't^I cm 
diameter laser beams overlap spatially. The net increase in flux is 
about a factor of 10 , depending on the actual geometry of the light 
trapping pattern we choose. The values of laser-induced OH reported 
in table G1 take this effect into account. 

The minimum flow velocity to clear the sample volume element 
between pulses at Np = 17,000 is 

V = 17,000 X 1 cm = 17,000 cm/sec = 380 mph 

The actual contribution from laser-induced OH to the observed 
signal must be measured directly when it is present. If it is negli- 
gible, it must be demonstrated as such. 

Fortunately, experimental susceptibility to the resonant inter- 
ference can be readily mapped by performing nighttime measurements , 
when the ambient OH concentration is extremely low. This strategy 
would require simultaneous H 2 O and O 3 density measurements during 
both day and night. 
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Alternatively, the enclosed flow technique allows for a chemical 
null experiment specifically tailored to separate the ambient OH siq~ 
nal from the laser— induced OH signal. A flow of C 3 Hg can be mixed into 
the sample flow at a concentration suitable to reduce ambient OH density 
a factor of ''^100 on a millisecond time scale. The laser— induced OH 
signal is unaffected since it is operative over a nanosecond time scale. 
The efficiency and presence of any possible side effects can be checked 
in the laboratory and by a suitable strategy of complementary daytime 
and nighttime measurements. The comparison between the chemical 
addition null experiment and the conventional null experiment, tuning 
the laser off the absorption line, provides another diagnostic check 
of system performance. 


TABLE G1 


Background 

Source 

Current 

Stratospheric 

System 

Proposed 

Tropospheric 

System 

Comments 

S3090 

nonresonant 

background 

fluorescence 

<5 


No prediction made. 
Presence of many more 
compounds in the tropo- 
sphere makes any predic- 
tion unreliable. 

S3090 

chamber 

scatter 

<5 

<5 


s3090 

Raman 

Rayleigh 

<5 

<5 


S3090 

solar scatter 

70 

<5 

Improved baffling and a 
narrowed PMT gas width 
will decrease this 

S3090 

dark count 

<1 

<1 

source. 

S3090 
laser- 
induced OH 

<5 


See text. 

S2820 

chamber 

scatter 

<1 

<1 


S2820 

dark count 

<1 

<1 


S2820 

solar scatter 

<1 

<1 
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APPENDIX H 


REDUCTION OF DATA FOR RADIOCHEMICAL METHOD: 

WASHINGTON STATE UNIVERSITY GROUP 

We postulate the following very general model for the 
chemistry of the air before and after the addition of the carbon- 
14 monoxide. Beginning with the ambient, which we assume to be 
in steady state, we write 



where n is the OH concentration as a function of time, t 

Po is the production rate of OH radicals, assumed in- 
dependent of n 

Tq is the lifetime of the radicals against chemical loss, 
a process assumed to be only first-order in n 

The production term includes several components, some of 
which will vanish promptly when the air is transferred out of the 
light, for example into the instrument inlet. We write 

P = P + P (H2) 

o p 

where Pp is the photolytic (promptly light-dependent) part, and P 
IS the residue. Furthermore, let us write 

P = f P = fn/To (H3) 

O 

where f is the non— photolytic fraction of the source term. The 
value of f can be estimated from models and can be measured, but 
will often not be known with any great accuracy. 

After the admixing of the carbon-14 monoxide the kinetic 
equation (HI) will be changed by the addition of another loss term, 
"go = “o/t', where k is the rate coefficient for the OH + 

CO reaction under the prevailing conditions, n^Q is the concen- 
hi^3tion of carbon— 14 monoxide and t' is the lifetime of hydroxyl 
against reaction specifically with carbon-14 monoxide. If the 
initial (injection point) concentration of hydroxyl is n°, we may 
solve to find n(t) after the carbon-14 monoxide is added: 

n(t) = (n° - Pt) e“^/T + Px (H4) 

where x, the overall lifetime of hydroxyl, is given by 

1/x = 1/xo + 1/x' (H5) 
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(H6) 


This equation can be rewritten as 

n(t) = n°[(l - fx/x^) + fT/T^] 

Given n(t) we can calculate the concentration ^c02 ^*^2 arising 

from the reaction 

+ OH ^'^C02 + H 

which we assume to be the dominant source. Thus 

npp, - n° = [(1 - tx/x) (1 - e"'^/^) + fT/x ] (H7) 

2 CO^ ^ 

where n^Q is the initial ^^C02 concentration and T is the 

duration of the reaction from mixing to quenching. A useful 
parameter is the normalized sensitivity S*, 

T ' 

S* = — - (nco - n°co ) (H8) 

Tn° 2 2 

or 

S* = - [(1 - fx/x ) (1 - e"’^^'^) + fT/r^] (H9) 

tji o o 

Values of S* can be calculated from equation (H9), and tend 
to cluster just below unity. Similarly, the strength of the de- 
pendence of S* on parameters, such as the reaction time T, can 
be assessed by computing 91n S*/31n T, and so on. Typical 
values are shown in tables Hi and H2 for representative sea-level 
and 8-km-altitude conditions respectively. It should be noted 
that S* shows a generally weak dependence on the operating condi- 
tions, indicating that the measurement results are not critically 
dependent on measurement of any parameters other than the activity 
of the carbon-14 dioxide, the mixing ratio of carbon-14 monoxide, 
the reaction time and the flow rate. In particular, the depen- 
dence upon f is very weak, with typical values of 91n S*/9f 
ranging downward from a maximum of 0.14, an important point inso- 
far as f is poorly known. Total ignorance of f, combined with an 
assumption that f = 0.5, cannot produce an error in n° exceeding 
7%. 

Reduction of the data begins with the assay, the time-inte- 
grated aliquot of the carbon-14 monoxide admixed flow. From the 
volume of the assay container, its increase in pressure during 
the measurement, and the laboratory temperature, the number of 
moles of air in the container are calculated. The contents are 
then transferred to a proportional beta-counter and the total 
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carbon-14 content (corrected for counting efficiency) is obtained. 
From these data the concentration n^o of carbon-14 monoxide in the 
reactor can be computed, using the reactor temperature and pres- 
sure . 


From the net carbon-14 dioxide counts, again corrected for 
counter efficiency, the reactor effluent mass flow rate, the run 
duration and the reactor temperature and pressure we compute the 
net increase in carbon-14 dioxide concentration at the reactor 
exhaust, - n^Q 2 ) . The ratio of this quantity to n^Q, deter- 

mined from the assay, permits calculation of the hydroxyl concen- 
tration at the injector. 


n° = (nco - )/(S* k nco T) (HIO) 

2 2 

apart from a small correction for wall losses in the diffuser- 
reactor unit. 

The external hydroxyl concentration requires additional 
corrections for losses of radicals to the injector nozzles (1 to 
2%), to the walls of the inlet tube (believed to be about 4%) and 
due to chemical decay during the transit time (less than about 
2%). Overall errors of less than 20% appear to be currently at- 
tainable (provided the problem of turbulence due to misalignment 
of the external airfoil is resolved). With some improvement in 
flow and temperature monitoring the overall errors could be 
reduced to less than 15%. 
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TABLE HI FRACTIONAL ERROR COEFFICIENTS AT SEA LEVEL IN CLEAN AIR 


Coeff . 

Quantity f = 

0.0 

0.2 

0.4 

0.6 

0.8 

1.0 

9 In S* 
9 In T 

Reaction time 

-0.26 

-0.23 

-0.20 

-0.17 

-0.15 

-0.13 

9 In S* 
9 In Tq 

OH lifetime 

+0.11 

+0.09 

+0.06 

+0.03 

+0.01 

-0.01 

9 In S* 
9 In 

I'^CO concentration 

+0.14 

+0.14 

+0.14 

+0.14 

+0.04 

+0.04 

9 In S* 
9f 

Non-photoly t ic 
fraction 

+0.14 

+0.13 

+0.13 

+0.13 

+0.12 

+0.12 


TABLE H2 FRACTIONAL ERROR COEFFICIENTS AT 8 km ALTITUDE IN CLEAN AIR 


Coeff. 

Quantity f = 

0.0 




0.8 

1.0 

9 In S* 

Reaction time 

-0.14 

-0.13 

-0.12 

-0.12 

-0.11 

-0.10 

9 In T 








9 In S* 

OH lifetime 

-0.04 

+0.03 

+0.02 

+0.01 

+0.01 

0.00 

9 In Tq 








9 In S* 

I'^CO concentration 

+0.10 

+0.10 

+0.10 

+0.10 

+0.10 

+0.10 

9 In 








9 In S* 

Non-photoly t ic 

+0.05 

+0.05 

+0.05 

+0.05 

+0.05 

+0.05 

9f 

fraction 
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